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FOREWARD 


This  Is  the  final  report  on  research  In  Radiation/Catalytic 
Augmented  Combustion  conducted  with  partial  support  of  the  Air  Force  Office  of 
Scientific  Research  under  Contract  No.  F49620-81-C-0028.  The  research  under 
this  contract  continues  the  work  performed  under  previous  Contract  No.  F49260- 
77-C-0085  and  described  by  Final  Technical  Report  (AF0SR-TR-82-0132) . 

This  final  report  covers  the  progress  of  research  since  its 
inception  on  April  l,  1981  through  March  31,  1983.  The  work  was  performed 
under  the  direction  of  Exxon  Research  and  Engineering  Company  (ERSE),  Linden, 
New  Jersey.  Or.  M.  Lavid  was  Principal  Investigator  until  August  30,  1982, 
when  Mr.  I.  D.  Crane  assuned  responsibility  through  completion.  Dr.  A.  E. 
Cerkanowlcz  assisted  In  writing  this  final  report.  Experimental  studies  on 
catalytic  augmented  combustion  were  performed  by  Drs.  R.  J.  Gill  and  D.  B. 
Olson  under  a  subcontract  at  AeroChem  Research  Laboratories,  Inc.,  Princeton, 
New  Jersey. 


No  further  research  on  augmented  combustion  is  being  undertaken  at 
ERSE.  Therefore,  to  provide  a  perspective  on  the  total  effort,  this  report 
also  Includes  highlights  of  portions  of  the  work  done  previously. 
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ABSTRACT 

Two  novel  concepts  for  extending  aircraft  operational  range  have 
been  researched  under  this  contract.  They  are  radiative  and  catalytic 
augmentation  techniques.  The  objective  has  been  to  Investigate  the 
feasibility  of  utilizing  these  techniques  to  alter  (augment)  combustion 
Initiation  and  reaction  kinetics  which  restrict  combustor  operation  via  limits 
on  flammability,  flame  propagation.  Ignition  and  stability.  Both  techniques 
have  demonstrated  the  potential  to  Improve  combustion  processes.  The 
radiative  technique  under  laboratory  static  conditions  has  successfully 
Ignited  fuel-air  mixtures,  and  has  enhanced  combustion  processes,  utilizing 
pulsed  and  continuous  vacuum  ultraviolet  (VUV)  light  sources.  Similarly,  the 
catalytic  technique  has  broadened  stability  limits  under  normally  difficult 
fuel  lean,  low  temperature  conditions.  Complementary  efforts  have  focused  on 
the  development  of  analytical  capabilities  required  for  modeling  the  radiative 
and  catalytic  techniques. 

The  radiative  technique  Is  an  original  approach  which  can  Influence 
combustion  processes  by  the  employment  of  a  non-physical ly  Intrusive  beam  of 
light.  This  technique  utilizes  selected  wavelengths  In  the  VUV  region  to 
produce  oxygen  atoms  by  photodlssoclatlng  molecular  oxygen.  When  a  critical 
concentration  of  atomic  oxygen  Is  achieved  (about  10“  atoms/ m3) ,  combustion 
Initiation  occurs.  Subsequent  reactions  of  atomic  oxygen  with  fuel  molecules, 
as  well  as  with  other  combustion  species,  lead  to  Ignition  and  sustained 
combustion  via  chain  reactions.  Single  combustion  enhancement  without 
initiation  can  be  achieved  at  lower  oxygen  atom  concentration  levels. 

Radiative  Ignition  and  combustion  enhancement  tests  have  been  conducted  under 
various  static  conditions.  This  work  Included  VUV  excimer  laser  Ignition 
tests,  and  additional  Ignition  and  enhancement  experiments  employing  a 
continuous  light  source  (EIMAC).  The  VUV  output  of  the  continuous  light 
source  was  found  to  be  marginal  for  Ignition  as  well  as  for  enhancement. 
However,  It  was  successfully  used  to  demonstrate  the  feasibility  of  Igniting 
fuel-oxidant  mixtures  (propane/02)  using  continuous  or  steady-state  light 
sources. 

Catalytic  combustion  is  a  concept  wherein  combustion  reactions 
Initiated  by  a  heterogeneous  catalyst  play  an  iHportant  role  In  the  energy 
release  process  of  a  reacting  fuel-air  system.  This  part  of  the  program  has 
been  primarily  aimed  at  determining  the  feasibility  of  Improving  flame 
stabilization  and  reducing  pressure  loss  In  aircraft  afterburner  systems  by 
replacing  the  conventional  bluff-body  by  a  porous  catalytic  flameholder.  The 
bluff-body  stabilizes  the  combustion  by  forming  a  recirculation  zone  which.  In 
turn,  enhances  mixing  and  heat  transfer.  However,  It  also  causes  a 
substantial  pressure  loss.  By  using  a  porous,  catalytic  flameholder  the 
temperature  and  concentration  of  reactive  species  in  the  flameholder  wake  can 
be  "forced"  to  levels  favorable  for  stable  and  efficient  combustion. 
Experiments  In  a  contlnous  flow  combustion  test  facility  have  shown  that  a 
porous  catalytic  flameholder  can  provide  */*',e  stable  combustion  and  a  lower 
pressure  drop  than  a  bluff-booy  ...  -.oer.  This  indicates  that 
catalytlcally  coated  flameholders  may  be  a  promising  approach  to  enhanced 
stabilization  of  highly  diluted  partially  combusted  gas  streams  such  as  those 
exiting  the  main  burner  section  of  a  gas  turbine  engine. 
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SUMMARY 

Aircraft  operation  Is  limited  by  combustion  associated  phenomena 
Involving  Ignition,  flammability,  flame  propagation,  and  stable  combustion. 
Recent  combustion  research  has  Identified  a  nunber  of  potentially  promising 
techniques  for  rectifying  combustion-associated  aircraft  limitations.  The 
objective  of  this  work  Is  to  develop  an  understanding  of  basic  mechanisms  and 
processes  associated  with  two  new  techniques  which  effect  combustion 
augmentation  by  ultraviolet  Irradiation  or  catalytic  surfaces. 

Both  the  radiative  and  catalytic  augmented  combustion  concepts  have 
been  demonstrated  to  favorably  alter  combustion  processes.  The  radiative 
technique  offers  the  potential  for  dynamically  controlling  ignition  and  flame 
propagation  through  non-physlcally  Intrusive  means.  The  catalytic  technique 
Is  a  physically  Intrusive  but  passive  approach  which  can  broaden  flame 
stability  limits.  Additionally,  analytical  tools  have  been  developed  to  aid 
in  the  effective  transfer  of  these  concepts  to  practical  (nonlaboratory) 
combustion  systems. 

•  Radiative  Augmented  Combustion 

The  Ignition  and  combustion  enhancement  of  typical  gaseous  and 
liquid  fuel-air  mixtures  used  In  air-breathing  propulsion  engines  by  vacuum 
ultraviolet  (VUV)  irradiation  and  by  ultraviolet  (UV)  Irradiation  have  been 
the  principal  focus  of  the  radiative  augmented  combustion  concept.  Ignition 
and/or  enhancement  of  combustion  Is  achieved  by  selectively  producing  reactive 
combustion  Intermediary  species  via  a  photodlssociatl ve  mechanism  or  combined 
photo-thermal  methods  (1-3).  Implementation  of  the  technique  is  via  the  use 
of  a  non-lntrusive  (zero  fluid  dynamic  drag)  beam  of  light  In  a  selected 
wavelength  region. 

Figure  1  Illustrates  one  possible  application  configuration  where  a 
non-physical  flameholder  Is  generated  in  a  high  velocity  flow  by  focusing  the 
continuous  radiation  output  from  a  VUV  light  source  at  some  point  within  the 
fuel-air  flow  field.  A  "zero  pressure  drop”  flameholder  is  thus  established 
since  no  physical  obstruction  to  the  flow  exists.  Feasibility  of  this  concept 
has  been  established  In  the  laboratory  by  ignition  of  a  stationary  propane- 
oxygen  mixture  using  a  continuously  operated  light  source.  Previously  (1-3), 
only  pulsed  light  sources  were  sufficiently  powerful  in  VUV  radiation  to 
Ignite  fuel-air  mixtures. 

The  combustion  intermediate  species  of  major  importance  in  this 
approach  Is  atomic  oxygen.  The  photochemical  initiation  mechanism 
specifically  utilizes  the  efficient  photodissociation  of  oxygen  which  occurs 
for  radiant  energy  below  180  nm.  Subsequent  reaction  of  atomic  oxygen  with 
fuel  molecules  (as  well  as  other  combustion  species)  leads  to  sustained 
combustion  (Initiation)  when  oxygen  atom  concentrations  on  the  order  of  10 
atoms/m3,  or  greater,  are  generated.  Oxygen  atom  concentrations  below  the 
critical  level  will  not  result  In  combustion  initiation  but  can  lead  to 
combustion  enhancement.  By  photochemlcally  providing  an  excess  of  oxygen 
atoms  above  a  typical  level  of  about  10*2  atoms/m3,  the  thermal  reactivity  of 
a  fuel-air  mixture  can  be  altered  leading  to  a  change  in  flame  propagation 
rate.  Further,  since  reduced  amounts  of  oxygen  atoms  are  effective  for 
enhancement,  the  usable  radiant  wavelength  range  can  be  extended  to  about  200 
nm. 


Radiative  Flame-Stabilization  Geometry  (zero  pressure  drop  flame-holder) 


Previous  research  (1-3)  used  specially  designed  pulsed  plasma  short 
arcs  to  provide  necessary  VUV  radiant  energy  for  initiation  studies.  These 
earlier  studies  explored  the  influence  of  radiation  wavelength,  fuel  type, 
fuel -oxygen-nitrogen  mixture  ratios,  pressure,  tenperature,  and  mixture 
velocity  on  the  radiative  initiation  process.  More  recent  work  (10)  confirmed 
the  results  of  the  previous  work  with  respect  to  the  radiative  initiation 
mechanism  and  extended  light  source  considerations.  The  successful  ignition 
with  a  continuous  light  source  was  reported  for  the  first  time. 

In  the  most  recent  research,  new  ignition  experiments  with  pulsed 
VUV  laser  sources  and  additional  continuous  light  source  (EIMAC)  ignition 
tests  have  been  undertaken.  The  VUV  laser  ignition  work  was  conducted  by 
employing  an  excimer  laser  at  two  VUV  wavelengths;  157  rtm  (F~)  and  193  nm 
(ArF).  Radiative  ignitions  were  achieved  with  the  fluorine  faser  and  not  with 
the  argon  fluoride  laser,  even  though  the  latter  has  about  20  fold  greater 
fluence.  This  result  demonstrates  the  important  role  of  the  spectral 
selectivity  of  the  absorption  coefficient  of  molecular  oxygen.  Successful 
radiative  ignitions  with  a  continuous  light  source  (EIMAC)  were  achieved  on 
only  two  occasions  and  only  after  many  attenpts  due  to  the  marginal  VUV  output 
of  this  xenon  lamp. 

In  the  previous  efforts  a  comprehensive  computer  model  of  the 
photochemical  Initiation  mechanism  was  formulated.  This  model  includes 
extensive  kinetics  for  photon  absorption,  molecular  dissociation  and  reactant 
mixture  kinetics.  This  model  has  confirmed  and  guided  experimental  effort  and 
provides  a  base  for  the  designer  to  examine  new  ignition  methods  and 
applications.  In  the  most  recent  research,  this  model  has  demonstrated  good 
phenomenological  agreement  with  the  observed  ignition  behavior.  In  particular 
it  has  shown  good  quantitative  agreement  with  the  laser  ignition  results. 

This  reinforces  the  adequacy  of  the  model  for  elucidation  of  photon-radical 
Interaction  in  combustion  initiation. 

Experiments  exploring  the  radiative  enhancement  of  combustion  were 
conducted  with  propane-air  mixtures  and  a  continuous  UV  source.  Scoping 
studies  of  radiative  effects  on  flame  propagation  velocity,  extinction 
distance,  and  extinction  time  indicated  that  extension  of  the  normal 
flammability  limit  should  be  possible.  If  successful,  the  use  of  continuous 
VUV/UV  light  sources  would  offer  a  potential  opportunity  to  broaden  the 
current  flammability  envelope,  extending  aircraft  operational  range  beyond 
presently  available  flight  corridor. 

Although  its  VUV  output  is  marginal,  EIMAC  is  currently  the  best 
available  continuous  VUV  light  source  and  therefore  it  was  used  for  the  most 
recent  combustion  enhancement  work.  A  large  number  of  tests  investigating 
propane-air  flames  were  conducted.  Due  to  experimental  difficulties  and  the 
marginal  VUV  output,  the  data  obtained  were  limited  and  not  conclusive. 
Nevertheless,  the  observation  of  some  dramatic  changes  did  indicate  an  effect 
of  Irradiation  on  the  combustion  process.  Considerably  more  work  is  required 
to  elucidate  the  processes  involved. 

•  Catalytic  Augmented  Combustion 

In  comparison  with  the  radiative  technique,  catalytic  combustion 
augmentation  is  physically  intrusive  on  the  fuel-air  flow  field  but  does  not 


require  a  direct  energy  supply.  Combustion  is  promoted  by  use  of  a  solid 
catalytic  surface  which  accelerates  the  fuel-air  reaction  through  physical 
contact.  The  catalyst  is  deposited  on  the  surface  of  a  high-temperature 
ceramic  substrate  which  Is  in  the  form  of  a  honeycomb  placed  in  the  reactant 
flow  stream  in  the  manner  of  a  conventional  flameholder.  In  contrast  to  solid 
conventional  bluff-body  flameholders,  however,  a  portion  of  the  reactive 
mixture  will  pass  through  the  catalytic  unit.  Oxidation  of  the  fuel,  and 
energy  release  is  promoted  in  this  portion  of  the  mixture  due  to  surface 
activated  reaction.  The  increased  temperature  and  concentration  of  reactive 
combustion  intermediaries  provided  by  the  flow  exiting  the  catalytic  unit  is 
favorable  to  stable,  efficient  gas  phase  combustion  of  the  bulk  of  the 
reactant  mixture. 

Figure  2  illustrates  one  possible  catalytic  flame-stabilization 
geometry.  Only  a  portion  of  the  reactant  mixture  flows  through  the  catalytic 
flameholder  thereby  avoiding  the  concommitant  energy  loss  due  to  a  high 
pressure  drop.  This  approach  represents  a  significant  departure  from  other 
catalytic  combustion  techniques  where  the  catalytic  monolith  occupies  the  full 
flow  path. 

Flame  stabilization  tests  by  AeroChem  Research  Laboratories 
(Appendix  I)  in  a  continuous  flow  combustion  facility  have  demonstrated  that 
porous  catalytic  flameholders  may  be  a  promising  approach  to  augment 
combustion  stability  limits.  During  the  contract  period,  three  flameholders 
of  the  same  cross  sectional  area  -  a  solid  bluff  body,  a  porous  catalytic 
structure,  and  a  porous  noncatalytlc  structure  -  were  characterized 
experimentally.  Results  at  low  reactant  velocity  (below  13  m/s)  show 
conclusively  that  the  fuel  lean  blowoff  limits  are  wider  with  a  catalytic 
flameholder  than  with  noncatalytic  but  otherwise  Identical  flameholders.  More 
Importantly,  for  velocities  above  11  m/s,  the  porous  catalytic  flameholder 
exhibited  a  wider  lean  blowoff  limit  than  a  solid  flameholder.  Additionally, 
the  porous  catalytic  flameholder  exhibited  a  lower  pressure  drop.  High  output 
combustors  or  afterburners  require  flame  stabilization  in  a  high  velocity 
fuel -air  mixture  and  their  performance  would  benefit  from  the  increased 
blowoff  limit.  Extension  of  the  experimental  data  for  velocities  up  to  200 
m/s  is  needed,  however,  to  verify  and  quantify  the  departure  in  blowoff  limit. 

Two  models  have  been  developed  for  specific  aspects  of  the  catalytic 
flame  stabilization  concept.  The  first  provides  for  description  of  the 
aerothermochemlstry  within  the  catalytic  unit  and  has  been  demonstrated  to 
correctly  predict  the  catalytic  performance.  The  second  deals  with  the  flow 
split  through  and  bypassing  the  catalytic  unit  and  the  wake  mixing  region 
downstream  of  the  catalytic  unit.  A  first  order  description  of  catalytic 
flameholding  Is  possible  using  a  combination  of  these  models.  Comparison  with 
experimental  data  Is  required  to  establish  the  model  parameters. 


Figure  2.  Catalytic  Flame-Stabilization  Geometry. 


DISCUSSION 


RADIATION  AUGMENTED  COMBUSTION 


Various  radiative  sources  were  experimentally  characterized  or 
evaluated  with  regard  to  photochemical  combustion  initiation  and  enhancement 
over  a  wide  range  of  operating  conditions.  Three  different  types  of  vacuum 
ultraviolet  (VUV)  igniters  were  effectively  utilized:  (1)  pulsed  plasma  short 
arcs,  (2)  focused  continuous  short  arcs,  and  (3)  pulsed  lasers.  Table  I  lists 
the  optical  features  of  these  light  sources.  A  continuous  ultraviolet  (UV) 
short  arc  source  was  also  employed  in  combustion  enhancement  experiments. 
Photochemical  mitigations  were  evaluated  by  measuring  the  stored  energy 
required  for  successful  combustion  initiation.  Measurements  of  flame  velocity 
and  reactant  burn-out  time  were  used  to  characterize  photochemical 
enhancement. 

•  Photochemical  Ignition  Experiments 

Detailed  spectroscopic  data  was  obtained  on  the  pulsed,  plasma  short 
arc  in  the  vacuum  ultraviolet  and  ultraviolet  wavelength  regions  (5).  This 
particular  source  has  been  used  for  the  majority  of  the  photochemical  ignition 
parameterization  experiments.  Figure  3  illustrates  the  spectral  distribution 
of  intensity  for  this  source  which  exhibits  pronounced  peaks  in  the  145  to  175 
nm  wavelength  region.  This  Is  a  region  of  photon  energy  which  is  extremely 
efficient  in  initiating  combustion  reactions.  Spectral  data  also  indicated 
that  the  temporal  behavior  of  the  light  pulse  can  be  effectively  modeled  as  a 
critically  damped  pulsed  of  the  form 

I/Imax  =  (t/t*)exp(l-t/t*) 

where  I  =  pulse  Intensity,  Imax  =  maximum  intensity,  t  »  time,  and  t*  *  time 
to  peak  intensity.  Propane-air  ignition  experiments  with  this  radiation 
source  reconfirmed  the  previously  observed  photochemical  ignition 
characteristics. 

The  second  type  of  light  source  successfully  employed  in  radiative 
augmented  combustion  experiments  was  a  commercially  available  EIMAC  lamp. 

This  was  a  continuous  xenon  point  source  with  a  low  grade  sapphire  window  and 
an  elliptical  reflector.  Blackbody  radiant  energy  output  in  the  150-400  nm 
spectral  region  was  developed  at  input  power  levels  of  up  to  600  watts. 
Ignition  of  propane-oxygen  mixtures  at  ambient  temperature  and  atmospheric 
pressure  was  achieved  at  equivalence  ratios  of  0.8  and  1.0  with  power  inputs 
of  615  watts  and  500  watts.  Although  the  light  source  was  destroyed  in  each 
case,  this  represents  the  first  demonstrations  of  radiative  ignition  using 
continuous  light  source  irradiation.  Continued  development  of  this  capability 
would  provide  the  option  of  designing  an  optical-radiative  flameholder  with 
zero  pressure  drop  due  to  physical  obstruction. 

Laser  ignition  experiments  were  conducted  by  the  principal 
Investigator  (Dr.  Lavld)  at  Los  Alamos  Scientific  Lab/Applied  Photochemistry 
Division  using  their  Lumonlcs  TE861  laser  facility.  The  lasers  were  used  at 
two  wavelengths;  193  nm  (ArF)  and  157  nm  (Fg)  with  measured  fluence  at  about 
300  to  400  mj/cnr  and  20  to  40  mJ/cnr  respectively.  Hydrogen-oxygen-nitrogen 
mixtures  at  varying  equivalence  ratio  and  subatmospheri c  pressure  were  used. 
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Figure  3  -  Spectral  Distribution  of  Intensity 
for  Pulsed  VUV  Light  Source 


Photochemical  ignition  was  not  achieved  using  the  ArF  laser.  This  is 
consistent  with  previous  results  (1)  which  indicated  that  wavelengths  below 
160  to  170  nm  were  the  most  effective  in  producing  radiative  ignition.  Use  of 
the  Fo  ^as«r  proved  successful  although  the  fluence  was  about  1/20  that  of  the 
ArF  laser. 

Two  characteristics  of  the  F2  laser  need  to  be  discussed  before  the 
results  of  Ignition  experiments  can  be  understood.  First,  laser  pulse  energy 
decreases  with  use  as  shown  in  Figure  4,  necessitating  that  the  F--He  mixture 
in  the  laser  chamber  be  refreshed  on  a  regular  basis.  More  importantly,  the 
focal  length  of  the  laser-optical  system  is  wavelength  dependent.  This  is 
clearly  Illustrated  by  the  results  in  Figure  5,  which  show  relative  pulse 
energy  vs.  focal  distance.  Exclusion  of  VUV  Irradiation  is  achieved  by 
purging  the  optical  path  with  oxygen.  The  focal  point  when  the  total  energy 
flux  Is  present  is  about  24.6  cm  while  the  optimum  focal  point  for  the  VUV 
occurs  at  about  19.7  cm.  When  operated  properly,  laser  ignition  of  hydrogen- 
oxidant  mixtures  were  achieved  on  a  regular  basis.  These  results  confirm  the 
radiative  ignition  concept  and  extend  the  repertoire  of  available  light 
sources  usable  for  photochemical  combustion  initiation. 

•  Combustion  Enhancement  Experiments 

The  ultraviolet  capabilities  of  the  EIMAC  source  were  also  utilized 
In  exploring  combustion  enhancement.  In  these  experiments,  propane/air 
mixtures  at  various  equivalence  ratios  and  pressure  were  subjected  to  UV 
Irradiation  simultaneous  with  flame  propagation  (9,11).  The  effect  of 
propagation  toward  and  away  from  the  light  source  was  studied.  Since 
comparison  with  literature  quoted  flame  speeds  Is  highly  dependent  on  specific 
conditions  such  as  direction  of  gravity,  heat  loss,  confinement,  etc., 
complementary  tests  In  the  absence  of  light  were  also  performed.  The 
clnemagraphlc  technique  employed  precluded  accurate  measurement  at  near 
stoichiometric  ondltlons  due  to  Intense  light  from  the  combustion  reactions. 
Tables  II  and  III  highlight  some  of  the  results  of  this  work.  A  more 
comprehensive  treatment  of  the  results  Is  given  in  Refs.  9  and  11. 

Table  II  Indicates  that  a  substantially  positive  effect  on  both 
extinction  time  and  distance  occurs  due  to  light  irradiation  in  front  of  a 
propagating  flame.  This  effect  was  most  pronounced  for  fuel  lean  mixtures 
where  up  to  an  81%  Increase  In  extinction  time  and  25%  increase  in  extinction 
distance  were  measured.  These  measurements  were  only  taken  with  the  light 
passing  through  the  unburned  mixture  toward  the  propagating  flame. 

Additional,  more  extensive  enhancement  experiments  have  been  done 
Investigating  propane-air  flames  using  the  EIMAC  light  source.  Although  its 
output  Is  marginal,  the  EIMAC  Is  currently  the  best  available  continuous  light 
source.  Emphasis  was  placed  on  measuring  flame  propagation  velocity  to 
evaluate  combustion  enhancement.  While  no  extinction  distances  or  times  were 
recorded  In  these  more  recent  experiments,  the  scope  of  the  work  was  broadened 
to  Include  two  different  optical  paths  between  the  flame  front  and  the  light 
source.  Flame  propagation  velocities  were  primarily  affected  when  the  source 
was  behind  the  flame  front  and  Irradiation  was  passing  through  the  combustion 
products,  ftje  to  experimental  difficulties  and  the  marginal  VUV  output,  the 
data  obtained  were  not  conclusive.  As  Indicated  In  Table  III,  some  conditions 
lead  to  inhibition  rather  than  enhancement.  However  the  observation  of  some 
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dramatic  changes  does  Indicate  an  effect  of  Irradiation  on  the  combustion 
process.  Considerably  more  work  needs  to  be  done  to  define  the  processes 
Involved. 

•  Analytical  Modelling 

Analytical  work  was  directed  toward  elucidating  the  fundamental 
Interaction  of  photon  absorption,  subsequent  dissociation,  the  chemical 
kinetics  of  the  reactant  mixtures,  and  possible  heat  loss  to  the 
surroundings.  Preliminary  modelling  results  discussed  in  Ref.  4  and  Appendix 
II  are  summarized  with  additional  details  here.  Currently,  the  chemical 
kinetics  Included  allow  treatment  of  hydrogen-oxygen -mixtures  with  the 
possible  presence  of  an  Inert  diluent.  The  reactant  species  which  are 
considered  are  presented  in  Table  IV.  The  presence  of  electronically  excited 
state  species,  not  typical  of  most  hydrogen-oxygen  schemes,  is  dictated  by 
their  photochemical  production,  and  although  their  Inclusion  Is  not 
Indlspenslble,  0(1D)  will  be  seen  to  play  a  detectable  role.  Some  90  chemical 
reactions  which  occur  between  these  species  are  Included  in  the  model.  Rate 
constants  in  the  form  AT"8  exp  (-C/T)  were  obtained  from  the  extensive 
literature  on  the  subject,  with  emphasis  on  data  In  the  Initiation  and  early 
Ignition  temperature  range  (300-1200  K).  The  rate  data  tabulated  In  Ref.  1 
were  used,  with  some  values  updated  to  reflect  those  used  In  recent  reviews 
(7,15). 

Radiant  absorption  Is  assumed  to  follow  a  Beer-Lambert  law.  The 
absorbing  species  Include  O*.  0-,  H-0,  HO-  and  H-O*.  Upon  absorbing  a  photon, 
these  species  are  assumed  to  dissociate  to  products  which  depend  on  the  energy 
of  the  photon.  Thus,  quantum  yields  as  a  function  of  wavelength  have  been 
taken  to  be  one  over  a  given  Interval  and  zero  otherwise.  The 
photodissociation  reactions  considered  are  given  In  Table  IV,  together  with 
the  relevant  wavelength  Intervals.  Photodissociation  can  thus  be  handled  as  a 
uni -molecular  reaction  In  which  the  dissociation  rate  Is  obtained  as  the 
product  of  the  concentration  of  the  absorbing  species  and  a  pseudo- first-order 
rate  constant,  which  depends  on  the  Incident  spectral  Irradlance,  I(x,t,Y),  at 
a  given  point  In  the  mixture  and  at  a  given  time.  Because  the  objective  Is  to 
gauge  the  robustness  of  the  Ignition  kernel,  which  lies  close  to  the  combustor 
window,  the  spatial  dependence  Is  not  considered  and  the  spectral  Irradlance 
used  Is  that  which  would  be  measured  at  the  mixture  side  of  the  window.  The 
number  of  Incident  photons  of  wavelength  X  to  X+dy  per  unit  area  time  Is  thus 
I(x/hc)l(t,x).  When  multiplied  by  the  wavelength-dependent  cross-section  of 
the  absorber  and  Integrated  over  the  relevant  wavelengths,  we  thus  obtain  the 
desired  rate  constant  for  photodissociation  Is  obtained.  The  integral  is.  In 
fact,  approximated  as  a  sum  with  one  nanometer  wavelength  Increments,  unlike 
the  chemical  rate  constants.  It  is  time -dependent  through  the  temporal 
dependence  of  I(t,X).  This  dependence  was  taken  to  be  of  the  form  I(t,x)  * 
Imax(X)P(t);  where  P  (t)  Is  the  pulse  shape  function,  generally  taken  to  be 
that  of  a  critically  damped  pulse.  Tabulated  absorption  coefficients  for  all 
absorbing  species  are  used;  absorption  of  0?  In  the  Schumann-Runge  bands 
follows  the  treatment  of  G.  Kocharts  (7).  For  each  species,  a  rate  balance 
was  formed  by  consideration  of  all  chemical  and  photochemical  reactions  in 
which  It  participates  (4). 


A  heat  balance  equation  accounts  for  heat  liberated  or  removed  by 
chemical  reaction  and  possible  heat  loss  to  surroundings.  The  heat  loss  has 
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been  taken  proportional  to  the  difference  between  gas  and  ambient 
temperatures,  with  heat  loss  coefficient  calculated  from  he  *  16  k/SCDd2;  k  * 
thermal  conductivity.  The  diameter  (d)  is  estimated  for  the  ignition^kernel 
by  using  the  Irradiated  cross-sectional  area.  The  model  is  finally 
constituted  as  an  initial  value  problem  for  the  resulting  system  of  nonlinear 
ordinary  differential  equations.  Solution  techniques  are  documented  in  Ref. 

4.  Diffusion  and  convective  effects  have  been  neglected. 

The  resulting  system  of  differential  equations  is  autonomous  and,  as 
recognized  by  Yang  and  Gray  (6),  it  is  profitable  to  consider  the  trajectories 
determined  by  the  solutions  In  phase  space  and  the  singular  points  of  the 
system.  Reference  to  multi-dimensional  phase  space  is  cumbersome,  and  for 
combustion  Initiation  In  the  iWO*  system,  largely  unnecessary.  It  suffices 
to  present  the  results  in  the  temperature  vs.  oxygen  atom  concentration  plane, 
with  the  oxygen  radical  acting  as  a  surrogate  for  the  other  radicals.  When  so 
done,  the  trajectories  reveal  a  phase  plane  typical  of  a  simple  critical  chain 
branching  system,  as  shown  In  Figure  6.  The  unlabeled  solid  trajectories 
represent  solutions  for  varying  initial  temperature  and  oxygen  atom 
concentration.  The  dashed  line  Indicates  the  separatrix  between  the  fast- 
reaction  region,  C,  and  the  noncombusting  region,  N.  This  separatrix  is,  of 
course,  composed  of  the  two  trajectories  which  lead  to  the  saddle  point,  S, 
heurlstically  the  critical  point  arising  from  the  balance  of  chain  branching 
and  chain  termination.  Trajectories  which  lie  below  this  separatrix  lead  to 
the  stable  node.  A,  at  ambient  temperature.  Trajectories  which  lie  above  the 
separatrix  are  characterized  by  rapid  chain-branching  and  Increase  in 
temperature  and  thereby  Indicate  Ignition. 

Thus,  In  the  absence  of  Irradiation,  the  reactant  mixture  which 
begins  at  a  point  on  a  given  trajectory  will  tend  by  chemical  action  to  follow 
that  trajectory.  Gonslder  now  the  effect  of  Irradiation  and  the  consequent 
production  of  oxygen  atoms.  The  resulting  solution  curve  is  no  longer  one  of 
the  reaction -determined  trajectories  but  departs  from  these  trajectories  in 
the  direction  of  Increasing  oxygen  atom  concentration.  Such  a  path  is  shown 
as  curve  1  In  Fig.  6.  It  Is  the  resultant  of  oxygen  atom  increase  due  to 
photodissociation  and  temperature  increase  due  to  chemical  reaction.  In  this 
case,  the  separatrix  Is  crossed  and  Ignition  results.  Curve  2,  however, 
represents  a  solution  In  which  a  radiant  pulse  has  effectively  ceased  before 
the  separatrix  has  been  crossed;  the  curve  thus  rejoins  a  trajectory  leading 
back  to  the  stable  node  and  no  Ignition  results.  In  contrast,  other  ignition 
approaches  (e.g. ,  "pure"  thermal,  spark)  rely  slgni flcantly  on  departure  from 
the  trajectories  In  the  direction  of  Increasing  temperature  and  follow  paths 
which  can  be  expected  to  differ  from  that  of  radiative  initiation.  7his 
analysis  provides  a  qualitative  framework  for  understanding  photochemical 
Initiation  and  Indicates  the  value  of  the  phase  plane  for  presenting  the 
quantitative  modeling  results  to  be  discussed. 

•  Comparison  of  Model  and  Experiment 

Application  of  the  model  to  pulsed  plasma  short  arc  initiated 
reactions  has  been  extensively  discussed  before  (Refs.  1,  4,  and  Appendix 
II).  The  model  has  proved  particularly  useful  in  understanding  the  laser 
Ignition  experiments.  Analysis  of  laser  Ignition  characteristics  indicates 
that  with  the  currently  available  fluences  the  shorter  wavelength  laser  (F?) 

Is  capable  of  achieving  photochemical  Ignitions,  while  the  ArF  laser  is  not. 
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LOG  OF  OXYGEN  ATOM  CONCENTRATION 


FIGURE  6.  REPRESENTATIVE  PHASE  PLANE  FOR  PHOTOCHEMICAL 
INITIATION  OF  COMBUSION. 

S:  SADDLE  POINT 

A:  NODAL  POINT  AT  AMBIENT  TEMPERATURE 
— :  TRAJECTORIES  IN  THE  ABSENCE  OF  IRRADIATION 
— :  SEPARATRIX  OF  COMBUSTING  AND  NON-COMBUSTING  REGIONS 

- 1:  SOLUTION  CURVE  FOR  IRRADIATION  LEADING  TO  IGNITION 

- 2;  SOLUTION  CURVE  FOR  IRRADIATION  WITHOUT  IGNITION 


This  finding  Is  in  agreement  with  our  numerical  predictions,  as  depicted  in 
Fig.  7.  The  bars  shown  Indicate  the  minimum  fluence  predicted  to 
photochemlcally  Ignite  a  premixed  hydrogen-oxygen  mixture  at  an  equivalence 
ratio  of  0.6  and  at  20  kPa  for  the  wavelengths  of  the  (157  nm)  and  ArF  (193 
nm)  lasers.  These  values  were  calculated  from  the  radiative  ignition  model  by 
Increasing  the  value  for  the  fluence  (i.e.  irradiance,  as  the  pulse  width 
remains  the  same)  until  a  solution  indicating  ignition  was  obtained.  The 
resulting  fluences  are  36  and  2.7  x  105  mj/crtr,  respectively.  This 
significant  difference  in  fluences  (and  therefore  energies)  is  primarily  due 
to  the  spectral  selectivity  of  the  absorption  coefficient  of  the  oxygen 
molecule.  At  157  nm,  which  Is  in  the  Schumann-Runge  continuum,  the  absorption 
coefficient  is  167  cm-1,  and  at  193  nm,  which  is  in  the  S-R  bands,  it  is  0.027 
cm"1,  about  four  orders  of  magnitude  lower  (7).  Experimental  results 
indicated  by  the  symbols  for  Ignition  and  no-ignition  events  agree  with  these 
predictions.  Although  the  fluence  of  the  F2  laser  is  considerably  smaller 
than  that  of  the  ArF  laser.  It  successfully  achieved  ignition,  while  the  ArF 
failed  in  photochemlcally  igniting  the  same  gaseous  mixtures. 

Minimum  ignition  energies  for  combustible  gaseous  mixtures  are 
typically  given  as  a  function  of  equivalence  ratios  at  various  pressures. 

Since  it  was  very  difficult  experimentally  to  vary  accurately  the  energy 
emitted  by  the  laser  and  absorbed  by  the  mixture,  this  was  kept  constant  and 
the  pressure  of  the  mixture  was  varied  Instead.  Figure  8  depicts  the  minimum 
pressure  at  which  the  mixture  is  photochemical ly  Igni table  by  the  F-  laser  vs. 
equivalence  ratio  for  hydrogen-oxygen  and  hydrogen-air  mixtures.  Tne  curves 
In  the  figure  are  derived  from  analytical  predictions  for  minimum  ignition 
pressures  which  were  obtained  by  fixing  other  model  inputs  to  values 
corresponding  to  those  of  the  experimental  conditions.  In  particular,  the 
value  used  for  the  heat  loss  coefficient  was  based  on  the  Irradiated  cross- 
section;  It  was  then  varied  only  according  to  Its  pressure  dependence.  A 
search  was  performed  until  two  pressures  were  obtained  which  differed  by  1 
kPa,  with  Ignition  Indicated  for  the  higher  and  not  the  lower.  For  and 

equivalence  ratio  of  0.6,  such  solution  curves  are  displayed  on  the  phase 
plane  in  Fig.  9.  It  Is  interesting  to  compare  this  figure  to  Fig.  6  and  to 
observe  the  expected  cusp-like  behavior  in  the  neighborhood  of  the  saddle 
point.  This  procedure  was  carried  out  for  a  sufficient  number  of  equivalence 
ratios  to  allow  the  curves  to  be  drawn  accurately. 

The  experimental  results  in  Fig.  8  are  indicated  with  bars  to 
represent  the  uncertainty  In  the  results  obtained  due  to  fluctuations  in  laser 
output.  The  agreement  between  analytical  predictions  represented  by  the  solid 
lines  and  experimental  results  shown  as  data  points  is  very  satisfactory.  As 
expected,  the  results  show  that  it  Is  easier  to  ignite  Hj/C^  than  H2  air 
mixtures.  The  primary  reason  for  the  increase  in  pressure  needed  for  the  air 
system  to  Ignite  is  the  necessity  to  increase  molecular  oxygen  concentration 
so  as  to  achieve  the  required  production  of  oxygen  atoms  upon  irradiation. 
Indeed,  the  critical  O2  concentration  somewhat  exceeds  that  for  the 
corresponding  H2/0?  mixture  due  to  the  presence  of  the  diluent,  N2,  resulting 
In  an  Increase  fn  thermal  inertia,  and  to  the  deleterious  effects  of 
Increasing  pressure,  e.g.  Increased  rate  for  certain  chain  termination 
paths.  Also,  It  Is  Important  to  note  that  experimental  and  numerical  results 
Indicate  that  fuel-lean  mixtures  with  some  excess  oxygen  are  more  favorable 
for  photochemical  Ignitions  than  other  mixtures,  with  the  most  favorable 
mixture  at  an  equivalence  ratio  of  0.6.  This  finding  is  quite  different  from 
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OXYGEN  ATOM  CONCENTRATION  (particlas/cm3) 


Figure  9.  Solution  curves  for  the  determination  of  minimum  ignition 
pressure  with  fluorine  laser. 
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that  for  spark  ignition  in  which  minimum  ignition  is  found  to  occur  at  about 
equivalence  ratio  of  unity  (12).  Again,  the  explanation  is  that  increased 
molecular  oxygen  results  in  higher  atomic  oxygen  concentration.  However,  when 
equivalence  ratio  is  further  reduced,  the  rate  of  chemical  heat  release  is 
decreased  due  to  the  lessened  availability  of  fuel  molecules,  which 
necessitates  higher  minimum  pressure  for  ignition. 

Oetailed  analysis  of  the  combustion  enhancement  experiments  would 
have  required  a  major  modificatin  of  the  code  to  allow  for  flame 
propagation.  Therefore  a  quite  different  nunerical  model  (14,16),  developed 
at  Lawrence  Livermore  Laboratory,  was  used  to  examine  the  enhancement  results 
obtained  in  a  propane-air  flame.  This  model  allows  prediction  of  laminar 
flame  properties.  It  consists  of  a  one-dimensional,  time -dependent  treatment 
of  the  conservation  equations,  together  with  detailed  chemical  kinetic 
mechanism  including  photodissociation  reactions.  This  attempt  proved 
unsuccessful  (11). 
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CATALYTIC  AUGMENTED  COMBUSTION 


The  objective  of  this  research  was  to  evaluate  the  feasibility  of 
using  catalytic  honeycomb  flameholders  to  expand  flame  stabilization  limits 
for  potential  applications  such  as  in  aircraft  afterburners.  This  was 
performed  by  experimentally  evaluating  flame  stablization  using  catalytic  and 
noncatalytic  flameholders  in  a  research  scale  continuous  flow  tube 
combustor.  Experiments  were  conducted  using  premixed  propane/air  to  determine 
the  effects  of  fuel/air  ratio,  flow  velocity,  and  flameholder  type  on  the 
flame  stability  limits.  Several  other  diagnostics  were  also  used  to  gauge  the 
efficiency  of  the  flame  stabilization  and  combustion  process.  These 
diagnostics  Included  gas  tenperature  profiles  in  the  combustor,  flameholder 
temperatures,  pressure  drop  across  the  flameholders,  and  CO  and  C02 
concentrations  in  the  burned  combustion  gases. 

•  Flame  Stabilization  Experiments 

The  experiments  on  catalytic  flame  stabilization  were  carried  out  by 
AeroChem  Research  Laboratories.  A  report  covering  their  research  is  provided 
in  Appendix  I.  The  combustor  used  for  these  experiments  is  shown  In  Figure 
11.  The  test  section,  which  contained  the  flameholder,  had  5.3  cm  internal 
diameter. 


Three  types  of  cylindrical  flameholders  were  used:  solid  brass, 
noncatalytic  ceramic  honeycomb  and  catalytic  ceramic  honeycomb.  The  ceramic 
flameholders  were  fabricated  from  Cordlerlte  monolithic  cylindrical  shaped 
honeycomb  plugs.  The  catalytic  flameholders  were  custom  prepared  by 
impregnating  the  Cordlerlte  with  palladium  catalyst.  Several  different 
diameter  and  cell  densities  were  screened.  The  catalytic  flameholder  for  the 
majority  of  the  experimental  runs  had  a  2.0  cm  outside  diameter  ard  a 
honeycomb  density  of  31  cells/cm  . 

It  was  expected  that  a  larger  pressure  difference  would  exist  across 
a  solid  body  flameholder  than  across  a  more  porous  body,  such  as  ceramic 
honeycomb.  Also  higher  gas  flow  velocities  were  expected  to  increase  the 
pressure  difference  across  the  bodies.  Both  of  these  effects  were  confirmed 
in  room  temperature  flow.  Experimental  data  for  a  solid  bluff  body,  a  fully 
open  honeycomb,  and  a  partially  blocked  honeycomb  are  shown  in  Figure  11. 

While  the  magnitude  of  the  pressure  drops  were  small  due  to  the  scale  being 
tested,  a  substantial  difference  existed  between  the  flameholders  (e.g.  the 
pressure  drop  for  a  solid  body  was  approximately  65X  greater  than  for  an  open 
honeycomb  at  12  m/sec  velocity). 

Fuel  lean  blowoff  limits  were  measured  to  assess  the  effect  of  the 
porous  catalytic  flameholders  upon  the  flame  stability  limits.  Reactant  flow 
velocities  were  varied  from  9  to  13  m/s  with  equivalence  ratios  varying  from 
0.6  to  0.9  on  the  fuel  lean  side  of  stoichiometric.  The  experimental  results 
illustrated  in  Figure  12  compare  the  lean  blowoff  limits  for  the  three  types 
of  flameholders  tested.  Blowoff  limits  are  significantly  wider  for  the 
catalytic  honeycomb  than  for  the  noncatalytic.  In  addition  the  limiting 
equivalence  ratio  for  the  catalytic  flameholder  decreases  (i.e.  more  lean 
conditions)  at  higher  flow  velocities.  This  was  an  unexpected  result  since 
normally  solid  body  limits  shift  toward  stoichiometric  with  Increasing 
velocity.  For  these  experiments  the  catalytic  body  showed  a  wider  limit  than 
the  solid  body  above  a  velocity  of  about  11  m/sec. 


FIGURE  10.  FLOW  TUBE  COMBUSTOR 
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O  ■  solid  bluff  body;  O  ■  partially  blocked  ceramic 
honeycomb  (31  cells  cm”2)  with  89  cells  open;  A  ■ 
ceramic  honeycomb  with  157  cells  open.  All  results  for 

300  K  air  flows. 


FIGURE  11.  PRESSURE  DROPS  FOR  2.79  cm  o.d.  FLAMEHOLDERS 
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FIGURE  12 
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.  COMPARISON  OF  BLOWOFF  LIMITS  FOR  2.79  cm  diam 
CERAMIC  HONEYCOMB  FLAMEHOLDERS  AND  2.87  cm  diam 
SOLID  BODY  FLAMEHOLOER 


a  *  catalytic 
o  *  noncatalytic 
Dashed  line  *  solid  body 


The  findings  of  this  program  show  that  catalytically  coated  ceramic 
(or  other  material)  flameholders  may  be  a  promising  approach  to  enhanced 
stabilization  of  highly  diluted  partially  combusted  gas  streams  such  as  those 
exiting  the  main  burner  section  of  a  gas  turbine  engine.  The  conditions 
necessary  for  efficient  combustion  while  employing  catalytic  flameholders 
must,  however,  be  more  thoroughly  investigated  before  practical  catalytic 
combustors  can  be  manufactured. 

Primary  practical  disadvantages  of  the  catalytic  flameholders  used 
In  this  feasibility  study  appear  to  be  associated  with  the  ceramic  substrates, 
which  In  the  case  of  Cordlerlte,  restrict  operation  to  rather  low 
temperatures.  In  the  case  of  some  other  ceramics  suitable  for  higher 
temperatures,  the  substrates  are  susceptible  to  fracture  due  to  thermal  shock 
and  are  less  porous.  The  palladium  catalyst  Investigated  in  this  program  is 
also  probably  not  practical  In  view  of  its  significant  vapor  pressure  at 
temperatures  above  1500  K. 

In  view  of  the  progress  achieved  on  this  program  and  the  potential 
advantages  of  catalytic  enhanced  flame  stabll Izations,  the  following 
recommendations  are  made: 

1.  A  similar  program  should  be  performed  to  compare  catalytic  and 
noncatalytlc  honeycomb  flameholders  In  a  larger  diameter 
combustor  and  at  higher  velocities.  The  larger  combustor  would 
reduce  wall  Interactions  and  a  wider  range  of  velocities  would 
better  define  the  variation  of  blowoff  limits  with  velocity. 

2.  Further  studies  of  the  CO,  CO?,  and  hydrocarbon  emissions  should 
be  performed  along  with  (1)  above. 

3.  Catalytic  flameholder  substrates  other  than  the  Cordierite  used 
In  this  work  should  be  tested.  These  substrates  would  include 
other  coated  ceramics,  coated  metals,  and  ceramics  incorporating 
the  catalysts  Into  their  lattice. 

4.  Catalytically  active  coatings  other  than  palladium  should  be 
tested. 

5.  Tests  should  be  performed  under  conditions  of  pressure, 
velocity.  Inlet  temperature,  and  fuel  air  ratio  which  more 
closely  simulate  the  practical  combustors  of  interest,  and  jet 
fuels  should  be  used  for  the  tests. 


The  above  recommendations  define  the  direction  of  future  work  but 
are  not  meant  to  be  conprehensl ve.  Whether  this  technique  can  eventually  be 
practical  Is  not  clear  at  this  time  but  an  evaluation  can  probably  be  made 
after  some  of  the  recommended  additional  experiments  have  been  performed. 

•  Analytical  Modelling 

Mathematical  modeling  of  catalytic,  porous  flameholders  was  divided 
Into  two  parts.  In  the  first  part,  conditions  within  the  porous  flameholder 
were  modeled.  Here  catalytic  reactions  are  dominant.  The  second  part 
considers  the  slow  split  around  and  through  the  flameholder  as  well  as  the 
critical  wake/recirculation  region  where  the  actual  flame  stabilization 
occurs. 


Description  of  the  internal  catalytic  combustion  processes  is  the 
most  advanced  and  tested  part.  Flexibility  and  utility  is  provided  by  the 
simultaneous  consideration  of  the  relevant  chemical  and  physical  processes 
using  a  one-dimensional,  steady  state,  plug  flow  model  with  global  reaction 
rates.  Heterogeneous  (catalytic)  and  homogeneous  (gas-phase)  reactions, 
velocity  variation,  and  axial  heat  conduction  in  the  catalytic  substrate  are 
Included  as  are  diffusion  and  temperature-dependent  surface  kinetics  for  the 
heterogeneous  reactions.  Combustion  initiation  (light-off)  and  extinction 
limits  are  particularly  dependent  on  the  processes  included  in  the  model. 
Appendices  II  and  III  provide  a  more  detailed  description  of  this  part  of  the 
model.  In  particular.  Appendix  III  demonstrates  the  excellent  agreement 
between  the  catalytic  model  and  experimental  results. 

The  second  part  of  the  modeling  effort  considers  the 
aerothermochemlstry  occurring  external  to  the  catalytic  flameholder  as 
depicted  conceptually  In  Fig.  13.  This  model  is  based  on  the  usual  wake 
recirculation  formulation,  as  influenced  by  the  work  of  Zukoski  and  Marble 
(17,18)  on  critical  ignition  times.  The  split  of  flow  through  and  around  the 
catalytic  monolith  Is  determined  by  pressure  drop  considerations.  The  flow 
Into  the  recirculation  zone  is  thus  the  sum  of  the  flow  through  the  monolith 
and  that  resulting  from  turbulent  inclusion.  Even  in  the  absence  of  catalytic 
effects,  the  need  to  account  for  flow  through  the  flameholder  represents  a 
departure  from  the  usual  solid  body  flameholder  description.  The  distance 
needed  for  a  flame  to  develop  Is  predicted  as  a  product  of  mean  convective 
velocity  and  reaction  time  obtained  from  the  analysis  of  Kunder  et  al.  (8)  and 
compared  to  a  recirculation  zone  length  calculated  on  the  basis  of  flow 
parameters.  Comparison  of  these  lengths  allows  the  determination  of  blowout 
or  flameholding. 

In  this  way,  this  model  allows: 

1.  Description  of  the  influence  of  overall  geometry,  blockage 
ratio,  and  substrate  configuration  on  aerodynamic  and 
thermodynamic  variables,  especially  pressure  drop. 

2.  Determination  of  the  degree  of  approach  flow  penetration  through 
the  catalytic  section,  recirculation  zone  size  and  shape,  and 
shear  layer  characteristics. 

3.  Description  of  the  concept's  ability  to  widen  stability  and 
ignition  limits  and  to  widen  the  cross-sectional  area  occupied 
by  the  stabilization  devices. 

Additional  details  on  this  flameholder  model  development  is  provided  in 
Appendix  IV. 

If  further  research  on  the  porous  catalytic  flameholder' s  concept  is 
undertaken,  a  full  model  combining  the  two  models  (part  1  and  2)  should  be 
developed.  Together  with  experimental  results,  this  complete  model  would  form 
the  basis  for  performance  optimization  and  extension  of  the  work  to  practical 
combustor  systems. 


CONCLUSIONS  AND  FUTURE  WORK 


Both  the  radiative  and  catalytic  techniques  have  demonstrated  the 
capability  to  enhance  combustion  processes  and  to  broaden  normally  encountered 
stability  1 imits. 

The  work  on  radiative  ignition  and  combustion  enhancement  is 
providing  fundamental  information  on  two  novel  combustion  processes.  Concepts 
which  represent  a  new  departure  and  extension  of  conventional  combustion 
practice  can  evolve  from  the  experimental  data  being  obtained.  Aspects  of  the 
radiative  ignition  and  enhancement  concept  have  been  demonstrated  under  static 
(no  flow)  conditions.  Successful  pulsed  light  source  ignition  experiments 
confirm  the  radiative  augmented  concept  and  demonstrate  the  technical 
feasibility  of  designing  an  advanced  optical -radiative  igniter.  Successful 
ignition  with  a  continuous  light  source  implies  the  possibility  of  using  the 
light  as  an  optical  radi ati ve  f 1 ame  stabi 1 i zer  (with  zero  pressure  drop) 
instead  of  the  conventional  f landholders.  The  ability  to  use  laser  radiative 
sources  further  broadens  the  future  potential  of  the  photoaugmentation 
technique.  Rather  than  point  ignition  or  local  enhancement,  line  ignition  or 
area  ignition  (by  using  optical  scanning)  may  be  possible  with  lasers.  TTTe 
results  of  combustion  enhancement  in  terms  of  higher  flame  propagation 
velocities  and  extended  extinction  times  and  distances  demonstrate  a  potential 
opportunity  to  extend  the  combustor  operating  limits  by  utilizing  the 
radiative  technique.  It  is  construed  that  the  enhanced  flame  propagation  can 
be  translated  into  higher  combustion  rate  and  extended  flammability  limits. 

From  the  experimental  results  reported  here  we  gain  confidence  that 
radiative  augmented  combustion  is  a  potentially  viable  technique  to  extend 
current  aircraft  operating  limits.  Eventual  application  to  gas  turbine  engine 
systems  is  envisioned  both  for  improved  combustor  operation  and 
flameholding.  Some  future  areas  of  potential  application  are:  high  altitude 
combustor  reignition  following  flame-out,  drag-free  flame  stabilization  in 
supersonic  combustors,  and  added  flexibility  for  conventional  combustors  to 
use  future  alternate  fuels.  To  this  end,  radiative  ignition  and  combustion 
enhancement  experiments  under  flow  conditions  are  required  as  well  as 
continued  VUV  light  source  development  In  the  direction  of  improved  beam 
optics. 


The  catalytic  flame  stabilization  concept  is  particularly  important 
to  aeropropulslon  combustion:  turbopropulsion  mainburners,  afterburners,  duct 
burners,  and  ramjet  dump  combustors.  Potential  benefits  include  improved 
ignitablllty,  stability  as  well  as  efficiency,  and  combustion  design 
flexibility  for  alternate  fuel  usage.  In  the  afterburner  application,  the 
conventional  bluff-body  flameholder  can  be  replaced  by  a  porous  catalytic 
device  resulting  In  less  pressure  drop  than  a  solid  flameholder  of  equal 
cross-sectional  area.  It  can  broaden  afterburner  stability  range  by  allowing 
for  operation  at  inlet  velocities,  temperatures,  and  fuel  mixtures  where 
conventional  flameholders  begin  to  fail.  In  addition,  it  has  the  advantage  of 
being  a  passive  autoignition  device.  Additional  catalytic  flameholding 
experiments  should  be  performed  at  higher  reactant  flow  rates  where  the 
maximum  benefit  of  the  catalytic  flameholders  can  be  realized.  Further 
development  of  high  temperature  substrate  material  is  also  required. 
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STATEMENT  OF  WORK 


The  work  statement  covering  the  contract  period  April  1981  through 
March  1983  is  presented  below.  The  contractor  shall  furnish  scientific 
effort,  together  with  all  related  services,  facilities,  supplies  and  materials 
needed  to  conduct  the  research. 

(1)  Radiative  Augmentation 

A.  Static  Experiments 

~T.  Continue  radiative  static  combustion  experiments.  Use  pulsed 
ILC  and  continuous  EIMAC  light  sources  to  determine  and  measure 
photochemical  ignition  and  combustion  enhancement. 

2.  Select  and  employ  appropriate  vacuum  ultraviolet  lasers  to 
conduct  radiative  static  experiments.  Investigate  the  dependence  of 
radiative  ignition  and  enhancement  on  light  wavelength.  Identify 
the  most  effective  wavelength  widths(s)  for  augmentation  and 
applications. 

B.  Flow  Experiments 

1.  Provide  a  combustion  test  facility  to  conduct  radiative  and 
catalytic  experiments  under  flow  conditions.  Use  advanced 
diagnostics  where  necessary  to  obtain  high  quality  experimental 
data. 

2.  Design  and  conduct  radiative  experiments  under  flow 
conditions.  Investigate  radiative  effects  on  ignition  and 
combustion  enhancement  in  premixed  flowing  reactive  mixtures 
subjected  to  pulsed  or  continuous  vacuum  ultraviolet  light.  The 
objective  is  to  attempt  to  demonstrate  radiative  augmentation  under 
flow  conditions. 

C.  Analytical  Modeling 

1.  Incorporate  the  emission  characteristics  of  the  light  sources 
used  into  the  current  combustion  initiation  model.  Compare  the 
experimental  data  with  the  model  result.  Use  the  model  to  elucidate 
the  photoabsorption-kinetics  interaction  and  to  assist  in  system 
optimization.  Include  heat  loss  mechanism,  species  diffusion  and 
wall  recombination  in  the  model.  Undertake  additional  model 
revisions  and  refinements  as  deemed  necessary. 

D.  Evaluation 

1.  Analyze  and  assess  the  experimental  and  analytical  results, 
particularly  as  they  apply  to  application  in  practical  combustion 
hardware. 

2.  Develop  a  test  plan  for  larger  scale  combustion  tests  of  the 
most  promising  radiative  augmentation  technique(s). 


(2)  Catalytic  Augmentation 

A.  Design  and  Acquistion 

1.  Design,  fabricate,  and/or  acquire  catalytic  flameholders. 
Evaluate  and  screen  chemical  and  physical  characteristics  of 
flameholders  to  determine  best  candidates  for  flow  experiments. 

B.  Flow  Experiments 

1.  Perform  experimental  flame  stabilization  tests  on  a  standard  V- 
gutter,  and  non-catalytlc  and  catalytic  monoliths  under  both  non- 
reacting  and  reacting  flow  conditions. 

2.  Perform  experimental  flame  stabilization  test  on  best 
flameholder  candidates  determined  in  Work  Statement  (2)A.l. 

C.  Analytical  Madelinq 

1.  Review  and  revise  existing  catalytic  flameholding  model  to 
obtain  a  sufficiently  detailed,  balanced  engineering  model.  Compare 
model  predictions  to  experimental  data  to  ascertain  agreement  and 
resolve  any  discrepancies  in  regimes  of  interest.  Use  the  model  to 
explore  the  Interaction  of  aerodynamic  effects  and  catalytic  and 
homogeneous  combustion.  Utilize  the  model,  as  an  adjuct  to  the 
experimental  effort  in  initially  analyzing  the  data  and  ultimately 
in  performance  optimization. 

0.  Evaluation 

1.  Analyze  and  assess  the  experimental  and  analytical  results, 
particularly  as  they  apply  to  application  in  practical  combustion 
hardware. 


2.  Oevelop  a  test  plan  for  larger  scale  combustion  tests  of  the 
most  promising  catalytic  augmentation  technique(s). 
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I.  INTRODUCTION 


The  objective  of  this  subcontract  was  to  evaluate  the  feasibility  of  using 
catalytic  honeycomb  flameholders  to  expand  flame  stabilization  limits  for  poten¬ 
tial  applications  such  as  in  aircraft  afterburners.  This  was  performed  by 
experimentally  evaluating  flame  stabilization  using  catalytic  and  noncatalytic 
flameholders  in  a  research  scale  flow  tube  combustor.  Experiments  were  con¬ 
ducted  using  premixed  propane/air  to  determine  the  effects  of  fuel/air  ratio, 
flow  velocity,  and  flameholder  type  on  the  flame  stability  limits.  Several 
other  diagnostics  were  also  used  to  gauge  the  efficiency  of  the  flame  stabiliza¬ 
tion  and  combustion  process.  These  diagnostics  included  gas  temperature  pro¬ 
files  in  the  combustor,  flameholder  temperatures,  pressure  drop  across  the 
flameholders,  and  CO  and  COj  concentrations  in  the  burned  combustion  gases. 

High  output  combustors,  such  as  in  turbojet  engines,  require  flame  stabil¬ 
ization  in  a  high  velocity  flow  of  fuel  air  mixture  within  a  duct.  Since  the 
overall  flow  velocity  is  usually  much  greater  than  the  burning  velocity  of  the 
mixture,  recirculation  regions  behind  solid  bodies  placed  in  the  flow  are  often 
used  to  create  stable  pilot  flames  which  serve  as  continuous  ignition  sources 
for  the  total  flow.  The  theory  of  flame  stabilization  due  to  recirculation 
zones  set  up  behind  bluff  bodies  has  been  very  thoroughly  treated  in  the 
literature.  Both  experimental l_s  and  theoretical*-*  discussions  of  this  pro¬ 
cess  are  available.  At  sufficiently  high  velocities  the  flame  will  no  longer 
be  stable  and  flame  blowoff  occurs.  The  conditions  at  which  this  happens  are 
dependent  on  the  flameholder  shape,  size,  pressure,  temperature,  and  fuel  air 
ratio.  Maximum  stabilization,  i.e.,  highest  blowoff  velocities,  occur  at  fuel 
air  equivalence  ratios  near  stoichiometric. 

In  the  case  of  catalytic  combustion,  the  efficiency,  product  gas  composi¬ 
tion,  and  temperature  profiles  can  be  modeled  and  compared  with  experimental 
results. However,  there  exist  no  theoretical  or  experimental  results  which 
detail  the  effect  of  porous  catalytic  flameholders  on  the  flammability  limits 
of  gas  turbine  or  research  combustors.  It  was  the  purpose  of  this  feasibility 
study  to  determine  whether  cataly tically  active  flameholders  have  wider  flame 
stability  limits  than  do  noncatalytic  flameholders,  and  to  gather  sufficient 
diagnostic  information  to  permit  a  preliminary  analysis  of  the  results. 
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II.  EXPERIMENTAL  FACILITIES 


In  this  section  the  three  major  components  of  the  combustor  facility  are 
discussed.  These  components  are  the  flow  tube  combustor,  the  f lameholders ,  and 
the  combustion  diagnostics. 

A.  FLOW  TUBE  COMBUSTOR 


The  atmospheric  pressure  flow  tube  combustor  is  shown  in  Fig.  1.  The  three 
section  brass  flow  tube  is  1.88  m  long.  The  first  section  (.0.39  m  length) 
serves  as  a  mixing  chamber  for  the  metered  fuel/air  mixture.  The  fuel  is  com¬ 
mercial  grade  propane  and  is  metered  through  a  0.13  cm  diam  calibrated  critical 
orifice.  The  air  is  delivered  from  dried  compressed  "house"  air  (rated  capacity 
of  compressor:  550  kPa  at  2.35  m*(STP)min-1)  and  metered  through  a  calibrated 
orifice  (0.68  cm  diam). 

The  fuel/air  mixture  is  radially  injected  into  a  7  cm  long,  10.2  cm  i.d. 
section  of  the  flow  tube  to  enhance  mixing  of  the  fuel  and  air.  Slightly 
downstream,  a  perforated  disc  (0.32  cm  diam  holes  with  0.48  cm  spacing)  per¬ 
forms  additional  mixing.  The  remainder  of  this  first  section  (0.32  m)  was  in¬ 
sulated  with  alumina  cylinders  (Zircar  Products  Inc.,  Type  ALC)  which  reduced 
the  cross  section  diameter  to  5.1  cm.  At  the  entrance  to  the  insulated  portion, 
a  3  cm  long  flow  straightener  fabricated  from  aluminum  honeycomb  (roughly  circu¬ 
lar  0.32  cm  diam  cells)  was  situated.  Near  the  end  of  this  first  section,  a 
pressure  port  was  installed  for  connection  to  the  reference  side  of  a  differen¬ 
tial  pressure  transducer.  Unless  otherwise  specified,  all  feedthrough  connec¬ 
tions  to  the  combustor  were  made  using  brass  1/4  in,  Swagelok  fittings. 

The  second  or  middle  section  of  the  combustor  housed  the  flameholder.  This 
0.26  m  long,  water  cooled  section  (Fig.  2)  was  fabricated  from  5.3  cm  i.d. 
brass  tubing  with  no  internal  alumina  insulation.  This  water  cooled  section 
was  used  to  reduce  the  operational  temperature  of  the  flameholders  and  to  prevent 
flame  stabilization  on  the  hot  alumina  insulation.  Several  feedthroughs  were 
provided.  Two  1/8  in.  Swagelok  feedthroughs  (1.7  cm  apart)  provided  support 
and  cooling  water  to  the  flameholder  support.  Two  Swageloks  (a  1/4  in.  and  a 
3/8  in.  fitting)  were  equipped  with  cemented  glass  windows  for  observing  the 
flamer ilder  and  the  flame  front.  These  two  ports  were  water  cooled  externally 
to  prevent  the  epoxy  cement  1  rum  degrading.  Thu  (1/4  in.)  Swagelok  port  alter¬ 
natively  served  as  a  thermocouple  wire  feedthrough  when  flameholder  temperatures 
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were  measured.  Just  downstream  of  the  flameholder  another  feedthrough  was 
equipped  with  a  ceramic  tube  insulated  spark  ignition  wire  which  could  In- 
inserted  and  withdrawn  from  the  combustion  gas  stream.  A  feedthrough, located 
8.4  cm  downstream  from  the  flameholder,  provided  a  port  for  a  thermocouple 
probe  or  pressure  measurements. 

The  final  combustor  section,  1.23  m  in  length,  was  lined  with  5.1  cm  diam 
alumina  insulation  for  its  entire  length.  Several  (1/4  in.)  Swagelok  feed¬ 
throughs  were  provided  at  15  cm  spacings  for  various  diagnostic  probes,  such  as 
thermocouples  and  combustion  gas  sampling.  The  exterior  tube  wall  was  water 
cooled. 

B.  FLAMEHOLDERS 


Three  types  of  flameholders  were  used  in  this  work.  These  cylindrical 
f landholders,  either  catalytic  or  noncatalytic  ceramic  honeycomb,  or  solid  brass 
were  held  within  a  (2.8  cm  o.d.,  2.5  cm  i.d.,  5.1  cm  long)  tubular  stainless 
steel  housing  supported  by  a  loop  of  0.32  cm  diam  stainless  steel  tubing  as 
shown  in  Fig.  3.  The  "U-tube"  support  was  silver  soldered  to  the  flameholder 
housing  and  was  water  cooled.  Swagelok  feedthroughs  supported  the  whole  flame- 
holder. 

The  ceramic  flameholders,  both  catalytic  and  noncatalytic,  were  fabricated 
from  Cordierlte  monolithic  cylindrical  shaped  honeycomb  plugs  (2.5  cm  o.d.  and 
5.1  cm  length;  Corning  No.  9475)  having  31  square  cells  per  cm1  cross  section 
and  cell  walls  0.03  cm  thick.  Cordierite*7  is  a  magnesia,  alumina,  silica 
composite  which  has  a  maximum  operating  temperature  of  %  1570  K,  a  porosity  of 
about  30%,  and  a  mean  pure  size  of  3.5  pm.  Its  thermal  conductivity  is  about 
twice  as  high  in  a  parallel  direction  to  the  cells  as  it  is  in  a  perpendicular 
direction.  Also,  it  has  a  low  (1  *  10-6  per  °C)  thermal  expansion  coefficient. 
The  catalytic  ceramic  flameholders  were  custom  prepared  by  a  vendor18 
by  repeated  inpregnation  of  Cordierite  with  an  organic  palladium  solution  and 
drying.  The  palladium  mass  loading  was  2.0  g  per  bulk  liter  of  monolith  (total 
volume  of  honeycomb) . 

In  the  experiments  reported  here,  some  of  the  honeycomb  flameholder  ceils, 
catalytic  and  noncatalytic,  were  sealed  at  each  end  using  Ceramabond  503  ad¬ 
hesive  as  shown  in  Fig.  3.  A  total  of  89  cells  remained  open  to  the  fuel/air 
gas  stream  with  68  cells  blocked.  This  configuration  was  used  to  hold  the 
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FIGURE  3  FLAMEHOLDER  WITH  HONEYCOMB  INSERT 
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honeycomb  in  the  steel  cylinder  and  to  prevent  melting  of  the  Cordierite  under 
combustion  conditions  (vide  infra) . 

C.  COMBUSTOR  DIAGNOSTIC  INSTRUMENTATION 


The  apparatus  was  equipped  with  three  combustion  diagnostics.  Temperature 
measurements  were  made  using  0.025  cm  diam  Pt/Pt-10%  Rh  (S-type)  thermocouples 
and  a  calibrated  Omega  Model  199  digital  temperature  readout.  Dual  bore  alumina 
tubing  was  used  to  insulate  the  thermocouple  leads.  The  thermocouple  junctions 
were  protected  with  a  thin  (less  than  1  pm  thick)  beryllia-yt tria  coating,  and 
were  in  direct  contact  with  the  combustion  gases.  Flameholder  temperatures 
were  measured  using  similar  thermocouples  buried  in  the  ceramic  adhesive  at 
either  end  of  the  honeycomb. 

The  pressure  drop  across  the  flaraeholders  was  measured  using  a  differential 
pressure  transducer  (Validyne  Engineering  Corp.  Model  DP15-22).  The  reference 
side  of  the  transducer  was  connected  upstream  of  the  flameholder  (see  Fig.  1) 
while  the  other  side  of  the  transducer  was  connected  to  a  port  8.4  cm  downstream 
of  the  flameholder.  The  connections  were  made  with  copper  tubing.  The  differ¬ 
ential  pressure  range  for  this  transducer  was  1.3  kPa  (10  Torr)  with  a  linearity 
of  ±  0. 52  of  full  scale  sensitivity. 

At  a  position  0.25  m  downstream  of  the  flameholder,  a  gas  sampling  port 
was  attached  to  the  flow  combustor.  The  undiluted,  continuously  sampled  combus¬ 
tion  products  were  analyzed  for  CO  and  C02  content  using  a  nondispersive  infra¬ 
red  (NDIR)  analyzer  manufactured  by  Horiba  Instruments  Inc.  (Model  MEXA-300B). 
This  instrument  measures  CO  in  the  range  of  0-7%  (volume  percent)  and  C02  in  the 
range  of  0-15%.  Reproducibility  is  quoted  as  ±  2%  of  maximum  CO  or  C02  concen¬ 
tration.  This  analyzer  was  also  tested  for  sensitivity  to  hydrocarbons.  The 
C02  analysis  was  found  to  be  weakly  sensitive  to  a  gas  stream  of  pure  methane 
(100%  CH,.  gave  an  apparent  C02  concentration  of  ca.  1%).  Both  CO  and  C02 
detectors  were  found  to  give  no  significant  response  to  pure  acetylene,  pure 
propane,  or  water  saturated  air  streams. 

These  CO  and  C02  measurements  (especially  the  CO  concentrations!  were 
expected  to  be  sensitive  to  details  of  the  flame  stabilization  process  near 
blowoff.  Preliminary  chemiluminescence  N0X  measurements  were  made  but  were 
discontinued  because  of  strong  water  interferences. 
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til-  RESULTS  AN  DISCUSSION 


The  basic  experimental  objective  of  this  work  was  to  establish  the  effect 
of  catalytic  flameholders  upon  the  flame  stability  limits  in  a  research  scale, 
flow  tube  combustor.  Several  different  flameholder  diameters  and  honeycomb  cell 
sizes  were  originally  scheduled  for  detailed  testing.  However,  it  was  determined 
experimentally  that  most  of  the  ceramic  flameholders  (only  Cordierite  honeycomb 
was  available  for  tests)  melted  before  combustion  diagnostics  could  be  performed. 
Monoliths  of  5.1  cm  length,  2.79  cm  and  1.44  cm  overall  diameters,  and  cell 
densities  of  31,  62,  and  93  cells  per  cm2  were  tested.  All  combinations  gave 
only  short  combustion  service  lifetimes.  Only  the  2.79  cm  diam  honeycombs 
having  31  cells  per  cm2  and  modified  by  sealing  off  some  of  the  outer  cells 
(Fig.  3)  gave  reasonable  lifetimes.  Therefore,  all  the  results  presented  below 
are  for  this  type  of  flameholder. 

A.  PRESSURE  DIFFERENCE  MEASUREMENTS 


One  method  of  determining  if  a  specific  flameholder  or  a  bluff  body  causes 
a  significant  mechanical  power  loss  ("insertion  loss")  is  to  measure  the  pressure 
drop  across  the  body.  It  is  reasonable  to  expect  a  larger  pressure  difference 
across  a  solid  body  than  across  a  more  porous  body,  such  as  the  ceramic  honey¬ 
combs.  Higher  gas  flow  velocities  are  expected  to  increase  the  pressure  diff¬ 
erences  across  the  bodies.  Both  of  these  effects  are  fully  manifest  in  the 
room  temperature  flow  experimental  data  obtained  for  a  solid  bluff  body,  a  fully 

open  honeycomb,  and  a  partially  blocked  honeycomb  as  shown  in  Fig.  4.  It 

should  be  noted  that  the  magnitudes  of  the  pressure  drops  are  rather  small. 

Even  for  the  highest  air  flow  velocity  (ca.  12  m  s_1),  the  drop  is  only  25  Pa 

(ca.  3  x  lO-1*  atm).  Also  note  that  the  velocities  given  are  for  300  K  gas  and 

for  flow  through  the  total  combustor  cross  section.  In  the  annular  region 
around  the  flameholder  the  velocities  are  larger  clue  to  partial  blocking  of 
the  combustor  cross  section.  [tie  solid  bluff  bodv  blocked  287,  of  the  total 
combustor  cross  section. 

It  is  of  interest  to  compare  the  experimental  results  for  the  solid  bluff 
body  with  a  theoretical  result.  Using  simple  arguments  from  fluid  mechan ics 1 9 
one  obtains  the  following  expression  based  on  conservation  of  flow  field  momen- 
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where: 


Ai 

Aj 

CD 

AP 

P 

V 


flow  tube  area  upstream  of  flameholder  (22.1  cma) 

cross  section  area  of  flameholder  (6.1  cm1) 

empirical  drag  coefficient  (-^  0.8  for 

Reynolds  numbers  of  %  10“,  Ref.  19) 

pressure  drop  (Pascals) 

muss  density  of  air  (1.2  *  IQ"3  g  cm  3 ) 

average  upstream  air  velocity  (cm  s"1) 


Using  this  expression,  we  estimate  for  the  solid  bluff  body  AP  =  11,  14,  17, 
and  20  Pa  for  8.9,  10.2,  11.2,  and  12.4  m  s" 1  average  upstream  air  velocities. 
At  10.2  m  s~‘  the  experimental  result  is  18.1  Pa  or  a  232  difference.  Thus  for 
solid  bluff  bodies,  the  pressure  drop,  AP,  can  be  approximated  using  the  above 
expression.  Probably  a  much  better  comparison  could  be  obtained  using  an 
experimental  CD-  Porous  bodies  will,  of  course,  offer  a  smaller  resistance  to 
the  flow  and  will  result  in  smaller  pressure  drops.  However,  since  drag  coef¬ 
ficients  for  honeycomb  structures  are  not  available,  a  similar  calculation  of 
pressure  drop  cannot  be  made. 


B.  blowoff  LIMITS 

In  this  program,  the  eftect  of  catalytic  honeycomb  flameholders  as  compared 
to  noncatalytic  honeycomb  and  solid  body  flameholders  were  of  interest.  Fuel- 
rich  blowoff  limits  were  found  to  be  nunreproducible.  Also  the  flames  oscillated 
between  flameholder  attachment  and  lift-off  over  a  relatively  wide  range  of 
fuel/air  ratios  causing  difficulties  in  determining  the  blowoff  limit.  There¬ 
fore,  our  efforts  were  directed  towards  measurement  of  fuel  lean  blowoff  limits 
which  were  more  reproducible,  and  are  of  more  practical  interest. 

Before  discussing  the  honeycomb  flameholder  results,  the  lean  blowoff 
limits*  for  three  solid  E lameholders  will  be  presented  (Fig-  5).  The  flow 
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velocity  refers  to  the  unburned  gas  velocity  upstream  of  the  flameholder.  These 
data  indicate  that  the  lean  blowoff  limits  are  dependent  on  the  diameter  of  the 
cylindrical  flameholder  (as  would  be  expected)  but  in  the  inverse  order  to  that 
expected.  One  anticipates  wider  blowoff  limits  (i.e.,  smaller  fuel  lean  blowoff 
equivalence  ratios)  for  the  larger  solid  body  flameholder.  The  experimental 
results  are  in  the  opposite  order — the  smallest  diameter  flameholder  shows  the 
lowest  lean  blowoff  limits.  However,  if  the  lean  limits  are  replotted  as  a 
function  of  the  (unburned)  gas  velocity  in  the  annulus  between  flameholder  and 
combustor  inner  wall,  the  results  become  independent  of  flameholder  diameter  as 
shown  in  Fig.  6. 

This  result  can  be  qualitatively  understood  in  terms  of  two  competing  pro¬ 
cesses.  The  flames  should  be  more  stable  for  larger  flameholder  diameters20 
(blowoff  velocity  «  d0'1*5)  but  the  annular  flow  velocity  increases  for  larger 
flameholders,  decreasing  the  flame  stability.  Plotting  the  data  as  a  function 
of  annular  velocity  should  show  the  true  variation  of  blowoff  velocity  with 
flameholder  diameter,  but  Fig.  6  shows  no  significant  systematic  variation. 

For  a  d0’45  dependence,  we  estimate  relative  blowoff  velocities  of  1.00,  0.91, 
and  0.66  for  3.56,  2.87  and  1.44  cm  o.d.  flameholders  (normalized  to  d  =  3.56 
cm).  The  102  difference  between  the  3.56  and  2.87  cm  diam  flameholders  would 
not  be  visible  considering  the  scatter  in  these  data.  For  the  d  *  1.44  cm 
flameholder  only  two  points  were  measured  and  the  accuracy  of  the  blowoff  limit 
is  estimated  to  be  significantly  poorer,  so  even  the  absence  of  an  observed 
trend  in  these  data  is  not  surprising.  Also  in  this  work  the  flameholders 
were  of  significant  size  relative  to  the  combustor  diameter.  In  this  case 
it  would  not  be  unexpected  to  see  that  flame  interactions  with  the  combustor 
wall  affect  the  blowoff  limits. 

To  measure  the  effect  of  a  catalytic  flameholder  upon  the  blowoff  limits, 
two  types  of  partially  blocked  Cordierite  honeycombs  were  tested.  One  type 
was  palladium  coated  whereas  the  ocher  was  uncoated.  They  were  otherwise 
identical.  When  tested  in  the  flow  tube  combustor,  the  results  shown  in  Fig. 

7  were  obtained.  Not  only  does  the  palladium  coated  flameholder  show  a  consis¬ 
tently  leaner  blowoff  limit  (smaller  £>),  but,  in  addition,  the  limiting  equiva¬ 
lence  ratios  decrease  at  higher  flow  velocities.  This  is  an  unexpected  result 
since  normally  solid  body  limits  shift  toward  stoichiometric  with  increasing 
velocity.  At  the  relatively  small  velocities  utilized  here,  this  shift  is 
expected  to  be  slight,  so  the  fact  that  the  noncatalytic  blowoff  limits  show  no 
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FIGURE  5  BLOWOFF  LIMITS  FOR  SOLID  BODY  FLAMEHOLDERS 
O  =  3.56  cm  dlam;  Q  **  2.87  cm  diatn;  A  =  1.44  cm  diam. 


FICURE  6  BLOWOFF  LIMITS  FOR  SOLID  BODY  FLAMEHOLDERS 

O  *  3.  56  cm  diam;  □  =  -’-H7  cm  diam;  A  =  1.44  cm  diam. 
Gas  velocity  for  annular  region  between  flameholder 
and  cumlustor  wall. 
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UPSTREAM  QAS  VELOCITY,  m  s'1 

FIGURE  7  COMPARISON  OF  BLOWOFF  LIMITS  FOR  2.79  cm  diatn 
CERAMIC  HONEYCOMB  FLAMEHOLDERS  AND  2.87  cm  diam  SOLID  BODY  FLAMEHOLDER 

A  =  catalytic;  O*  noncatalytic. 

Dashed  line  =  solid  body.  Data  as  in  Fig.  5. 

slope  is  not  surprising.  Both  of  these  observations  strongly  suggest  that  a 
catalytic  flameholder  should  be  more  resistant  to  "flame-out"  caused  by  fuel 
lean  blowoff  as  compared  with  a  similar  noncatalytic  flameholder. 

The  comparison  of  blowoff  limits  from  a  solid  body  flameholder  and  honey¬ 
comb  flameholders  of  nearly  the  same  diameter  is  not  strictly  valid  since  the 
honeycomb  exhibits  a  lower  pressure  drop  for  a  given  gas  velocity.  Tf  the  limits 
were  identical,  the  honeycomb  flameholder  would  be  more  efficient  (smaller  pres¬ 
sure  drop).  Data  at  higher  velocities  would  be  useful  in  clarifying  this  compari¬ 
son,  but  the  available  apparatus  had  a  limited  velocity  span. 

C.  COMBUSTOR  TEMPERATURE  PROFILES 

having  established  the  fuel  lean  blowoff  limits  for  the  catalytic  and 
noncitalytic  flameholders,  several  parameters  were  measured  under  this  range 
of  conditions.  The  first  parameter  measured  was  the  axial  temperature  profile 
along  the  flow  tube  combustor  downstream  of  the  flameholder.  Because  the 
coated  and  uncoated  flameholder  blowoff  limits  are  different,  the  coated  flame¬ 
holder  was  operated  at  both  its  blowoff  limit  and  at  the  blowoff  limit  of  the 
uncoated  flameholder  in  order  to  make  comparisons.  Note  that  it  is  not  possible 
to  operate  the  noncatalytic  f lame holder  at  both  sets  of  conditions.  Thus  there 
are  three  sets  of  data  to  compare:  (1)  the  coated  flameholder  at  its  blowoff 
limit  conditions;  (2)  the  coated  flameholder  at  the  uncoated  blowoff  limit  con¬ 
ditions;  and  (3)  the  unrooted  ll. une holder  at  its  blowoff  limit  conditions. 
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Actually  in  each  case  tht  equivalence  rat  was  0.02  unit  fuel  richer  than  the 
reported  limit.  This  results  in  extremely  stable  combustor  operation. 

In  Figs.  8  through  11,  the  temperature  profiles  for  four  (unburned)  flow 
velocities  tested  are  shown  The  origin  of  the  x  axis  is  placed  at  the  upstream 
face  of  the  flameholder  so  the  first  two  points  are  upstream  and  downstream 
flameholder  temperatures.  These  flameholder  temperatures  were  measured  by 
embedding  thermocouples  in  the  ceramic  adhesive  at  each  end  of  the  honeycomb, 
near  the  edge.  At  8.4  cm,  25  cm,  and  40  cm  downstream  of  the  origin,  gas  tem¬ 
peratures  were  measured.  No  corrections  were  made  for  radiation/convection 
cooling  of  the  thermocouple  data  for  either  gas  or  flameholder  temperatures. 
Since  the  insulated  combustor  walls  were  hot,  radiative  losses  were  much  lower 
than  in,  for  example,  an  unenclosed  flame  where  the  surroundings  are  cool. 

Thus  the  reported  gas  temperatures  are  probably  slightly  lower  than  actual  by 
approximately  50  K. 

With  the  catalytic  and  noncatalytic  honeycomb  flameholders  operated  under 
identical  conditions  (4)  =  0.89),  the  temperature  profiles  at  low  velocities 
(8.9  and  10.2  ms-1,  Figs.  8  and  9)  are  essentially  identical,  both  in  the  gas 
phase  and  in  the  ceramic  flameholders.  Both  flameholders  ran  quite  hot,  T  > 
1000  K,  and  were  hotter  at  the  downstream  end  than  at  the  upstream  end.  At  the 
two  higher  velocities.  Figs.  10  and  11,  the  gas  phase  temperatures  remain  nearly 
identical  but  the  flameholder  temperatures  are  observed  to  be  different  under 
these  conditions.  The  noncatalytic  honeycomb  is  much  cooler,  with  an  upstream 
temperature  of  about  700  K  compared  to  a  catalytic  flameholder  temperature 
which  is  above  1000  K. 

The  temperature  rise  observed  between  the  upstream  and  downstream 
end  of  the  ceramic  flameholders  is  shown  in  Fig.  12.  These  data  show  that  at 
1  =  0.89  both  the  catalytic  and  noncatalytic  flameholders  exhibit  similar  be¬ 
havior  with  a  temperature  rise  of  about  300  K  at  velocities  greater  than  10  m 
s_t  and  a  temperature  rise  of  only  about  100  K  at  the  lowest  velocity.  So  we 
see  that  at  the  higher  flow  velocities,  the  catalytic  flameholder  runs  much 
hotter  but  with  approximately  the  same  difference  in  temperature  as  the  non¬ 
catalytic  flameholder. 

When  the  catalytic  flameholder  was  operated  at  conditions  near  its  fuel 
lean  blowoff,  4>  =  0.7  to  0.8,  significantly  different  results  were  obtained. 
Figures  8  through  11  show  that  the  final  gas  phase  flame  temperature  is  ca. 

200  K  lower  than  for  the  other  conditions.  Calculated  adiabatic  flame  tem¬ 
peratures  with  the  two  sets  of  conditions,  <J>  *  0.75  and  <J>  =  0.89  are  1960  K 
and  2180  KJ1  compared  with  the  experimental  data  taken  at  25  cm  and  10.2  m  s- 1 
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FIGURE  8  COMBUSTOR  TEMPERATURE  PROFILES  USING  HONEYCOMB  FLAMEHOLDERS 


Upstream  gas 
-  0.89;  D 
holder,  $  = 
5. 1  cm  Ion  3 


velocity  »  8.9  ra  s*‘-  O  =  noncatalytic  flameholder,  0 
«  catalytic  flameholder,  <J>  *  0.89;  A  =  catalytic  flame- 
0.-78.  The  distance  origin  is  at  the  upstream  end  of  the 
flameholder  so  the  first  points  on  each  curve  are  honey¬ 
comb  temperatures. 
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FIGURE  9  COMBUSTOR  TEMPERATURE  PROFILES  USING  HONEYCOMB  FLAMEHOLDERS 

Upstream  gas  velocity  =  10.2  m  s-1.  O  =  noncatalytic  flameholder,  1 
*  0.89;  O  =  catalytic  flameholder,  f-  0.89;  A  =  cataLvtic  flame- 
holder,  *)  =  0.76.  The  distance  origin  is  at  the  upstream  end  of  the 
5.1  cm  long  flameholder  so  the  first  points  on  each  curve  are  honey¬ 
comb  temperatures. 
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FIGURE  10  COMBUSTOR  TEMPERATURE  PROFILES  USING  HONEYCOMB  FLAMEHOLDERS 

Upstream  gas  velocity  =  11.2  m  s~l.  O  =  noncataly t  ic  f  lameholder ,  <J> 
=  0.88;  □  *  catalytic  f lameho ider ,  *  0.88;  &  *  catalytic  flame- 

holder,  <J>  =  0.75.  The  distance  origin  is  at  the  upstream  end  of  the 
5.1  cm  long  flameholder  so  the  first  two  points  on  each  curve  are 

honeycomb  temperatures. 
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FIGURE  LI  COMBUSTOR  TEMPERATURE  PROFILES  USING  HONEYCOMB  FLAMEHOLDERS 

pa  cream  ;as  velocity  =  12.4  m  s~  1 .  O  =  noucatalvtic  tlameholder,  ; 
0.88;  Q  =  catalytic  flameholder,  t  -  0.88;  A  =  catalvtic  flame- 
holder,  5  =  0.71.  The  distance  origin  is  at  the  upstream  end  of  the 
5.  1  cm  long  !  I aiin. ! in  1  di‘ r  so  tin'  i  irsi  iwo  point  s  on  each  curve  arc 

honeycomb  temperatures. 
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FIGURE  12  TEMPERATURE  RISE  ACROSS  FLAMEHOLDERS . 

O  =  noncatalytic  flameholder;  Q  =  catalytic 
flameholder ,  $  =  0.88-0.89;  A  =  catalytic  flame- 
holder,  !  =  0.71-0.78. 

which  show  the  experimental  temperatures  to  be  about  250  K  below  adiabatic  for 
both  sets  of  conditions.  Thus  the  difference  in  the  gas  phase  temperature  mea¬ 
surements  in  Figs.  8  through  11  for  the  two  sets  of  conditions  is  simply  due 
to  the  difference  in  fuel/air  equivalence  ratios.  Figure  12,  however,  shows 
that  the  temperature  rise  in  the  flameholder  is  much  larger  under  the  more 
fuel/ lean  conditions  even  though  the  upstream  temperature  is  lower  in  every 
case.  This  is  interpreted  to  mean  that  a  larger  amount  of  catalytic  reaction 
is  occurring  at  <t>  =  0.7  to  0.8  than  at  tj>  *  0.89. 

D.  C0/C02  GAS  ANALYSIS 

Combustion  of  fuel  lean  propane/air  premixed  flames  at  atmospheric  pressure 
is  expected,  thermodynamically,  to  almost  quantitatively  convert  the  luel  to 
carbon  dioxide  and  water.  I’sing  the  data  from  Ref.  21,  the  equilibrium  concen¬ 
trations  of  GO j  and  CO  arc  plotted  in  Kin.  I*.  I  ho  C02  ooncontrat ion  varies  from 
6-101,  essential!/  due  to  dilution  by  excess  air,  while  the  CO  concentration 
only  exceeds  J . 17'  at  equivalence  ratios  greater  than  0.8.  The  experimental 
results  in  this  work  ire  given  in  Figs.  It  anti  15.  One  difficulty  with  these 
concent  rat  ion  measurements  is  the  large  quantity  ot  water  vapor  in  the  exhaust 
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FIGURE  13  CALCULATED  EQUILIBRIUM  CONCENTRATIONS  OF  CARBON 
MONOXIDE  AND  CARBON  DIOXIDE  FOR  PROPANE  AND  AIR  COMBUSTION 
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FIGURE  15  CARBON  MONOXIDE  CONCENTRATIONS  FOR  HONEYCOMB  FLAMEHOLDERS 

The  noncatalytic  (O)  and  catalytic  (Q)  flameholder  data  for  non- 
catalytic  flameholder  biowoff  limits  (0  %  0.9);  the  catalytic  (A) 
flameholder  data  at  catalytic  flameholder  biowoff  limits  ($  %  0.7 

to  0.8). 


22 


TP-432 


gases,  about  15%,  which  was  assumed  to  have  condensed  in  the  metal  lines  and 
filter  before  the  C0/C02  concentrations  were  measured  in  the  ND1R  instrument. 

This  introduces  additional  uncertainty  in  the  absolute  value  of  these  concen¬ 
trations.  A  better  technique  would  have  been  to  keep  the  lines,  filter,  and 
instrument  above  the  dew  point,  but  this  was  not  done.  The  noncatalytic  honey¬ 
comb  flameholder  results  show  the  highest  C02  concentration  and  the  lowest  CO 
concentration.  The  reason  for  the  higher  efficiency,  i.e.,  higher  CO  to  C02 
conversion  is  not  known.  The  catalytic  flameholder  under  the  same  conditions, 

<j>  =  0.89,  shows  roughly  similar  results  with  only  slightly  higher  CO  concentra¬ 
tions  at  velocities  up  to  11  m  s~ 1 ,  but  with  drastically  different  results  at 
the  highest  flow  velocity  where  the  C02  concentration  decreases  sharply  with 
a  corresponding  increase  in  CO  concentration.  Unfortunately,  when  the  sum  of 
the  CO  and  C02  concentrations  under  these  conditions,  <}>  =  0.89,  is  compared 
with  the  equilibrium  concentration.  Fig.  13,  a  larger  yield  is  observed  than 
theory  predicts.  For  example,  at  10.2  ms-1,  the  CO  +  C02  concentration  is 
observed  to  be  about  11%,  slightly  larger  than  the  calculated  10.2%.  This  is 
considered  indicative  of  the  uncertainties  in  the  overall  calibration  of  the 
NDIR  instrument  and  in  the  correction  for  the  change  of  gas  volume  due  to  water 
condensation.  No  reasonable  amount  of  smaLl  hydrocarbon  would  have  a  signifi¬ 
cant  effect.  Under  the  leaner  conditions  corresponding  to  blowoff  limits  for 
the  catalytic  flameholder,  C02  concentrations  are  observed  to  be  considerably 
smaller  and  the  CO  concentrations  somewhat  larger,  approximately  1%.  These 
data  were  obtained  between  $  =  0.71  and  0.78  where  the  calculated  CO  concentra¬ 
tion  is  negligible.  The  sum  of  the  CO  and  C02  experimental  measurements  at 

-  0.74  and  a  velocity  of  10.2  m  s- 1  gives  an  approximate  value  of  7%  (compared 
to  the  calculated  value  of  9%).  No  experimental  reason  for  this  difference  was 
found  and  it  may  indicate  the  presence  of  unburned  hydrocarbons  under  these 
condit ions . 

In  order  to  interpret  the  experimental  results,  we  have  attempted  to  remove 
the  discrepancy  between  the  experimental  and  theoretical  predictions  by  normal¬ 
izing  the  experimental  CO  concentrations.  Figure  16  shows  the  relative  fraction 
of  CO,  normalized  to  the  sum  of  CO  and  C02,  vs.  the  flow  velocity.  These  data  again 
show  that  the  catalytic  and  noncatalytic  flameholders  at  it  =  0.89  are  similar 
except  at  the  highest  velocities,  where  the  catalytic  flameholder  near  its 
blowoff  limits  yields  much  more  CO.  This  can  be  qualitatively  understood  if  we 
assume  that  the  catalytic  honeycomb  produces  a  Larger  quantity  of  CO  in  its 
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FIGURE  16  RELATIVE  CARBON  MONOXIDE  FRACTION  FOR  HONEYCOMB  FLAMEHOLDERS 

The  noncataly tie  (O)  and  catalytic  (□)  flameholder  data  for  noncatalyt- 
ic  flameholder  blowoff  limits  (4>  ^  0.9);  the  catalytic  (A)  flameholder 
data  at  catalytic  flameholder  blowoff  limits  (<f>  ^  0.7  to  0.8). 

exhaust  gases  at  ^  =  0.75  compared  to  <p  =  0.89.  This  carbon  monoxide  would  mix 
with  the  unburned  propane  and  air  to  yield  a  CO/CjHg/Oj  mixture  which  might  be 
less  efficient  towards  conversion  (i.e.,  combustion)  to  C02  than  a  simple 
C 3 H e / 0 j  mixture.  The  temperature  rise  data  in  Fig.  12  support  the  assumption 
that  much  more  catalytic  activity,  i.e.,  CO  production  occurs  under  the  leaner 
conditions.  Also  supporting  this  thesis  is  the  fact  that  CO/air  mixtures  support 
a  stable  flame  only  for  CO  concentrations  greater  than  1 2 . 5% 2 2  compared  to 
propane/air  mixtures  which  are  stable  for  propane  concentrations  down  to  2.2%. 
Therefore,  it  seems  reasonable  to  assert  that  large  quantities  of  CO  are  pro¬ 
duced  earlier  in  the  combustor  by  the  catalytic  flameholder  and  that  it  wouLd 
take  longer,  or  different  operating  conditions,  for  complete  conversion  of  CO 


to  CO 
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IV.  SUMMARY  AND  RECOMMENDATIONS 

The  objective  of  this  program  was  to  investigate  the  feasibility  of  com¬ 
bustion  enhancement  in  the  form  of  increased  flame  stability  by  the  use  of 
catalytically  active  f lameholders.  The  results  of  this  work  were  obtained  in 
an  atmospheric  pressure,  5.1  cm  diam,  research  combustor  and  were  not  intended 
to  simulate  actual  operating  conditions  of  any  particular  device.  The  results 
show  conclusively  that  the  fuel  lean  blowoff  limits  are  wider  with  a  catalytic 
flameholder  than  with  noncatalytic  but  otherwise  identical  f lameholders.  The 
combustion  efficiencies,  i.e.,  conversion  of  C0a,  for  both  catalytic  and  non¬ 
catalytic  flameholders  when  operated  under  identical  conditions  are  similar 
except  at  the  highest  velocities.  For  the  catalytic  flameholder  at  the  extended 
(leaner)  conditions,  however,  the  combustion  efficiency  appears  to  be  lower 
than  for  either  flameholder  at  higher  equivalence  ratios,  i.e.,  the  ratio 
CO/(COa  +  CO)  is  higher.  A  larger  amount  of  heat  release  is  observed  in  the 
catalytic  flameholder  under  the  leaner  conditions  and  the  resulting  gas  phase 
mixture  probably  requires  a  longer  residence  time  for  complete  combustion  than 
for  noncatalytic  flameholders. 

The  amount  of  catalyst  required  to  achieve  significant  widening  of  the 
blowoff  limits  is  small;  therefore  larger  catalyst  loading  or  more  efficient 
catalysts  may  widen  the  blowoff  limits  even  further  although  this  hypothesis 
was  not  investigated  in  the  program. 

Primary  practical  disadvantages  of  the  catalytic  flameholders  used  in  this 
feasibility  study  appear  to  be  associated  with  the  ceramic  substrates,  which  in 
the  case  of  Cordierite,  restrict  operation  to  rather  low  temperatures.  In  the 
case  of  some  other  ceramics  suitable  for  higher  temperatures,  the  substrates 
are  susceptible  to  fracture  due  to  thermal  shock  and  are  less  porous.  The 
palladium  catalyst  investigated  in  this  program  is  also  probably  not  practical 
in  view  of  its  significant  vapor  pressure  at  temperatures  above  1500  K. 

The  findings  of  this  program  show  that  catalytically  coated  ceramic  (or 
other  material)  f lameholders  may  be  a  promising  approach  to  enhanced  stabiliza¬ 
tion  of  highly  diluted  purtiallv  combusted  gas  streams  such  as  those  exiting 
the  main  burner  section  of  a  gas  turbine  engine.  The  conditions  necessary  for 
efficient  combustion  employing  catalytic  flameholders  must,  however ,  he  more 
thoroughly  investigated  before  practical  catalytic  combustors  can  ho  manufactured. 


In  view  of  the  progress  achieved  on  chis  program  and  the  potential  advan¬ 
tages  of  catalytic  enhanced  flame  stabilizations,  the  following  recommendations 
are  made: 

1.  A  similar  program  should  be  performed  to  compare  catalytic  and  non- 
catalytic  honeycomb  flameholders  in  a  larger  diameter  combustor  and 
at  higher  velocities.  The  larger  combustor  would  reduce  wall  inter¬ 
actions  and  a  wider  range  of  velocities  would  better  define  the 
variation  of  blowoff  limits  with  velocity. 

2.  Further  studies  of  the  CO,  C02,  and  hydrocarbon  emissions  should  be 
performed  along  with  (1)  above.  If  the  catalyst  produces  enough  CO 
to  significantly  reduce  the  burning  velocity  of  the  mixture,  the 
advantages  of  the  catalytic  flameholder  may  be  reduced. 

3.  Catalytic  flameholder  substrates  other  than  the  Cordierite  used  in 
this  work  should  be  tested.  These  substrates  would  include  other 
coated  ceramics,  coated  metals,  and  ceramics  incorporating  the 
catalysts  into  their  lattice. 

4.  Catalytically  active  coatings  other  than  palladium  should  be  tested. 

5.  Tests  should  be  performed  under  conditions  of  pressure,  velocity, 
inlet  temperature,  and  fuel  air  ratio  which  more  closely  simulate  the 
practical  combustors  of  interest,  and  jet  fuels  should  be  used  for 
the  tests. 

The  above  recommendations  define  the  direction  of  future  work  but  are  not 
meant  to  be  comprehensive.  Whether  this  technique  can  eventually  be  practical 
is  not  clear  at  this  time  but  an  evaluation  can  probably  be  made  after  some 
of  the  recommended  additional  experiments  have  been  performed. 
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Case  Studies  in  the  Simulation  of 
Novel  Combustion  Techniques' 

A.  E.  Cerkanowicz  and  J.  G.  Stevens 


This  paper,  along  with  the  one  provided  In  Appendix  III  presents 
details  on  the  modeling  of  catalytic  combustion  processes  within  porous 
honeycomb  catalysts  through  which  a  reactant  mixture  is  passed.  Since  the 
paper  in  Appendix  III  had  a  very  limited  outside  distribution,  both  papers 
have  been  included  to  provide  a  complete  set  of  references  on  this  modeling. 


CASE  STUOIES  IN  THE  SIMULATION  Of  NOVEL  COMBUSTION  TECHNIQUES' 


A.  E.  Cerkanowlcz 
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John  6.  Stevens 

Deportment  of  Mathematics  and  Computer  Science 
Nontclatr  State  College,  Upper  Montclair,  New  Jersey  07043 


ABSTRACT 

Digital  simulation  Is  an  area  of  considerable  activity 
and  growing  Importance  In  combustion  research.  This  paper 
presents  considerations  Involved  In  the  simulation  of  two 
novel  combustion  augmentation  techniques  —  catalytic  and 
radiative  (photochemical)  enhanced  combustion.  Both 
simulations  are  documented  as  to  the  Interrelationships 
between  the  physical  and  chemical  bases  and  assumptions  used 
In  deriving  the  model  equations,  the  program  objectives,  and 
the  Influence  of  these  objectives  on  the  approaches  utilised 
and  the  results  obtained.  The  goals  of  the  simulations 
differed  from  their  inception  and  the  Implications  of  these 
differences  are  explored.  In  our  approach,  catalytic 
combustor  operation  Is  described  by  a  two-point  boundary 
value  problem  for  a  non-linear  systma,  which  Is  solved  using 
an  adaptive,  variable-order,  finite  difference  scheme.  The 
numerical  stubbornness  encountered  and  the  occurrence  and 
numerical  treatment  of  bifurcation  are  presented.  The 
radiative  augmented  combustion  Initiation  model  Is  constituted 
as  an  Initial  value  problem  for  a  large  non-linear  system  of 
differential  equations  and  Is  solved  numerically  using  a 
software  package  Implementing  the  popular  backward  differ¬ 
entiation  formulas.  The  results  of  these  simulations  are 
compared  with  the  experimental  data  available  and  show 
generally  excellent  agreement.  Specific  applications  of 
these  simulations  are  reported. 

INTRODUCTION  AMO  OBJECTIVES 

Combustion  systems  are  limited  by  combustion-associated 
phenomena  such  as  flaamablllty,  flame  propagation.  Ignition, 
and  stable  combustion  and  by  the  formation  of  combustion- 
related  pollutants.  The  mechanisms  and  processes  associated 
with  promising  techniques  for  combustion  Initiation  and 
enhancement  and  for  reduction  of  pollutant  formation  are  not 
sufficiently  well  understood.  He  shall  discuss  and  contrast 
from  a  simulation  viewpoint  two  advanced  combustion  augmenta¬ 
tion  techniques  evaluated  In  the  Government  Research 
Laboratories  of  Exxon  Research  and  Engineering  Company. 

One  approach  Involves  combustion  assistance  via  the  use 
of  an  oxidizing  catalyst.  In  such  system,  the  fuel-air 
mixture  is  passed  over  a  catalytic  surface  which  heteroge¬ 
neously  promotes  combustion.  A  number  of  potential 
applications  have  been  Identified.  Conventional  gas  turbine 
combustors  produce  unacceptable  levels  of  the  oxides  of 
nitrogen  (N0„),  due  to  the  high  temperatures  developed  In 
the  primary  combustion  zone  before  the  Injection  of  secondary 
air.  This  high  temperature  zone  Is  necessary  to  Insure 
stable  combustion.  Catalytic  combustors  can  achieve  stable 
combustion  without  these  high  peak  temperatures ,  resulting  In 
a  substantial  reduction  In  NO*  awtsslons  (1).  Fuel-air 
mixtures  which  cannot  be  combusted  stably  by  ordinary 
homogeneous  means  can  be  burned  catalytlcally  and  thereby 
the  energy  content  of  many  types  of  low  heating  value  streams 
can  be  utilized  (2).  Automotive  catalytic  converters  are  a 
much  studied  application  of  catalytic  combustion,  one  which 
differs  from  the  applications  In  which  we  are  interested  by 
the  absence  of  significant  gas  temperature  rise  In  the 
catalytic  reactor  (3,  4). 


The  operation  of  a  catalytic  combustor  Involves  several 
chemical  and  physical  processes.  The  modeling  of  such 
systems  was  undertaken  to  study  the  Interaction  of  these 
processes,  to  provide  a  more  detailed  understanding  of 
ex  per  1  menu  1  data,  to  predict  trends  In  and  levels  of 
combustor  performance  as  a  function  of  operating  characteris¬ 
tics,  and  to  assess  new  areas  of  application. 

A  second  novel  combustion  technique,  radiative  (photo¬ 
chemical)  augmented  combustion.  Is  being  Investigated  as 
part  of  an  Air  Force  sponsored  program.  In  this  approach, 
combustion  of  on  unsensitized  fuel-air  mixture  Is  Initiated 
or  enhanced  by  Irradiation  with  ultraviolet  light  which 
results  In  the  photodissociation  of  oxygen  molecules  and 
combustion  Intermediate  species  (5).  The  overall  program 
encompasses  both  experimental  and  analytical  approaches. 

The  goal  of  the  modeling  effort  is  the  elucidation  of  the 
fundamenUl  Interaction  of  photon  absorption,  subsequent 
dissociation,  and  the  chemical  kinetics  of  the  reacUnt 
mixture.  After  comparison  of  model  results  with  experlmenUl 
dau,  the  model  Is  to  be  used  to  provide  guidance  for  the 
experimental  effort  and  as  a  basis  for  the  InterpreUtlon 
and  comparison  of  experlmenUl  results.  The  model  Is 
currently  limited  to  photochemical  Initiation  of  combustion 
In  hydrogen-oxygen  mixtures.  Component  modules  of  this 
model  are  being  used  In  a  succession  of  more  complex  models 
which  will  allow  the  examlnetlon  of  enhancement  effects,  the 
results  of  which  will  be  reported  elsewhere  (&). 

THE  MODELS 

This  section  provides  a  brief  description  of  the  models 
used  In  the  simulation  of  catalytic  combustion  and  radiative 
augmented  combustion.  The  derivation  of  the  governing 
equations  for  caUlytlc  combustion  Is  Included  In  (7),  with 
references  to  earlier  <n>rk.  Only  the  underlying  assumptions 
are  sunsnarlzed,  with  emphasis  on  their  relation  to  our 
objectives. 

Catalytic  Combustion 

In  caUlytlc  combustion  a  fuel -air  mixture  Is  passed 
through  a  honeycomb  monolithic  substrate  coated  with  an 
oxidation  caUlyst,  typically,  finely  dispersed  platinum. 
Combustion  results  from  heterogeneous  reaction  at  the 
caUlytlc  surface  and  from  homogeneous  oxidation  within  the 
gas  phase.  Figure  1  exhibits  the  relationship  of  some  of 
the  physical  and  chemical  phenomena  underlying  catalytic 
combustion.  In  addition  to  the  oxidation  reactions,  these 
Include  heat  and  mass  transport  between  the  surface  and  the 
gas,  momentum  transport  within  the  gas  stream,  and  axial 
heat  conduction  In  the  substrate. 

As  previously  stated,  our  aim  is  to  relate  catalytic 
combustor  performance,  which  Is  measured  by  overall  fuel 
conversion  efficiency  (and  to  a  lesser  extent  by  conversion 
and  temperature  profiles),  to  operating  characteristics. 

The  variables  specifying  the  operating  conditions  are 
listed  In  Table  1.  In  order  to  accomplish  this  objective. 

It  Is  necessary  to  model  each  of  the  relevant  processes  at 
an  appropriate  and  balanced  level  of  detail.  Our  approach 
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Is  that  of  a  system  study  tn  which  etch  of  the  components 
(l.e.,  the  physical  phenomena)  of  the  system  is  modeled  as 
simply  as  possible  In  order  to  arrive  at  a  system  model 
which  is  tractable  and  yet  does  not  neglect  Important 
physical  and  chemical  events. 

The  type  of  mode)  employed  Is  a  plug-flow,  steady- 
state,  one-dimensional  model,  frequently  used  (n  the 
simulation  of  ehemlca'  reactors.  This  choice  Is  appropriate 
In  that  we  do  not  propose  to  study  the  transient  behavior 
of  the  system  nor  the  details  of  the  fluid  dynamics  or 
transport  phenomena.  Rather,  a  one-dimensional  model 
Incorporating  heat  and  mass  transport  coefficients  Is  In 
line  with  the  level  of  detail  desired.  This  choice  Is  also 
indicated  by  the  desire  to  simulate  conditions  leading  to 
both  laminar  and  developing  turbulent  flow. 

ceramic  or  metallic  substrate 


**  of  tlf  Catalytic  1 


T  mass  transport 
I  of  fuel  and  0? 
I  to  surface 
fuel -air  mixture  ' 
flowing  through  channel 


-  typically  0.16  to 
0.32  cm  In  diameter 


rheat  transport 
between  surface 
and  gas 

homogeneous  combustion 
may  occur  in  gas  phase 


,  heterogeneous  combustion 
—  heat  liberated  at  sur 


surface 


1.  Cmm-pfimmm  propertloa,  auatrato 
taapmrttuft  tod  final  comcentrat  Ite 
•t  tJW  coulytt  nrftce  «r« 
ueifora  tAxow«bo«t  •  raactbi  crota- 
••ctioa .  talal  diapwrmton  ia 
noglwtad. 

2.  turn l  amidst tm  cut  occur  at  the 
eatalyat  martacm  and  In  Uia  gaa 
pltad*.  Clobal  lonvatap)  brrtkMli** 
(■t«  uprtHioma  arm  aaployad  for 
caapltti  oaidatioi  reactions  an  each 
pftaao.  rao  notoroooaioowa  react  a  or 
li  talma  to  be  first  ordor  in  foot 
ao4  toro  ordor  Is  osyoso. 

1  ■  Bsst  snd  SIS!  ttaaapen  are  oodslod 

diwf  tronofor  eooffieionts  wnoro 
tns  dxiviay  foress  sro  tbs 

is^siatwi  sod  cwnsatmiqa 
Ilf fsroncss  beta—  tko  beaU  fso 
pftsso  snd  tb«  catalytic  awrfscs. 

*•  TriMpsn  or  tbo  fsol  ia  coosldorod 
to  occur  toy  diffuaion  through  tb* 
bo»4*ry  layer  snd  to  M  OoiOMOd 
by  tin  catalytic  raoctiom  rsto  at 
tbo  a or f sco. 

S.  Tho  solocity  is  allawad  to  vary 

axially,  auCjoct  to  s  - - 

balance  wttieh  includoo  tbo  uoo  of 
•  ronciof  friction  factor . 


*•  Mill  boot  oondoctton  in  tbo 
substrate  Is  incited  is  tbo 


’■  ui  tbs  boat  released  is  tbo 
roocto*  lo  a a ana  ad  to  salt  tb 
CtsSaator  in  tbo  fos  men. 


fbe  uaesun  of  pltifflaa  la  disotoood 
is  tbo  tost.  Is  tbs  prosonco  of  laryo 
drsdisscs,  omtsi  dioporoiss  cm  effect 
tbo  atais  of  tbo  prof  lias. 


Tbo  aoo  of  yJosol  klootlcs  oo»  pmi d 
saofwl  in  oadoUny  cnnatien  oyatan* 

(|i  and  is  consist«it  «itt  tbo  ovorsil 
tml  of  dossil  and  tbo  stats  of 
Uoaladff  of  catalytic  machaniaam  for 
coa^ios  fools-  Tbo  aeaiaotiom 
coo cars lay  react lam  ordor  is  rossdwsSlo 
for  fool-loos  cosdltlsBU  of  lotorost . 

Tbis  ssouspcian  allow*  tbo  ttoo  of  o 
ono^dibsssicnol  oodal  and  po emits 
troocsMnt  of  both  lunar  and  turbulent 
flos  by  tbo  cbolos  of  opyroprma 
sorrolotloso  for  tbo  transport 

Tblo  laiopeio*  permits  tbo  seedy  of 
ray  loos  is  obltb  tbs  catalytic 
oososrsidn  css  b#  trasoport  or 
klsstically  limit ad. 


fori  lor  eaocb  ottos  aamnal  cosstaat 
So  loci  ty.  tbs  tanpor  stars  noo  ottbls 
tbo  rooctor  ibcrossoo  tbs  yos  volocity. 
thus  docroosiny  tbs  rssidoncs  tibo  and 
tbstsby  oiyslf Icobtly  affoctiny  tbo 
eososrsiM  prsdlctod  111 . 

Tbio  is  ispnrtont  os  o  (sodbsck 
•Mhaala  asd  far  eorroct  conpon  -  rr 
Of  corislc  and  aotolllc  eiSairatea. 

bo*  wo  ayata  cost  lyurot  ions .  looa 
of  od lob otic  ops ratios  by  rwlloat  boot 
loss  f  ns  tbo  easoator  ess  play  ss 
l^irtMt  role  (101  . 


catalytic 

surface 


axial  heat  conduction 
In  the  substrate 


Figure  1 

Chemical  and  Physical  Processas 
Modeled  In  Catalytic  Combustion 


Table  1 

Tarlebl es  Specifying  Operating  Conditions 


Fuel  Composition 
Equivalence  Ratio 
Heat  of  Combustion 
Dlffuslvlty  of  Futl 
Inlet  Gas  Temparature 
Inlet  Pressure 
Inlet  Velocity 


Catalytic  Pra-expontntlal 
Factor  and  Activation  Energy 
Homogeneous  Pre-exponential 
Factor  and  Activation  Energy 
Substrate  Conductivity 
Channel  Diameter 
length  of  Bed 
Foroslty 


The  assumptions  underlying  the  model  aro  summarized  In 
Table  2.  together  with  comments  relating  then  to  our 
objectives.  The  dimensional  form  of  the  governing  equations, 
boundary  conditions,  and  nomenclature  are  presented  in 
Table  3.  Simplification  of  the  equations  leads  to  a  two- 
point  boundary  value  probloa  for  a  non-linear  systen  In  fual 
flux,  substrate  temperature,  gas  temperature  and  velocity. 

The  equations  were  set  In  dimensionless  form  before 
solution.  The  dimensionless  groupings  (e.g.,  the  retlo 
of  characteristic  residence  time  to  characteristic  catalytic 
reaction  time)  were  studied  to  Increase  our  understanding  of 
the  Interaction  of  the  physical  processes  Involved  (7). 


Equations  and  Nomenclature  for 
the  Catalytic  Combustion  Model 


(1)  Fuel-species  balance  (gas-phase) 

«(*cf)  /  Elh\  «f  "Q 

s  -jr  ♦  «h  «»(-  *r)  <cf>  V  4  (cf  •  V  * 0 

(2)  Fuel-Species  balance  (surface) 
l(C,-Cfs)  .Ac..p(-^)  Cfl 

(3)  Gas-phase  energy  balance 

-SpvC,,  g  ♦  oh(Tt-T)  ♦  S(-4H)Ah  exp  J-  (Cf jV^)"0  -  0 

(A)  Solid-phase  energy  balance 

V,  -  oh(Ts  -  T> +  exp(-  srj)  ■ 0 


s‘;‘  37 


oh(T$  -  T)  +  o(-AH)Ac 


(5)  Momentum  equation 

-  <">  57  -  SI '  S  <*>  v  •  ° 

(6)  Continuity  equation 

57  (or)  •  0 

(7)  Equation  of  state 
p  OT  °ovoR  T 

p  ■  #  RT  *  “V”  V 

(8)  Stoichiometry 

[Jig;  V  i  'I. I 

(9)  Boundary  conditions 

at  z  •  01  cf  .  (Cf)0.  c0z  -  (C02)o.  v  -  v  t  -  V  . 


at  Z  *  l: 


Til  hit  3  (continued) 


Nomenclature 

A  •  pre-exponential  factor; 

Ac.  m/s;  A|,.  (mol/m3)1-n/s, 
where  n  Is  total  homogeneous 
reaction  order 

C  •  concentration  (associated 
with  catalyst  layer  If 
subscripted  by  s;  otherwise, 
gas-phase),  mol/m3 

Cp  *  specific  heat  capacity  of 
gas  at  constant  pressure, 
J/(kg-K) 

E,  •  activation  energy,  J/aol 

f  •  fanning  friction  factor 

h  •  heat  transfer  coefficient, 
J/(mZ.s-IC) 

(-OH)  •  heat  of  reaction,  J/aol 

k  *  mass  transfer  coefficient, 
a/s 

H  •  molecular  weight,  kg/aol 

n  *  reaction  order 


Radiative  Initiation 


R  •  gas  constant,  J/(mol-K) 
S  •  void  cross  section,  a2 

S$  ■  solid  cross  section, 

,2 

T  •  temperature,  t 
u  •  velocity,  a/s 
t  •  axial  dlaenslon,  a 

X  •  thermal  conductivity, 
J/(m-s-K) 

p  ■  density,  kg/a3 
o  •  wetted  perlaeter,  m 
i  *  equivalence  ratio 

Subscripts 
0  •  Inlet 
c  ■  catalytic 
f  •  fuel 
h  •  homogeneous 
02  or  0  •  oxygen 
s  •  surface  (solid) 


Self-sustained 
combustion  region 


Based  on  qualitative  and  quantitatlva  analysis  of  and 
subsequent  agreement  with  experimental  results,  a  mechanism 
for  the  photochemical  Initiation  of  sustained  coabustlon  in 
unsensltlxed  fuel-air  systems  has  been  developed  (5).  The 
Initiation  mechanlsa  follows  a  path  In  which  the  critical 
Intermediate  species,  atomic  oxygen.  Is  generated  by  photo- 
dissociation  of  oxygen  molecules  via  radiant  energy 
absorption  at  wavelengths  below  245  na 

02  +  hv  *  0  ♦  0. 

The  most  efficient  production  of  oxygen  atoas  occurs  at 
wavelengths  below  180  na  where  one  of  the  atoas  produced 
Is  in  the  excited  electronic  state,  0('D). 

When  a  critical  concentration  of  oxygen  atoms  is 
exceeded ,  combustion  Initiation  occurs.  Subsequent  reaction 
of  the  atomic  oxygen  with  the  fuel,  as  well  as  other  Inter¬ 
mediates  and  products,  leads  to  Ignition  and  sustained 
combustion.  Conceptually  Illustrated  In  the  phase-plane 
diagram  of  temperature  versus  the  logarltha  of  oxygen  atoa 
concentration.  Figure  2,  the  photochemical  Initiation 
process  follows  path  NP  (5).  Ideally,  the  photochemical 
path  need  not  Involve  a  temperature  rise.  In  reality,  some 
thermal Izatlon  does  occur  as  will  be  evident  from  results 
subsequently  presented.  The  normal  thermal  (spark)  path 
Is  shown  along  NT  and  can  be  expected  to  show  different 
Ignition  characteristics  from  the  photochemical  path.  It 
should  be  noted  that  other  combustion  Intermediates  do  play 
a  role  and  that  a  multidimensional  phase  analysis  Is 
necessary  for  a  complete  description  of  the  processes 
Involved.  Nevertheless,  oxygen  atom  concentration  provides 
a  satisfactory  surrogate  representation  of  the  general 
behavior  of  the  Intermediate  species  Involved  and  the  two- 
dimensional  analysis  gives  a  useful  framework  for  under¬ 
standing  photochemical  initiation. 

The  Initiation  model  Includes  the  reactant  mixture 
kinetics,  photodissociation  of  light  absorbing  species,  and 
temperature  rise  due  to  reaction  heat  release.  The  species 
considered  for  hydrogen-oxygen  systems  ire  listed  In 
Table  4.  It  Is  Important  to  observe  the  Inclusion  of 
several  excited  state  species,  the  effect  and  necessity  of 
which  will  be  described.  Some  ninety  chemical  reactions 
that  occur  between  these  species  ere  included  In  the  model . 
An  extensive  literature  review  was  undertaken  end  best 
values  for  rate  constants  in  the  form  AT-*exp(-C/T)  were 
chosen  for  temperature,  T,  In  the  Initiation  and  early 
Ignition  range  (300-1200K). 
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Table  4 

Hydrogen-Oxygen  Photochmalcel  System 
REACTANT  SPECIES 


H  H02 

VN>  H?02 

PHOTODISSOCIATION  REACTIONS 
°2  ♦  hv  (X  <  245)  -  0 

<*2  ♦  hv  (x  <  175)  -  0 


♦  oro) 


Oj  ♦  hv  (x  «  1140)  ♦  02  *0 

°3  ♦  hv  (l  <  310)  -  02(hg)  ♦  Of1©) 

°3  ♦  hv  (X  <  260)  *  Ojt1^,*)  ♦  0(’o) 

HjO  +  hv  (x  <  242)  -OH  +  H 

h°2  ♦  hv  (x  <  456)  »  OH  +0 

H202  ♦  hv  (X  <  365)  *  OH  ♦  OH 


Radiant  absorption  by  reactant  and  combustion  Inter¬ 
mediate  species  Is  considered  to  follow  a  Beer-Lambert  Law. 
The  dissociation  products  formed,  which  may  Include  excited 
states,  are  dependent  on  the  energy  of  the  photon  absorbed. 
The  photodissociation  reactions  which  can  occur  are  listed 
in  Table  4. 

Photochemical  dissociation  Is  handled  as  a  unlmolecvlar 
reaction  In  which  the  rate  of  production  of  dissociation 
products  Is  given  by  the  product  of  the  concentration  of  the 
given  absorbing  species  and  a  pseudo  first-order  rate 
constant,  which  depends  on  time,  the  distance  from  the 
light  source  window,  and  the  concentration  of  all  absorbing 
species  to  that  point.  The  expression  for  the  rate  constant 
Is 


I(x,t,A)T{X)eU)[A/hc]a(X)exp£-^0l(A)/j|ni(x.t)dtJ  dx 


where:  t  *  time 

x  •  distance  from  light  source  window 
X  >  wavelength 

n<  *  concentration  of  light  absorbing  species 
af  *  absorption  cross-section  of  light  absorbing 
species 

4  ■  quantum  yield  of  dissociation 
T  •  window  transmission 
I  •  source  Irradlance 

Quite  general  light  source  characteristics  can  be  considered 
In  the  specification  of  the  source  Irradlance.  Intensity  Is 
considered  to  be  dependent  on  wavelength,  time.  Input  energy, 
and  optical  arrangement  (f.e.,  parallel  beam,  focused  beam, 
point  source).  As  well,  the  transmission  characteristics  of 
window  material  such  as  quartz  or  sapphire  can  be  Incor¬ 
porated  by  specification  of  window  transailsslon.  Absorption 
coefficients  for  all  light  absorbing  species  are  tabulated, 
and  the  treatment  of  molecular  oxygen  absorption  In  the 
Schumann-Runge  bends  follows  that  given  by  G.  Kockarts  (11)- 
In  the  current  Implementation  of  the  model,  all  absorbing 
species  are  reckoned  to  dissociate  with  the  rate  constant 
calculated  from  the  cited  equation,  but  the  attenuation  of 
light  Intensity  by  the  mixture  Is  taken  to  be  due  solely  to 
02.  This  assumption  Is  reasonable,  even  In  wavelength 
regions  In  which  02  does  not  absorb  strongly,  due  to  the 
low  column  density  of  the  other  absorbers  during  the 
Irradiation  period. 

For  each  of  the  species  there  Is  thus  a  resulting 
equation  for  the  rate  of  change  In  Its  concentration 

dY  1 
~dt  ' 

_  ,.  ■  .reactant  ,  ,  .  1  .absorber  1 

1  :  ^  chem, j  'concentrations'  'photo, 'concentration 

chemical  reactions  photochemical 

Involving  1*h  species  reactions  of 

1th  jpacles 

where  r*  *  concentration  of  1th  species,  particles/cm3; 

the  +  sign  Is  used  If  the  1th  species  Is  a  product; 

-,  If  a  reactant. 

To  these  equations  1$  added  a  simple  energy  balance 
of  the  form 


and  Implemented  In  subroutine  STSSOL  (IT).  Variable  order 
Is  provided  through  deferred  corrections.  An  asymptotic 
error  estimator  Is  used  to  keep  track  of  the  rate  of 
convergence  upon  successive  use  of  deferred  corrections. 

Nhen  this  rate  Is  not  sufficiently  fast  according  to 
built-in  criteria,  the  basic  mesh  for  the  finite-difference 
discretization  Is  automatically  refined  (the  adaptive 
feature),  the  last  Iterate  Interpolated  to  the  new  mesh, 
and  the  process  repeated.  The  solution  of  the  non-linear 
system  of  equations  required  for  the  deferred  corrections 
Is  accomplished  using  Newton’s  method  and  thus  requires 
the  Jacobian  of  the  system. 

In  solving  the  catalytic  combustion  model  equations 
using  this  method,  the  choice  of  Initial  iterates  is 
critical.  For  a  given  solution  and  a  given  mesh,  there 
extsts  a  tube  about  the  solution  In  which  the  Initial 
Iterate  must  lie  In  order  to  assure  convergence  of  the 
Newtonian  Iteration.  The  size  of  this  tube,  that  is,  how 
close  the  Initial  Iterate  must  lie  to  the  solution.  Is  an 
Indication  of  the  numerical  stubbornness  of  the  system 
for  a  given  set  of  operating  conditions. 

There  are  two  sources  of  numerical  difficulty  present 
In  this  systoa;  the  first  Is  the  stiffness  often  encountered 
In  the  modeling  of  chemically  reacting  systems.  A  good 
general  discussion  of  solving  stiff  equations  Is  given  In 
(13);  specific  attention  Is  given  to  chemical  systems  In 
(!♦).  One  source  of  stiffness  Is  the  disparity  between 
heterogeneous  and  homogeneous  reaction  rates  (dependent. 

In  fact,  on  substrate  and  gas  temperatures).  To  a  lesser 
extent,  even  in  cases  of  dominant  catalytic  activity,  the 
difference  In  rates  of  change  for  the  substrate  temperature 
and  the  other  system  variables  Is  also  a  source  of  stiffness. 
In  the  course  of  solving  the  equations,  we  have  performed 
several  parametric  studies  which  Illustrate  these  points 
vividly,  which  we  shall  net  present  here. 

The  system  also  has  stiffness  associated  with  very 
saull  values  for  the  substrate  conductivity,  As.  (Hore 
properly,  it  occurs  for  large  values  of  the  Peclet  number 
for  heat  transfer,  Pes  ■  Lp0v„CpSAsSs . )  After  simpli¬ 
fication  equation  4  of  Table  3  becomes 

Vs  (^T  *  S*|  z.0)  *  SovCp(T-T0)  ♦  (-AH)  ( vCf — v0 ( Cf  )0 ) - 

Recall  that  the  boundary  conditions  for  substrate  temperature 

are 


dT  .  -  (-AH)  ,k  ,  .reactant  , 
dt  1  Aypdp  '  Chem;  'concentrations' 

all  chemical 
reactions 

where  AH  *  heat  of  reaction,  J/mol 
p  •  mixture  density,  kg/caJ 
Av  •  Avagadro’s  number,  particles/mol 
Cp  •  specific  heat  capacity  at  constant 
pressure,  J/kg-K 

The  model  Is  finally  constituted  as  an  Initial  value  problem 
for  a  large  nonlinear  system  of  differential  equations. 
Diffusion  and  convective  effects,  as  well  as  heat  loss  to 
the  surroundings,  have  been  neglected.  He  shall  subsequently 
discuss  the  use  of  this  approach  for  the  study  of  initiation 
and  the  early  stages  of  Ignition,  l.e.,  to  describe  In  this 
time  period  the  photon-reactant  mixture  Interaction  In  the 
neighborhood  of  the  light  source  window  prior  to  the 
breakdown  of  the  assumptions  made. 


SOLUTION  TECHNIQUES 


The  equations  of  the  catalytic  combustion  and  photo¬ 
chemical  Initiation  models  were  solved  numerically  using 
general  purpose  software.  Our  experiences  with  the  methods 
used  are  detailed  In  the  following  sections. 


Catalytic  Combustion 


For  the  catalytic  combustion  model,  the  boundary  value 
problem  wes  solved  using  an  adaptive,  variable-order, 
finite-difference  scheme  developed  by  lentlnl  and  Pereyra 


£»|  -SI  -0. 

dz  I  z-0  dz  1  z*L 

For  Aj  very  small,  the  solution  Is  essentially  that  given  by 
setting  the  right  hand  side  equal  to  zero.  However,  the 
boundary  conditions  are  not  now  In  general  satisfied.  Thus 
the  substrate  temperature  as  given  by  the  solution  of  the 
full  equations  has  a  component  which  changes  rapidly  toward 
the  entrance  and  the  exit  of  the  combustor.  This  development 
thus  gives  rise  to  a  "boundary  layer"  phenomenon  (a 
mathematical  term,  here  not  to  be  confused  with  a  flow 
boundary  layer)  and  the  accompanying  stiffness  (15). 

One  approach  to  overcoming  these  difficulties  is  to 
choose  an  Initial  Interate  which  lies  quite  close  to  the 
solution.  This  Initial  Iterate  can  be  obtained  by  approx¬ 
imating  the  solution  analytically.  However,  a  further 
complication  Is  that  there  are  regimes  In  which  the  system 
has  multiple  solutions.  The  source  of  this  bifurcation  Is 
Important  In  understanding  the  operation  of  a  catalytic 
combustor,  as  will  be  shown  In  the  discussion  of  the  results. 
The  numerical  Investigation  of  multiple  solutions  was 
accomplished  by  examining  the  solutions  obtained  from  three 
different  types  of  Initial  Iterates.  These  were:  an 
•unignited"  Initial  Iterate  which  represents  no  fuel 
conversion,  a  "lower-ignited"  one  obtained  from  the 
assumption  of  mass-transport  limited  catalytic  conversion, 
and  an  "upper  Ignited"  one  reflecting  significant  homogeneous 
conversion.  Generally,  for  a  given  set  of  parameters,  use 
of  these  Initial  Iterates  sufficed  etther  to  give  consider¬ 
able  assurance  that  there  was  a  unique  solution  (l.e., 
different  Initial  Iterates  would  converge  to  the  same 
solution);  or,  when  there  were  multiple  solutions,  to 
determine  easily  the  uRper  and  lower  solutions. 


On  occasion,  for  choices  of  perimeters  for  which  the 
nu*er Ice!  stubbornness  mentioned  above  was  present.  It  was 
necessary  to  use  parametric  continuation  to  obtain 
convergence  to  within  the  absolute  error  tolerance  specified, 
usually  10*3  (12,  IS).  Thus  a  sequence  of  solutions  was 
obtained,  starting  with  a  value  for  the  particular 
"offendfng"  parameter,  for  which  convergence  could  be 
obtained.  This  solution  was  then  used  as  the  Initial 
iterate  In  solving  the  equations  for  a  value  of  the  parameter 
closer  to  that  desired.  This  process  was  repeated  until  the 
solution  was  obtained  for  the  desired  pa-ameter  value. 

SYSSOL  contains  an  option  to  perform  continuation  ‘automat* 
Ically",  though  we  found  the  authors'  choice  of  only  one 
Newton  step  for  each  of  the  continuation  steps  was  not 
sufficient  to  conquer  the  most  challenging  situations. 

Parametric  continuation  was  carried  out  for  large  values 
of  the  Peclet  number  and  successful  Integration  was  obtained 
for  values  of  Pes  up  to  .lOv  However,  as  was  observed 
earlier,  for  the  bulk  of  the  reactor  length  the  substrate 
temperature  profile  and  even  more  so  the  other  profiles  agree 
with  those  obtained  for  much  lower  Pe$  values.  Thus  high 
values  were  routinely  replaced  by  a  value  known  to  give  good 
agreement  and  which  resulted  In  more  easily  achieved 
successful  Integration.  Ue  chose  this  approach  rather  than 
solve  the  locatly  adiabatic  (infinite  Pe$)  case  by  Initial 
value  methods  because  we  were  often  Interested  In  Inves¬ 
tigating  the  use  of  metallic  substrates  which  can  give  Peg 
values  below  20  and  did  not  wish  to  implement  two  different 
solution  techniques. 

Finally,  we  had  only  limited  success  using  the  adaptive 
feature  for  stubborn  problems.  That  Is,  we  observed  cases 
which  we  could  successfully  Integrate  using  an  analytical 
Initial  iterate  with  an  Initial  choice  of  101  or  201  mesh 
points,  yet  the  same  Initial  Iterate  on  26  or  51  mesh 
points  would  fall  to  converge  even  after  the  mesh  had  been 
refined  to  201  points.  The  reason  for  this  behavior  Is 
that  the  Iterate  on  the  smaller  number  of  mesh  points  was 
not  In  the  tube  of  convergence  of  the  Newton  Iteration. 

The  result  of  deferred  correction  when  Interpolated  to  the 
refined  mesh  was  thus  not  a  sufficiently  good  Initial 
Iterate  to  permit  convergence  on  the  new  mesh.  This  Is  not 
a  serious  problem  and  can,  of  course,  easily  be  overcome  by 
using  the  finer  mesh  to  begin  with.  On  problems  not 
Involving  these  numerical  difficulties,  the  adaptive  feature 
permits  efficient  solution  both  at  high  and  "engineering* 
precisions.  Indeed,  SY5S0L  has  proved  to  be  remarkably 
versatile  and  capable  of  economically  Integrating  the 
equations,  even  when  the  problem  lies  In  the  truly  "horror* 
range. 

Radiative  Auaeented  Combustion 

The  equations  of  the  photochemical  Initiation  model 
were  solved  using  the  stiff  option  Incorporating  the  back¬ 
ward  differentiation  formulas  Implemented  In  one  of  the 
Gear  family  of  codes,  STIFF,  by  Kahaner  and  Sutherland 
(17,  16).  Because  of  the  excellent  state  of  the  art  in 
codes  for  the  solution  of  Initial  value  problems  for  stiff 
ordinary  differential  equations,  such  problems  can  be 
solved  today  quite  routinely.  The  stiffness  of  the  system 
seemed  readily  apparent  to  us  at  the  outset.  For  example, 
the  oxygen  singlet  0('0)  participates  In  very  fast  reactions. 
Is  produced  only  be  photodissociation,  and  virtually  tracks 
the  shape  of  the  vacuum  ultraviolet  radiant  pulse.  Other 
species  change  on  a  much  different  timescale.  Nevertheless, 
we  did  experiment  using  the  Adams-Moulton  method  In  STIFF, 
only  to  find  that  indeed  for  one  sample  problem  the  use  of 
the  chord  method  with  an  analytic  Jacobian  (a  rather  unusual 
combination)  more  than  quadrupled  the  execution  time  compared 
to  the  stiff  option  with  the  same  solution  technique,  while 
use  of  functional  Iteration  was  not  successful  within  the 
prescribed  criteria. 

To  aid  In  the  formidable  task  of  coding  the  photo¬ 
chemical  Initiation  model  equations,  an  equation-writer 
program  from  Los  Alamos  Scientific  Laboratory  was  used, 
after  modification  to  Include  the  heat  balance  (19).  The 
program  accepts  species  and  reaction  cards  as  Input  and 
produces  subroutines  coded  In  FORTRAN  for  both  the  species 
rate  expressions  and  the  Jacobian  of  the  system.  Although 
there  exist  several  chemical  kinetics  packages,  this 


approach  lends  Itself  well  to  the  photochemical  application 
end  by  producing  the  Intermediate  source  code  provides  the 
flexibility  to  use  the  modules  already  completed  In  a 
succession  of  more  complex  models,  treating,  for  example, 
diffusion  and  wall  recombination  effects.  The  complete 
program  package  tncludes  subroutines  to  compute  the  chemical 
and  photodissociation  rate  constants. 

RESULTS :  CATALYTIC  COMBUSTION 

The  model  was  used  to  simulate  experimental  results 
obtained  by  Pfefferle.  et  al .  (20).  Solutions  to  the 
model  equations  were  obtained  using  Input  values  correspond¬ 
ing  to  the  experimental  operating  conditions.  It  Is 
Important  to  observe  that  the  model  solutions  depend  on 
Input  values  not  only  for  quantities  which  can  be  specified 
precisely  (e.g.,  reactor  length),  but  also  for  quantities 
for  which  It  Is  more  difficult  to  assign  values.  These 
Include  the  catalytic  and  homogeneous  reaction  rate 
parameters.  Since  many  studies  have  dealt  with  the 
determination  of  these  quantities.  It  would,  however, 
clearly  be  Inappropriate  to  choose  values  contrary  to  the 
bulk  of  the  experimental  results.  The  same  considerations 
apply  to  the  choice  of  the  correlations  for  the  heat  and 
mass  transport  coefficients.  Initially  we  used  correlations 
recommended  by  Hawthorne  (A).  Subsequently,  Votruba,  et  al . 
(21)  have  pointed  out  that  the  correlations  given  by 
Hawthorne  tend  to  give  somewhat  low  values  for  the  transport 
coefficients.  Preliminary  coaparfMn- ef  model  results  with 
the  experimental  data  of  Pfefferle,  et  al.  was  given  In  (7). 
He  present  here  new  results  based  on  the  correlations 
given  by  Votruba,  et  al.  and  a  previously  unreported  feature 
of  the  model  solutions. 

Figure  3  shows  the  effect  of  varying  Inlet  gas 
temperature  on  conversion  for  combustor  operating  conditions 
given  In  the  legend.  Further  experimental  details  can  be 
found  In  (20).  The  experimental  data  points  are  Indicated 
by  solid  circles,  while  the  solid  curve  was  obtained  from 
computer  solutions  of  the  model  Including  both  hoangeneous 
and  heterogeneous  contributions.  Agreement  between  the 
experimental  and  modeling  results  Is  striking. 
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Experimental  Oats  (»)  and  Nodal -Generated 
Curves  Showing  the  Effect  of  Inlet  Temperature 
_ On  Conversion _ 

For  a  given  Intet  temperature  below  390*  (point  Bt), 
there  Is  a  unique  solution  which  predicts  that  substrate 
and  gas  temperature  throughout  the  bed  equal  the  Inlet 
temperature  and  that  there  Is  no  appreciable  conversion. 

In  the  range  390-500K,  there  are  at  least  two  solutions  for 


each  Inlet  temperature .  The  upper  and  lower  solutions  are 
shown  and  undoubtedly  represent  stable  operating  states. 

The  branch  labelled  I  would  be  followed  upon  Increasing  the 
Inlet  gas  temperature.  At  a  ’light-off"  temperature  of 
approximately  490K,  catalytic  Ignition  has  occurred  and 
there  Is  a  substantial  Increase  In  conversion.  The 
extinction  branch,  labelled  E,  would  be  followed  upon 
decreasing  Inlet  temperature  after  having  obtained  catalytic 
Ignition  and  having  passed  point  82  on  the  curve.  At 
approximately  390K,  a  sharp  extinction  Is  predicted.  This 
phenomenon  was  observed  by  Rfefferle,  et  al .  (20,  p.  6) 
and  observed  and  discussed  by  Hlavacek  and  Votruba  (22). 

The  values  used  for  the  catalytic  pre-exponential 
factor  and  activation  energy  directly  affect  the  temperatures 
predicted  for  light-off  and  extinction.  Values  of  914  m/s 
for  Ag  and  48  kJ/mol  for  Eac  were  used,  which  conform  gulte 
well  to  those  cited  In  the  literature,  especially  Hawthorne 
(4). 

For  temperatures  In  the  range  390-7S0K,  the  conversion 
predicted  Is  due  almost  completely  to  the  catalytic 
reaction  with  no  significant  contribution  from  the  homo¬ 
geneous  reaction.  At  the  higher  temperatures  In  this 
range,  the  catalytic  reaction  rate  Is  limited  by  the  rate 
of  mass  transport  to  the  surface.  At  lower  temperatures, 
some  kinetic  limitation  Is  present  as  well. 

A  second  hysteresis  loop,  which  Is  associated  with  the 
homogeneous  oxidation  reaction.  Is  encountered  In  the  range 
750-350K.  Similar  considerations  apply  to  Its  branches  as 
to  those  of  the  catalytic  hysteresis  loop.  In  our  first 
comparison  with  these  data,  this  feature  was  not  reported. 

Our  model  prediction  of  this  hysteresis  was  subsequently 
confirmed  by  unreported  experimental  results  (23). 

The  curve  In  Figure  3  labelled  ’catalytic  only*  Is 
obtained  by  setting  the  homogeneous  reaction  rate  to  zero 
and  Is  thus  due  to  the  catalytic  contribution  only.  For 
Inlet  temperatures  In  this  range,  no  homogeneous  combustion 
would  occur  In  the  absence  of  the  catalytic  reaction.  For 
this  case.  In  the  absence  of  homogeneous  reaction  Incomplete 
conversion  would  occur.  Thus  the  homogeneous  combustion, 
which  completes  fuel  conversion,  Is  supported  or  "boat- 
strapped"  by  the  catalytic  reaction.  These  results  depend 
on  tne  values  for  the  homogeneous  pre-exponential  factor 
and  activation  energy.  Consistent  with  (8),  An  was  taken 
to  be  2.87  x  10’  m3/mo1-s  and  Eah  To  be  167  kJ/mol . 

The  use  of  the  model  has  revealed  other  phenomena  that 
were  unexpected  for  which  the  physical  basis  has  subsequently 
been  provided  (e.g.,  the  occurrence  of  gas  temperatures 
exceeding  the  adiabatic  flame  temperature  even  for  a  Lewis 
number  of  unity  (7)).  The  model  has  been  employed  In  many 
applications  which  cannot  be  described  In  detail  here.  We 
have  used  It  to  examine  the  potential  for  the  catalytic 
combustion  of  several  low  heating  value  streams,  to  predict 
successfully  combustor  geometry  needed  for  catalytic 
light-off  for  a  U.S.  Air  Force  afterburner  study,  and  to 
provide  Information  needed  for  the  economic  component  of 
an  environmental  assessment  of  advanced  gas  turbine 
systems  (24).  in  all,  the  model  has  proved  to  be  a  useful 
tool  In  the  Interpretation  of  experimental  data  and  for 
simulations  used  to  screen  diverse  potential  applications 
of  this  promising  technology. 

RESULTS:  RADIATIVE  AUGMENTED  COMBUSTION 

The  equations  of  the  photochemical  Initiation  model  were 
solved  for  experimental  conditions  of  Interest.  Model 
solutions  for  the  irradiation  of  a  298K,  40  kPa  stoichio¬ 
metric  hydrogen-oxygen  mixture  are  shown  In  Figure  3.  Each 
curve  represents  a  photochemical  path  for  a  given  radiant 
Intensity  in  the  temperature  versus  oxygen  atom  concentration 
phase  plane,  with  time  Implicit.  The  radiant  pulse  was 
assumed  to  be  a  critically  damped  discharge  with  time  to 
peak  intensity  of  20  us  (see  Inset  curve  of  Figure  3).  The 
radiant  energy  was  assumed  to  be  uniformly  distributed  over 
the  wavelengths  140  to  360  nm.  The  numbers  labelling  each 
curve  Indicate  the  source  Intensity,  I,  (more  properly, 
Irradlance  per  unit  wavelength)  relative  to  the  value  I)  • 
0.0687  J/cm’.s.nm  (e.g.,  85  ■  85-1)). 


For  curves  with  relative  Intensity  80  and  below,  the 
radiant  pulse  generates  oxygen  atoms  and  some  temperature 
Increase.  Oxygen  atoms  are  then  consumed  by  reaction  with 
02*  Hj,  and  other  intermediates.  The  temperature  profile 
flattens  out  and  no  rapid  temperature  rise  characteristic 
of  ignition  is  obtained.  For  I  *  80,  the  Integration  was 
carried  out  to  one  hundred  milliseconds  without  further 
appreciable  temperature  rise.  At  radiant  intensities  of 
85  and  above,  the  pulse  generates  somewhat  higher  oxygen 
atom  concentrations  than  at  the  lower  intensities.  The 
oxygen  atoms  are  again  consumed  and  the  mixture  temperature 
Increases  in  the  process.  However,  along  this  path  the 
mixture  has  entered  an  unstable  region  eventually  leading 
to  oxygen  atom  generation,  more  rapid  temperature  rise,  and 
finally  rapid  burning. 

In  Interpreting  these  results,  it  is  important  to 
remember  that  no  heat  loss  mechanism  Is  Included  in  the  model 
In  a  physical  system,  heat  transport  and  heat  loss  to  the 
surroundings  do  occur.  An  Indication  of  whether  the  mixture 
would.  In  a  real  system,  obtain  Ignition  Is  provided  by  the 
time  required  to  obtain  the  rapid  temperature  rise  character¬ 
istic  of  tgnltlon.  Thus  If  the  model  predicts  rapid 
temperature  rise  within  a  time  period  for  which  heat  loss 
wuld  not  be  sufficient  to  overcome  the  rate  of  heat 
eneratlon,  then  Ignition  In  the  physical  system  Is  Indicated 
For  our  experimental  conditions,  a  characteristic  heat  loss 
period  of  300  us  Is  reasonable.)  Equivalently,  for  the 
phase  plane  analysis  of  Figure  2  which  includes  heat  loss, 
the  mixture  has  crossed  the  separatrlx  and  entered  the 
region  ef  self-sustained  combustion. 

The  curves  shown  in  Figure  4  are  Instructive  for  the 
further  study  of  the  phase  plane.  Curve  A  exhibits  the 
effect  of  the  complete  pulse  with  relative  Intensity  100. 
Curves  8  and  C  show  the  result  of  shutting  the  light  source 
off  at  70  and  65  us,  respectively.  Both  curves  show  the 
eventual  onset  of  rapid  temperature  rise.  However,  the 
time  to  1200K  Is  135  us  for  B  and  14  000  us  for  C. 

Comparison  of  these  times  reveals  that  In  a  system  with 
heat  loss.  If  the  light  were  shut  off  at  Bs,  ignition  would 
almost  surely  still  be  obtained;  but  by  doing  so  at  Ci ,  it 
would  not.  Hence,  curve  A  must  cross  the  separatrlx  at  a 
point  between  Bi  and  Cj .  In  this  fashion,  the  Initiation 
model  can  be  used  to  determine  details  of  the  phase  plane 
for  a  given  set  of  experimental  conditions. 

The  general  phase  plane  analysis  of  Figure  2  suggests 
that  the  trajectory  for  a  mixture  which  Just  manages  to 
cross  the  separatrlx  would  then  follow  a  path  of  decreasing 
oxygen  atom  concentration  and  Increasing  temperature  to 
slightly  above  the  thermal  Ignition  temperature.  The 
trajectory  would  then  turn  around  and  proceed  to  ignition. 

The  curves  in  Figure  3  are  In  excellent  agreement  with  this 
qualitative  analysis,  as  Is  clear  from  the  comparison  of  the 
behavior  of  curves  80  and  85  and  the  known  thermal  ignition 
temperature  (point  T  on  the  temperature  scale).  These 
results  are  consistent  with  the  existence  of  a  critical 
intensity  below  which  ignition  Is  not  achieved,  in  agreement 
with  earlier  experimental  results  (5). 

In  order  to  assess  the  Importance  of  Including  excited 
state  species,  we  obtained  solutions  with  0('D)  and  Its 
reactions  Included  and  then  with  the  assumption  that  the 
photodtssoclatlon  of  O2  produced  only  ground  state  oxygen 
atoms.  The  results  differed  greatly  --  the  pulse  energy 
required  for  ignition  and  the  time  to  ignition  without 
consideration  of  0('D)  was  twice  that  for  the  full  kinetic 
set.  Hence  It  Is  clear  that  the  metastable  species 
0('D)  (and  possibly  other  excited  state  species)  and  its 
reactions  must  be  Included  In  the  modeling  of  photochemical 
Initiation. 

Final  comparison  of  modeling  and  experimental  results 
cannot  be  made  yet.  The  radiant  energy  required  to  obt-in 
ignition  for  the  40  kPa  stoichiometric  hydrogen-oxygen 
mfxture  as  indicated  by  the  modeling  results  thus  far  is 
greater  than  earlier  experimental  results  would  have  led 
us  to  expect.  However,  many  factors  remain  to  be 
Investigated.  For  example.  If  ten  percent  of  the  radiant 
energy  is  peaked  In  the  144  to  148  nm  spectral  region, 
then  a  pulse  of  the  same  total  radiant  energy  as  that  of 
the  80  curve  In  Figure  3  reaches  1200A  four  times  faster 


than  does  the  85  curve  which  Is  based  on  a  flat  spectra) 
distribution.  A  detailed  Investigation  of  the  character¬ 
istics  of  the  light  sources  being  used  Is  In  progress. 

For  the  same  conditions,  another  Important  modeling  result 
(s  that  for  critically  damped  discharges  at  con-'ant  total 
pulse  energy  there  exists  a  pronounced  :r.1mal  time  to  peak 
Intensity,  as  measured  by  the  time  required  to  reach  1200k. 

The  existence  and  determination  of  this  optimal  point  Is 
of  evident  Importance  for  light  source  design.  In  order 
to  develop  a  consistent,  detailed  description  of  photochemical 
Initiation  and  enhancement  of  combustion,  the  Investigation 
Is  continuing  with  the  complementary  experimental  and 
modeling  efforts. 

SUMMARV/CONCLUSIONS 

Both  modeling  projects  have  enhanced  our  understanding 
of  the  Interactions  of  the  physical  and  chemical  processes 
Involved  In  the  combustion  augmentation  techniques  under 
consideration.  For  catalytic  combustion,  some  of  the 
phenomena  which  result  from  these  Interactions  and  have 
been  further  elucidated  In  the  course  of  the  modeling  are 
the  effect  of  velocity  Increase  In  the  reactor  on  fuel 
conversion,  the  ’bootstrapping’  of  homogeneous  combustion 
by  the  catalytic  reaction,  and  the  catalytic  and  homogeneous 
1 Ight-off/extlnctlon  hystereses.  In  addition,  the  model 
has  proved  valuable  as  a  tool  to  predict  combustor 
performance  for  given  operating  conditions  and  has  been 
used  In  this  capacity  to  guide  experimental  efforts  and 
examine  promising  areas  of  application. 

The  results  of  the  photochemical  Initiation  modeling 
have  already  Identified  Important  considerations  for  the 
experimental  effort.  The  Importance  of  careful  spectral 
characterization  of  the  light  source  has  been  confirmed. 

The  prediction  of  the  approximate  radiant  energy  needed 
to  produce  Ignition  for  a  given  system  Is  possible  with 
the  present  model,  allowing  comparison  with  experimental 
results.  Further,  for  pulses  of  constant  total  energy,  a 
pronounced  effect  of  pulse  width  on  the  production  of 
Ignition  has  been  predicted.  This  effect,  apparently  due 
to  partial  thermal Izatlon  allowing  more  efficient  use  of 
the  oxygen  atoms  subsequently  produced  by  photodissociation, 
can  be  of  considerable  Importance  for  light  source  design. 

As  well,  the  model  has  been  used  to  study  features  of  the 
fundamental  phase  plane.  Work  Is  continuing  on  the 
application  of  the  model  In  support  of  the  experimental 
effort. 

Both  projects  have  benefited  from  the  use  of  general 
purpose  software  for  the  Integration  of  the  model  equations. 
We  anticipate  that  the  laudable  trend  In  the  development 
and  availability  of  high  quality  software  for  these 
purposes,  which  Is  easy  and  economical  to  use,  widely 
applicable,  and  reliable,  will  continue  and  even  accelerate. 
Of  course,  such  programs  still  require  Informed  use,  as 
Is  likely  to  be  the  case  for  some  time.  In  the  meantime, 
we  find  generally  applicable  the  advice  of  Shamplne  and 
Gear  (13)  that  a  conservative  choice  of  tolerance,  careful 
use  of  the  error  criterion,  experimentation,  and  a 
thoughtful  examination  of  the  numerical  results  are 
Indlspenslble. 

Finally,  we  observe  that  the  two  models  differ  In  the 
level  and  placement  of  detail  In  the  consideration  of  the 
underlying  chemical  and  physical  phenomena.  The  objective 
of  the  catalytic  combustion  modeling  project  Is  the 
description  of  steady-state  system  performance.  Because 
the  details  of  the  combustion  mechanisms  or  the  formation 
of  trace  spec'es  were  not  In  the  scope  of  the  study, 
simplified  chemical  kinetics  were  employed.  The  model 
equations  are.  however,  quite  Inclusive  In  the  balanced 
treatment  of  the  processes  Involved.  The  photochemical 
Initiation  model,  on  the  other  hand.  Incorporates  detailed 
consideration  of  photon  absorption  and  chemical  kinetics, 
but  uses  a  simplified  heat  balance.  This  approach  Is 
consistent  with  the  objective  of  using  the  model  to  study 
photon/reactant  mixture  Interaction  In  combustion  Initiation, 
In  these  examples,  the  project  objectives  have  guided  the 
model  formulations.  Thus  In  the  light  of  modeling 
objectives,  a  balance  must  be  sought,  llnneeded  detail 
results  In  costly  consumption  of  comouter  time  and  can 
hinder  the  discovery  of  significant  solution  features. 

Lack  of  appropriate  detail  can  make  modeling  results 
useless  or  even  misleading. 
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A  model  for  catalytic  combustor  operation  is  described  for  a  well-mixed  plug- flow  reactor 
Heterogeneous  and  homogeneous  reactions  are  considered  with  heat  and  mass  transfer 
handled  using  transport  coefficients.  Gas  uelocity  is  allowed  to  vary  axially  and  heat  con¬ 
duction  in  the  catalyst  bed  is  included.  Computer  solutions  were  i  btained  for  axial  pro¬ 
files  of  conversion,  gas  and  substrate  temperature,  velocity,  and  pressure.  Model  predic¬ 
tions  are  compared  to  catalytic  combustion  data,  showing  excellent  agreement.  Catalytic 
light-off,  hysteresis,  and  homogeneous  contribution  are  discussed.  Data  from  a  'hybrid 
combustor  developed  in  this  laboratory  are  analyzed  using  a  simplified  model. 


Introduction 

■4 

Conventional  gas  turbine  combustors  produce  unacceptable  levels 
of  NO,  because  of  the  high  temperatures  developed  in  the  primary 
combustion  zone,  especially  under  high-pawer  conditions.  If  the  high 
flame  temperatures  can  be  significantly  reduced  without  comprom¬ 
ising  combustion  efficiency,  then  a  substantial  reduction  in  NO, 
emissions  can  be  achieved.  Catalytic  combustion  is  being  investigated 
experimentally  and  theoretically  in  this  laboratory  as  a  means  of 
achieving  this  end,  as  well  as  improvements  in  performance  and  cost 
such  as  improved  stability  and  product  gas  uniformity.  This  paper 
describes  the  theoretical  modeling  effort. 

The  objective  of  this  effort  was  to  develop  a  sufficiently  detailed 
model  of  the  operation  of  catalytic  combustors  to  enhance  the  un¬ 
derstanding  of  this  new  technology.  CurrenUy,  the  model  considers 
steady-state  operation  of  a  well  mixed  plug  flow  reactor  with  heat  and 
mass  transfer  between  gas  and  solid  phases.  Reaction  and  heat  release 
are  considered  to  occur  both  on  the  catalyst  surface  and  in  the  gas 
phase.  By  employing  the  appropriate  conservation  equations  for  the 
two  phases,  combustor  operation  can  be  described  in  terms  of  a  two- 
point  boundary  value  problem.  A  powerful  solution  technique  is 
employed  to  obtain  computer  solutions  for  temperature  (gas  and 
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substrate),  pressure,  velocity,  and  fuel  conversion  profiles  as  a  function 
of  reactor  geometry,  catalyst  and  substrate  properties,  and  inlet 
conditions.  Reactor  configurations  containing  screens,  pellets,  or 
monoliths  can  be  handled  by  selection  of  appropriate  heat  and  mass 
transfer  correlations. 

A  discussion  of  the  numerical  technique  mentioned  in  the  foregoing 
and  of  earlier  modeling  work  in  which  the  possibility  of  parallel  gas- 
phase  reaction  was  not  considered  was  reported  previously  [l].-  This 
initial  analysis  indicated  that  velocity  variation  through  the  reactor 
(and  hence  reactant  residence  time)  cannot  be  neglected  in  calculating 
combustor  operating  characteristics.  Indeed,  constant  velocity  models 
result  in  significant  overestimation  of  conversion.  Many  previous 
modeling  efforts  neglect  the  effect  of  velocity  variation  since  they  were 
concerned  primarily  with  catalytic  clean-up  of  exhaust  or  tail  gas 
where  insufficient  heat  release  occurs  to  affect  the  gas  velocity  sig¬ 
nificantly. 

The  effectiveness  and  practicality  of  the  model  were  demonstrated 
by  excellent  agreement  between  model  predictions  and  experimental 
catalytic  combustion  data.  Data  from  a  hybrid  combustor  are  pre¬ 
sented  and  analyzed  using  a  simplified  model.  Refinement  and  ex¬ 
tension  of  our  current  model  are  in  progress  to  provide  for  the  inclu¬ 
sion  of  multiple  fuel  species,  two-step  reaction  schemes,  and  internal 
heat  removal  from  the  catalytic  section. 

Physical  Assumptions  of  the  Model 

The  following  assumptions  were  made  in  the  derivation  of  gov¬ 
erning  mass,  heat,  and  momentum  balances  for  the  steady-state 
catalytic  combustor  model: 

Numbers  in  brackets  designate  References  at  end  ,»f  paper 
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1  Properties  for  the  gas-phase-are  uniform  throughout  a  reactor 
cross  section  (plug  flow).  Axial  dispersion  is  neglected. 

2  The  solid  temperature  and  fuel  concentration  associated  with 
the  substrate/catalyst  layer  are  similarly  uniform. 

3  Reaction  of  the  fuel  to  products  can  occur  at  the  catalyst  surface 
and  in  the  gas  phase.  Global  (one-step)  Arrhenius  rate  expressions 
are  employed  for  complete  oxidation  reactions  in  each  phase.  As¬ 
signable  orders  in  fuel  and  oxygen  are  incorporated  in  the  homoge¬ 
neous  expression,  while  the  heterogeneous  reaction  is  provisionally 
taken  to  be  first  order  in  the  fuel  and  zero  order  in  oxygen.  The  latter 
assumption  is  considered  adequate  for  the  fuel-lean  conditions  of 
interest  (2). 

4  Transport  of  the  fuel  species  is  considered  to  occur  by  diffusion 
through  the  boundary  layer  and  to  be  balanced  by  the  catalytic  re¬ 
action  rate  on  the  surface. 

5  Heat  released  by  the  catalytic  oxidation  reaction  is  liberated 
at  the  surface  and  subsequently  convected  into  the  gas  phase.  Heat 
and  mass  transport  are  modeled  using  transfer  coefficients  where  the 
driving  force  is  provided  by  temperature  and  concentration  differ¬ 
ences  between  the  bulk  gas  and  the  catalytic  surface. 

6  The  velocity  is  allowed  to  vary  axially,  subject  to  a  Fanning 
friction  factor.  For  a  given  application,  the  transfer  coefficients  are 
evaluated  using  appropriate  correlations  at  estimated  axially  mean 
conditions. 

7  Axial  heat  conduction  in  the  substrate  is  included  in  the  model, 
but  solid-to-solid  radiant  effects,  gas-phase  conduction  and  radiation 
are  neglected.  AU  the  heat  released  by  combustion  in  the  reactor  is 
assumed  to  be  carried  out  of  the  combustor  by  the  gas  stream. 

Conservation  Equations 

The  resulting  model  equations  are  given  below  in  a  form  applicable 
to  a  monolithic  (constant  channel  cross  section)  support.  Similar 
equations  have  been  written  for  packed-bed  or  screen  geometries  but 
will  not  be  presented. 

Fuel-species  balance  (gas-phase) 

f  S^7^  +  SA"  ”P  (-  ^p)  (C/)"t(Co,)"«  +  ah  (C,  -  C„ )  -  0 

(1) 

Fuel-species  balance  (surface) 


k(C/  -  C/„)  »  A<  exp  (-|£)C„ 


Gas-phase  energy  balance 
dT 

-SpuCp  —  +  ah(T.  -  T) 
dz 


+  S(-AH)A*  exp  (-  (C/)"'(Co,)"°  -  0 


S.\,  ^  -  ah(T,  -  T)  +  <r(-A H)A(  exp  ( -  Ct,  =  0  (4) 

dz  ~  '  A  7\ ' 


Momentum  equation 


du  dp  af 

-  (pu)  - - - - -  (pu)  o  =  0 

dz  dz  25 


Continuity  equation 


Equation  of  state 


d 

—  Ow)  ■  0 

dz 


M  Mu 


Stoichiometry 

vCp2 

(uCoj)o 

Boundary  conditions 


'  (vCi)J 


atz  »  0:  Cf  »  (C/)o,  Co,  *  (Co,)o.  u  *  o0,  T  =  To.  — —  ”  0 

dz 

at*  «  0 

dz 

Considerable  simplification  of  the  above  equations  is  possible. 
Equation  (6)  allows  po  to  be  replaced  by  pooo,  while  equation  (7)  and 
(8)  can  be  used  to  bliminate  pressure  and  oxygen  concentration  as 
explicit  variables. 


Dimensionless  Equations  and  Groupings. 

Techniques 

Rearrangement  and  introduction  of  the  dimensionless  .. 

_  z  .  T  -  T0  .  T, -T0  v 

~’um70' 

i  c'  i  .  c>"  .  uC'  .  i 

Off  sm  -  .  yh  f  3  ,  u )  m  _  m  try/ 

( Cf)o  Un(C/)o 

provides  the  following  set  of  catalytic  combustor  equations: 

+  do('h  -  tfr)  m  o  (9) 
,  , .  ___  /  \  , 


JoWf  —  “  D»r  exp  h- 


Solid-phase  energy  balance 


dt 

—  -</«(*,-«> 

df 


— _  —Nomenclature 

Note:  Dimensionless  groupings  found  in 
Table  1  are  not  included. 

A  *  pre-exponential  factor.  A*,  m/s;  Ah, 
(mol/m1)1  "'(s,  where  n  is  total  homoge¬ 
neous  reaction  c-der 

C  ■  concentration  (associated  with  catalyst 
layer  if  subscripted  by  *;  otherwise,  gas- 
phase),  mol/ nr 

Cf  "  specific  heat  capacity  of  gaa  at  constant 
pressure,  J/ikg-K) 

E,  •  activation  energy,  J/mol 
f  “  fanning  friction  factor 
h  “  heat  transfer  coefficient,  J/(m3-s-K) 
(-AH)  *  heat  of  reaction,  J/mol 


•  mass  transfer  coefficient,  m/s 
=  molecular  weight,  kg/mol 

•  reaction  order 

=  Prandtl  number 

•  pressure,  Pa 

•  gas  constant,  J/(mol-K) 

•  void  cross  section,  m3 
”  solid  cross  section,  m3 

•  Schmidt  number 

•  temperature,  K 
1  velocity,  m/s 

■  axial  dimension,  m 

■  thermal  conductivity,  J/(m-s-K) 
1  density,  kg/m3  , 

■  wetted  perimeter,  m 


<4  =  equivalence  ratio 
=  dimensionless  concentration  (see  text  for 
definitions) 

w  «  dimensionless  mass  fraction 
9  =  dimensionless  temperature 
v  =  dimensionless  velocity 
t  =  dimensionless  axial  dimension 

Subscripts 
0  •  inlet 
c  **  catalytic 
f  -  fuel 

H  =  homogeneous 
0„>  or  O  ■  oxygen 
s  •  surface  (solid) 


»«  h,.  a  cue 


• -  J*V-  -  *>  +  B°~  “P  (y~)  +t-  “  0  (12) 


Jeu-~  +  Jf»3 
av  af _ 

<«■"  •/*(!  +  »)-*>* 


subject  to 


— -0 

af 

dti 

f-  1:  — -0 
as 

Equation  ( 10)  and  the  relation  w  ■  uty  can  now  be  used  to  eliminate 
and  Moreover,  upon  substituting  for  the  last  two  terms  in 
equation  (12)  from  equations  (9)— (11)  and  integrating  one  obtains: 

“  •  Pe,  (d  —  B(1  —  «))  (12') 

“f 


subject  to  d  (1)  —  B  (l  -  u(l))  «  0. 

The  definitions  of  the  dimensionless  groupings  occurring  in  the 
above  equations  for  a  monolith  are  given  in  Table  1,  together  with  an 
indication  of  their  physical  significance. 

High-order  approximations  to  the  solutions  of  the  governing 
equations  were  obtained  using  an  adaptive,  variable-order,  finite- 
difference  scheme  developed  by  Lentini  and  Pereyra  (3|.  Much  of  the 
power  of  this  method  results  from  the  combination  of  Newton's 
method  with  iterated  deferred  corrections.  As  for  most  finite  differ¬ 
ence  methods  employed  to  solve  nonlinear  problems,  some  care  must 
he  taken  to  insure  that  the  initial  iterate  lies  in  the  neighborhood  of 
convergence  about  a  solution,  especially  in  regions  having  multiple 
steady  states.  Where  stiffness  and  partial  instability  were  encountered 


it  was  on  occasion  found  necessary  to  employ  parametric  continuation. 
Further  discussion  of  the  use  of  this  method  can  be  found  in  (1]. 

Parametric  Variations 

In  order  to  investigate  the  behavior  of  solutions  of  the  model 
equations,  several  parametric  variations  were  carried  out.  Our  aim 
was  to  study  the  relation  of  the  homogeneous  to  the  heterogeneous 
reaction.  Accordingly,  the  homogeneous  Damkonler  number,  D„a, 
was  varied  for  three  values  of  the  adiabatic  temperature  rise,  B,  with 
the  other  parameters  fixed  at  physically  consistent  values.  This 
variation  can  be  thought  of  as  increasing  the  characteristic  homoge¬ 
neous  reaction  rate  relative  to  a  fixed  catalytic  reaction  so  as  to  as¬ 
certain  the  effect  on  solution  profiles.  Fig.  1  exhibits  conversion  (l  - 
w)  profiles  for  a  case  which  Is  typical  of  a  combustor  operated  fuel-lean 
with  relatively  low  inlet  temperature  (500  K).  B  »  1  was  used,  indi¬ 
cating  moderately  high  heat  release.  The  curve  with  D„»  »  Ogives  the 
conversion  profile  for  the  case  where  only  catalytic  combustion  is 
taken  into  consideration.  The  effect  of  increasing  D,,/,  is  to  increase 
conversion  due  to  contribution  from  the  terms  accounting  for  ho¬ 
mogeneous  combustion.  The  onset  of  gas-phase  combustion  occurs 
at  a  fairly  distinct  point  in  the  bed.  At  this  point  the  gas  phase  has 
attained  a  reasonably  well-defined  homogeneous  light-off  (ignition) 
temperature.  For  combustors  utilizing  hydrocarbon  fuels,  the  values 
for  D„*  needed  here  to  show  significant  homogeneous  contribution 
are  unrealistically  high  (see  discussion  of  kinetic  parameters  in  the 
following).  Nevertheless,  the  clearly  delimited  homogeneous  com¬ 
bustion  zone  exhibited  in  this  case  should  be  typical  of  conditions  for 
which  the  dimensionless  homogeneous  activation  energy,  here  taken 
to  be  40,  is  high.  This  results  in  a  strong  temperature  dependence  of 
the  homogeneous  reaction  rate.  In  fact,  in  the  absence  of  any  catalytic 
reaction  (D0I  =•  0),  no  conversion  is  indicated  for  any  of  the  values  of 
D„(,  used.  Thus,  the  homogeneous  reaction  is  "bootstrapped"  by  the 
heat  released  due  to  catalytic  combustion.  This  strong  dependence 
also  gives  rise  to  numerical  “stubbornness"  as  can  be  appreciated  by 
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Meet  of  variation  of  homogonoous  DafnNoM  or  number  on  conversion 
pr«flto» forth*  cam  with  B  -  O.S,  7,  -  20,  y.  -  20,  D„  -  0.15,  J„  -  J0 
m  >•  +9  m  200,  J?  ■  0,  Pi,  ■  § 
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Flg.1  DMct  of  variation  of  homogeneous  Damkohler  numb*r  on  conversion 
profiles  tor  the  case  wfth  B  «  1. 7»  »  40, 7.  *  20, 0M  «  O.tS,  JM  =  Ja  = 
1,  Jg  ■  200,  J,  -  0.  Pa,  -  5 


inspecting  Pig.  1.  The  integration  scheme  employed  together  with 
continuation  on  DcA  successfully  overcame  this  difficulty. 

Conversion  profiles  obtained  upon  variation  of  D,*  with  B  «  O.S 
and  7a  ■  20  are  presented  in  Fig.  2.  These  values  suggest  fuel-lean 
operation  but  with  a  much  higher  gas  inlet  temperature  (1000  K).  The 
solid  curves  represent  a  family  of  solutions  showing  the  effect  of 
varying  D**.  With  the  lower  value  of  20  for  7 a,  the  temperature  de¬ 
pendence  of  the  homogeneous  rate  is  now  much  weaker.  Conse¬ 
quently,  the  onset  of  homogeneous  contribution  is  no  longer  sharp 
but  gradual.  Such  a  gradual  onset  could  be  expected  from  the  presence 
of  a  weak  but  non-negligible  homogeneous  reaction  occurring 
throughout  the  reactor  length  even  in  the  absence  of  any  catalytic 
reaction.  This  situation  is  indicated  by  the  conversion  profile  with  Dah 
«  0.10  and  Doc  ”  0  in  Fig.  2.  Again  in  this  case,  the  increase  in  gas 
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Rg.3  Unconverted  loaf  traction  (w),  gas  temperature  (9/*),  a 
lamparatur*  (9,/B)  profHas  showing  gat  overtomperature  lor  1 
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case  with  B 
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temperature  due  to  catalytic  combustion  greatly  enhances  the  ho¬ 
mogeneous  reaction  rate  with  the  combined  effect  being  significantly 
increased  conversion. 

Two  more  observations  about  this  series  of  solutions  should  be 
made.  First,  two  dashed  curves  are  shown  in  Fig.  2  which  represent 
additional  steady  states  that  show  better  conversion  than  those  rep¬ 
resented  by  the  solid  curves  (for  the  same  values  of  D„* ).  There  is 
strong  numerical  evidence  that  the  pure  catalytic  model  (i.e.,  with 
3  0)  has  a  unique  solution  for  these  parameter  values.  Hence,  for  this 
case  the  presence  of  multiple  steady  states  is  apparently  attributable 
to  inclusion  of  the  terms  related  to  homogeneous  combustion.  For 
other  combustor  conditions,  multiple  states  have  been  observed  for 
the  pure  catalytic  model  as  well.  Second,  the  phenomenon  of  gas  ov¬ 
ertemperature  can  be  observed  in  Fig.  3  showing  the  profiles  of  un¬ 
converted  fuel  fraction,  gas  temperature,  and  substrate  temperatures 
(both  temperatures  divided  by  B)  for  D.*  =  0.10.  From  the  boundary 
condition  9(1)  -  B(  1  -  u>(l)),  it  is  clear  that  the  dimensionless  gas 
temperature  at  the  reactor  outlet  can  not  exceed  the  adiabatic  tem¬ 
perature  rise.  Yet  for  this  case,  8  does  exceed  B(9/B  >  1)  within  the 
reactor.  Before  reaching  the  exit,  however,  heat  has  been  lost  by  the 
gas  to  the  substrate  thus  cooling  the  gas  below  the  adiabatic  flame 
temperature.  In  this  connection,  it  should  be  noted  that  the  value  used 
for  the  solids  Peclet  number,  Pe,  -  5,  corresponds  to  a  very  conductive 
substrate  (e.g.,  metallic).  For  higher  values  of  Pe.  (£100)  no  such 
overtemperature  is  observed.  The  explanation  of  this  occurrence  is 
thus  related  to  axial  heat  conduction  in  the  substrate.  Indeed  if  this 
conduction  were  not  allowed,  the  condition  9  -  B(1  -  w)  would  be 
imposed  throughout  the  bed,  i.e.,  local  adiabaticity,  with  the  conse¬ 
quence  that  9  could  never  exceed  B.  This  overtemperature  phenom¬ 
enon  is  different  from  substrate  overtemperature  and  transient  ov¬ 
ertemperature  which  have  been  previously  discussed  [4|. 

The  final  variation  of  V,h  to  be  described  was  carried  out  with  B 
*  0.03  and  7*  -  20  and  is  presented  in  Fig.  4.  These  values  are  typical 
of  a  combustor  being  operated  extremely  fuel-lean  at  high  gas  inlet 
temperature  (1000  K).  Thes*  conditions  can  arise  when  a  fuel-air 
mixture  is  partially  oxidized  in  a  preburner  before  entering  the  cat¬ 
alytic  reactor.  Experimental  data  obtained  by  Siminski  and  Shaw  (5| 
f.  jm  such  a  hybrid  combustor  is  discussed  below.  For  the  chosen 
values  of  D,,/,,  a  moderate  homogeneous  reaction  occurs  in  the  absence 
of  any  catalytic  reaction.  When  both  reactions  are  taken  into  account. 

Transactions  of  the  ASME 
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Fig.  4  Effect  at  variation  at  hontoganaaus  DemkaNar  number  on  converdon 
pratHea:  cate  with  B  •>  0.03, 7.  ■  5,  y„  »  20,  JM  *  Ja  «  1,  J,  -  200,  J, 
-  0.  Pe,  -  200 


the  homogeneous  reaction  contributes  significantly  throughout  the 
length  of  the  bed.  However,  in  contrast  to  the  “bootstrap"  phenom¬ 
enon  observed  earlier,  the  reactions  in  this  case  do  not  appear  to  en¬ 
hance  each  other.  Indeed  the  unconverted  fuel  fraction  is  observed 
to  equal  closely  the  product  of  the  fractions  that  would  be  obtained 
for  independent  catalytic  and  homogeneous  reactions. 

Correlations  and  Kinetics 

The  analysis  of  data  which  will  be  subsequently  discussed  was 
carried  out  with  an  option  in  our  computer  program  which  allows  for 
the  dimensionless  groupings  to  be  calculated  from  physical  input 
variables.  The  correlations  used  to  calculate  the  heat  transport 
coefficient,  and  thus  Jh.  were  similar  to  those  discussed  by  Hawthorn 
[6|.  The  mass  transport  number  was  then  computed  assuming  the 
constancy  of  the  ratio  of «/«( Pr)273  to  Jo  (Sc)27*  as  in  Satterfield  [7|. 
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Fig.  S  Experimental  data  (•)  and  modet-generaled  curves  (Hawing  the 
affect  of  Met  temperature  on  conversion 


Vortruba,  et  al.  [8|  have  pointed  out  difficulties  with  these  theoreti¬ 
cally  developed  expressions  and  have  presented  different  expressions 
for  laminar  flow.  Physical  properties  occurring  in  the  correlations  have 
been  taken  to  be  those  of  air  at  inlet  pressure-  and  prov  isionally  at  a 
temperature  equal  to  the  mean  of  the  inlet  tempers:  -ire  and  the  adi¬ 
abatic  flame  temperature.  The  reference  binary  dufusivity  of  the 
specific  fuel-air  system  of  interest  is  provided  as  an  input  variable. 

The  choice  of  suitable  global  homogeneous  reaction  rate  expressions 
is  difficult  A  major  reason  is  that  we  are  currently  treating  a  complex 
reaction  over  a  wide  range  of  physical  conditions  as  a  one-step  reac¬ 
tion.  For  initial  calculations,  we  have  employed  a  second  order  ex- 
presaion  (first  order  in  fuel  and  in  oxygen)  with  constant  pre-expo¬ 
nential  factor.  Values  for  Ah  and  £„*  have  been  chosen  within  ranges 
obtained  by  suitably  interpreting  the  empirical  expressions  given  by 
Longwell  and  Weisa  [9J. 

The  choice  of  an  appropriate  catalytic  rate  expression  is  further 
complicated  by  the  wide  variations  in  such  expressions  appearing  in 
the  literature  and  frequently  by  insufficient  information  being  given 
to  permit  their  use.  References  (2, 6,  I0J  do  give  rate  data  for  hetero¬ 
geneous  combustion  of  hydrocarbon  fuels.  Consonant  with  the  above 
sources,  the  catalytic  rate  expression  used  for  fuel-lean  applications 
was  first  order  in  fuel  and  zero  order  in  oxygen.  Care  was  taken  that 
values  for  Ac  and  £,c  were  within  ranges  suggested  by  the  literature 
and  that  the  overall  chemical  rate  thus  obtained  for  relevant  tem¬ 
peratures  reflected  experimental  values. 

Comparison  With  Catalytic  Combustor  Data 

A  comparison  of  the  model  with  experimental  catalytic  reactor  data 
reported  by  Pfefferle,  et  al.  [11]  was  made,  demonstrating  excellent 
agreement  Computer  solutions  were  obtained  using  physical  input 
data  corresponding  to  the  experimental  operating  conditions.  Thus, 
actual  values  for  reference  velocity,  pressure,  fuel  to  air  ratio,  bed 
length,  channel  diameter,  and  porosity  were  input  parameters.  With 
reference  to  substrate  conductivity,  we  have  successfully  solved  the 
equations  with  values  of  Pe,  up  to  I0\  but  experience  has  shown  that 
such  solutions  are  virtually  identical  to  those  obtained  for  Pe.  »  200, 
which  value  was  adopted  for  this  study.  Fuel  parameters  chosen  were 
those  for  propene.  The  choice  of  kinetic  parameters  is  discussed  in 
the  following.  Finally,  we  observe  a  probable  variance  between  con¬ 
version  as  predicted  by  the  model  and  the  hydrocarbon  conversion 
reported  in  [11].  The  former,  due  to  the  global  assumption  of  the 
model,  means  complete  oxidation,  while  the  latter  appears  to  be  the 
ratio  of  exit  hydrocarbon  (as  Cn)  concentration  to  inlet  fuel  concen¬ 
tration. 

Fig.  5  shows  the  effect  of  inlet  temperature  upon  conversion  with 
combustor  operating  conditions  given  in  the  legend.  Further  experi¬ 
mental  details  can  be  found  in  [11).  The  experimental  data  points  are 
indicated  by  circles,  while  the  solid  curve  was  obtained  from  our 
computer  solutions  for  the  full  model  including  heterogeneous  and 
homogeneous  contributions.  We  observe  the  striking  agreement  in 
behavior  of  the  theoretical  curve  with  the  data. 

For  a  given  inlet  temperature  below  390  K  (point  B),  there  is  a 
unique  solution  which  predicts  that  substrate  and  gas  '  emperatures 
throughout  the  bed  equal  the  inlet  temperature,  the  reaction  is 
kinetically  limited,  and  there  is  no  appreciable  conversion.  In  the 
range  390-520  K,  there  are  two  solutions  for  each  inlet  temperature, 
labeled  as  branches  I  and  E  between  points  Bi  and  B  >  in  Fig.  5.  It  is 
reasonable  to  assume  that  both  branches  represent  stable  solutions. 
The  ignition  branch,  /,  would  be  followed  upon  increasing  inlet 
temperature.  At  a  "light-off’  temperature  of  approximately  490  K, 
catalytic  ignition  has  occurred,  the  substrate  temperature  consider¬ 
ably  exceeds  the  gas  inlet  temperature  and  there  is  substantial  con¬ 
version.  The  values  used  for  the  two  catalytic  rate  parameters.  A,  and 
Eac,  directly  affect  the  temperature  at  which  this  activation  of  the 
catalyst  occurs.  Values  of  914  m/s  for  A,  and  50.2  kj/mol  for  were 
used  for  this  study.  These  conform  quite  well  to  those  cited  in  the 
literature,  especially  by  Hawthorne  [6|. 

The  extinction  branch,  labeled  E  in  Fig.  5,  would  he  followed  upon 
decreasing  gas  inlet  temperature  after  having  obtained  catalytic 
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PRESSURE:  101  kP» 

INLET  TEMPERATURE:  662  K 
FUEL:  PROPENE 

REFERENCE  VELOCITY:  18.3  «i/i 


Fig.  •  ExpartmanUJ  data  (•)  and  mad«l  gwwrateg  aava  Mwwtng  the  affect 
of  fuel/air  ratio  on  convarMon 


ignition  and  having  paaaed  the  point  Bj  on  the  curve.  At  approxi¬ 
mately  390  K,  a  sharp  extinction  is  predicted.  The  model  thus  indi¬ 
cates  the  presence  of  a  hysteresis  effect.  This  phenomenon  was  ob¬ 
served  by  Pfefferle.  et  at  (11.  p.  6]  and  has  been  observed  and  dis- 
cussed  by  Hlavacek  and  Vortruba  [12], 

For  temperetuna  in  the  range  of  390-660  K  the  conversion  pre¬ 
dicted  <  - 1  n  be  attributed  almost  totally  to  the  catalytic  reaction  with 
no  significant  contribution  from  homogeneous  reactions.  A  good 
measure  of  rate  domination  for  catalytic  oxidation  is  provided  by  the 
surface  fuel  concentration  to  gas  fuel  concentration  ratio  (CfJCf). 
When  this  ratio  is  dose  to  zero,  mass  transport  limitation  exists;  when 
it  is  dose  to  one,  the  rata  is  kinetically  limited.  For  420  K  this  ratio 
averages  about  0.14  throughout  the  bed  indicating  that  the  surface 
reaction  rata  is  still  a  factor  in  determining  the  overall  rate.  Hence, 
increased  conversion  is  observed  at  higher  temperatures  in  this  range 
as  a  result  of  enhanced  surface  reaction  rate.  Indeed  by  644  K  the 
model  predicta  that  the  ratio C/./C/  equals  approximately  0.009,  in¬ 
dicating  that  the  reaction  is  approaching  mass  transport  limitation. 
The  computer  generated  curve  clearly  predicts  the  behavior  of  the 
data  in  this  interval. 

A  second  pronounced  increase  in  conversion  is  predicted  in  the 
interval  750-850  K.  At  the  lower  temperature  the  conversion  is  72 
percent,  while  it  exceeds  99  percent  for  850  K.  The  curve  in  Fig.  5  la¬ 
belled  “catalytic  only”  is  obtained  by  setting  the  homogeneous  reac¬ 
tion  rate  to  zero  and  thus  considering  terms  arising  only  from  the 
catalytic  reaction.  Comparison  of  this  curve  to  the  solid  curve  reveals 
that  the  rise  in  conversion  in  this  range  (as  well  as  some  contribution 
from  650  K  on)  results  from  terms  accounting  for  homogeneous 
combustion.  Thus,  the  homogeneous  reaction  is  ignited  or,  perhaps 
better,  supported  by  the  catalytic  reaction,  at  observed  by  Pfefferle, 
etal.  [Ilf.  This  “bootstrapping”  phenomenon,  observed  in  the  para¬ 
metric  variations  discussed  previously,  is  further  illuminated  upon 
examining  the  curve  labelled  “homogeneous  only”  obtained  when 
terms  for  the  catalytic  reaction  are  neglected.  The  temperature  needed 
to  achieve  homogeneous  "light-ofT  in  the  absence  of  any  catalytic 
reaction  is  1050  K  (Fig.  5).  Indeed,  axial  conversion  profiles  obtained 
from  the  model  with  catalytic  and  homogeneous  terms  begin  to  di¬ 
verge  from  those  obtained  with  only  catalytic  terras  at  a  reactor 
location  in  which  the  gas  temperature  is  nearing  this  homogeneous 
ignition  temperature. 

The  gas  temperature  range  for  which  the  model  predicts  significant 
contribution  from  the  homogeneous  reaction  is  dependent  upon  the 


values  used  for  At,  and  £„>,  Consistent  with  the  literature  cited  pre¬ 
viously,  At,  was  taken  to  be  2.87  x  109m'/(mol-s)  and  £a*  to  be  167 
kJ/mol.  Experimental  data  suggest  the  onset  of  significant  homoge¬ 
neous  reaction  at  somewhat  higher  temperatures  than  does  the  model 
with  these  values.  No  attempt  was  made  to  vary  A *  and  £„*  to  obtain 
a  better  fit  Nevertheless,  the  agreement  of  model  and  data 
throughout  the  temperature  range  300-900  K  is  striking. 

Reactor  conversion  data  and  model-generated  curve  are  given  in 
Fig.  6  as  a  function  of  fuel-to-air  ratio  {FI A  wt)  with  other  operating 
conditions  fixed  as  indicated.  For  FI  A  values  below  0.015,  the  con¬ 
version  obtained  is  indicated  by  the  model  to  be  due  to  mass  transport 
limited  catalytic  oxidation.  The  inlet  temperature  of  662  K  is  suffi¬ 
ciently  high  to  allow  the  catalyst  to  remain  active  even  to  extremely 
lean  conditions.  With  the  kinetic  parameters  chosen,  both  data  and. 
modal  exhibit  a  marked  increase  in  conversion  around  a  FI  A  of  0.025. 
This  increase  is  due  to  the  contribution  from  the  terms  accounting 
for  homogeneous  combustion,  as  can  be  seen  from  comparison  with 
the  “catalytic  only”  curve.  The  rise  suggested  by  the  data  is  somewhat 
more  rapid  than  that  reflected  by  the  curve,  as  was  also  the  case  for 
Fig.  5.  In  that  the  sharpness  of  the  homogeneous  onset  is  strongly 
dependent  on  the  value  obtained  for  7*  and  thus  upon  £,*,  as  was 
noted  in  the  foregoing,  a  value  for  £„>,  higher  than  167  kJ/mol  is 
suggested.  Further  investigation  into  kinetic  values  and  sensitivities 
is  being  undertaken.  However,  excellent  agreement  with  the  data  is 
again  evident. 

Comparison  With  Data  From  a  Hybrid  Combustor 

The  results  of  a  catalytic  combustor  development  program  carried 
out  in  this  laboratory  are  discusaed  elsewhere  by  Siminski  and  Shaw 
(5).  In  order  to  achieve  stability  at  low-power  (idle)  conditions,  a  hy¬ 
brid  combustor  in  which  the  fuel-air  mixture  is  partially  oxidized  in 
a  preburner  prior  to  entering  the  catalytic  reactor  was  designed, 
constructed,  and  tested.  At  this  time,  detailed  calculations  with  the 
full  model  for  this  hybrid  system  have  not  been  completed.  We  shall, 
however,  provide  a  simple  analysis  of  the  data,  based  on  the  as¬ 
sumption  of  catalytic  rate  limitation  either  by  mass  transport  or  by 
surface  kinetics,  neglecting  homogeneous  contribution. 

The  oxidation  of  fuel  in  the  prebumer  results  in  extremely  lean 
operation  for  the  catalytic  reactor  and  high  gas  inlet  temperature.  For 
these  conditions,  the  assumption  that  substrate  and  gas  temperatures 
throughout  the  bed  approximately  equal  the  inlet  temperature  (and 
consequently  that  the  velocity  remains  essentially  constant)  appears 
reasonable.  On  this  basis  the  model  equations  simplify  and  can  be 
integrated  to  obtain: 

In  u  (1)  •  -  D„r  (kinetics  limited)  (13) 

Inu  (1)  “-«//>  (transport  limited)  (14) 

If  similarity  of  heat  and  mass  transport  and  turbulent  flow  in  the 
catalytic  reactor  is  assumed,  along  with  constant  pressure  throughout 
the  prebumer  section,  then,  accounting  for  the  various  dependencies 
of  Doc  snd  Jo  on  To,  we  have 

In  w  (1) - exp  {—E0,/RT„)/T„  (kinetics  limited)  (15) 

or 

In  w  (1)  ~  -  T0° 1  (transport  limited)  (16) 

The  Tu-dependence  of  the  unreacted  fuel  fraction  at  the  exit  &>  ( 1 )  is 
thus  much  stronger  (for  reasonable  values  of  E„ )  Cor  kinetics-limited 
than  for  mass-transport-limited  operation. 

The  subject  data  for  HC  (hydrocarbon)  and  CO  oxidation  are 
plotted  in  such  ways  as  to  lead  to  linear  graphs  if  the  data  were  to 
represent  mass-transport  limitation  according  to  equation  (16)  in  Fig. 
7  or  kinetics  limitation  according  to  equation  (15)  in  Fig.  8  and  9.  The 
experimental  data  for  different  monolith  geometries  were  corrected, 
before  plotting,  according  to  length  and  hydraulic  diameter  depen¬ 
dences  implied  by  the  equations. 

Fig.  7  shqws  that,  for  most  catalysts,  the  Tu-dependence  of  -In  w(l) 
is  much  too  high  to  be  consistent  with  a  mass-transport  limited  op- 
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eration.  For  CO,  the  values  —In  w(l)  for  these  catalysts,  however,  show 
a  decreasing  To-dependence  at*  high  T0,  suggesting  the  onset  of  the 
mass-transport  limited  regime,  and  thus  the  experimental  values  of 
-In  u(l)  should  nearly  equal  Jo-  The  indicated  range  of  Jo  values 
( 1.5- 1 .9)  is  very  close  to  that  calculated  for  a  CO-air  mixture  at  1200 
K  and  303  IcPa.  In  contrast,  the  values  of  -In  u(l)  for  the  HC  con¬ 
version  data  do  not,  for  any  catalyst,  achieve  low  enough  temperature 
sensitivity  to  indicate  mass-transport-limited  operation. 

Examination  of  the  test  plot  for  kinetic  limitation  of  CO  oxidation. 
Fig.  8,  shows  reasonably  consistent  linear  trends  especially  for  lOVTo 
between  1.0  and  1.4.  Interpreting  the  slope  of  the  visually  fitted  line, 
one  obtains  (E„)co  a  50.2  kJ/moL  This  value  is  very  close  to  that 
given  by  Hawthorn  (6)  for  alumina-based  noble  metal  catalysts.  For 

e 


HC 


higher  temperatures  (0.8  <  IOVT0  <  1.0),  the  graphs  appear  to  be 
flattening  out,  i.e..  tending  toward  positive  slope  as  would  be  expected 
for  the  onset  of  mass-transport  domination. 

For  the  HC  data  presented  in  Fig.  9,  the  linear  trend  is  exhibited 
throughout  the  entire  temperature  range.  This  result  is  consistent 
with  the  previous  observation  that  the  HC  conversion  did  not  appear 
to  be  mass-transport  limited.  The  slope  of  the  line  indicated  in  the 
figure  gives  (£oc)hc  a  29.3  kJ/mol.  This  value  is  well  below  the  range 
48-60  kJ/mol  suggested  by  Hawthorn.  However,  given  the  range  of 
values  reported  for  (£„c  )hc  in  the  literature  and  the  intrinsic  uncer¬ 
tainties  in  modeling  a  hybrid  system,  e.g.,  the  presence  of  significant 
amounts  of  CO?  and  HjO  and  the  absence  of  knowledge  of  the  specific 
fuel  species  present  in  the  feed  to  the  catalytic  reactor,  this  variance 
would  not  appear  significant.  The  assumption  of  no  contribution  from 
homogeneous  combustion  was  supported  by  the  experimental  ob¬ 
servation  that,  even  at  the  highest  operating  temperatures,  no  con¬ 
version  was  obtained  in  the  reactor  section  in  the  absence  of  a  catalyst. 
That  the  hybrid  combustor  appears  kinetically  limited  for  HC  oxi¬ 
dation  at  higher  temperatures  than  would  be  expected  from  work  with 
automotive  afterburners  and  “pure”  catalytic  combustors  is  most 
likely  due  to  the  considerably  higher  gas  velocities  involved.  It  is 
somewhat  surprising  that  CO  oxidation  appears  to  be  tending  to 
mass-transport  limitation  at  temperatures  at  which  the  HC  oxidation 
is  not  showing  this  trend.  Further  modeling  effort  is  planned  to  pro¬ 
vide  a  fuller  understanding  of  hybrid  combustor  operation. 

Summary/ Conclusions 

•  A  variable-order  finite-difference  method  was  employed  to  solve 
the  two- point  boundary  value  problem  which  constitutes  the  model. 
The  approach  successfully  overcame  the  nonlinearity,  stiffness,  and 
partial  instabilities  encountered. 

•  Homogeneous  combustion  must  be  considered  in  any  general 
model  which  attempts  to  describe  catalytic  combustor  operation. 
Under  certain  conditions,  homogeneous  combustion  is  supported  or 
“bootstrapped”  by  the  heat  released  from  catalytic  combustion. 

•  Comparisons  of  catalytic  combustor  data  and  model  predictions 
show  excellent  agreement.  When  the  model  is  used  to  examine  con¬ 
version  as  a  function  of  gas  inlet  temperature  for  typical  lean  operating 
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conditions  with  a  noble-metal  catalyst,  a  sharp  catalytic  light-off 
temperature  is  predicted,  as  is  a  light-off/extinction  hysteresis.  At 
higher  temperatures  the  onset  of  gas-phase  combustion  occurs,  re¬ 
sulting  in  further  increased  conversion,  which  rapidly  approaches  100 
percent. 

•  Data  from  a  hybrid  combustor  are  reasonably  consistent  with 
a  simplified  analysis  of  catalytic  combustion,  neglecting  homogeneous 
contribution.  Over  the  temperature  range  investigated,  CO  oxidation 
tends  from  surface  kinetics  limitation  with  (£ac)co  »  50*2  kJ/mol 
toward  mass  transport  limitation.  Hydrocarbon  conversion  data, 
however,  exhibit  kinetically  limited  behavior  throughout  the  tem¬ 
perature  range  with  (£„ c)hc  a  29.3  kJ/mol. 
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Catalytic  Flame  Holding  Model 
Douglas  Cornell 
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A  first  model  of  the  aerothermochemistry  associated  with  catalytic 
flameholder  devices  is  developed.  This  work  was  performed  as  part  of  an  NSF 
summer  research  participation  project  and  is  being  included  here  because  it 
was  not  published  in  the  open  literature. 


CATAL YT IC  FLAME  HOLDING  MODEL 


The  catalytic  flame  holder  whose  properties  are  modeled  in 
this  report  consists  of  a  honeycomb  support,  the  surface  of  which  is 
coated  with  a  catalyst  which  can  promote  the  oxidation  of  hydrocarbon 
fuels.  Its  properties  are  investigated  with  respect  to  use  in  a  jet 
engine  afterburner.  The  development  of  the  model  is  influenced  by  the 
finding  of  Zukowski  and  Marble1,  that  with  a  conventional  flame  holding 
device  a  critical  time  for  ignition  can  be  determined  (L/V^)  where  L 
is  the  length  of  a  recirculation  zone  behind  the  bluff  body  and  V.  is 
the  blow  off  velocity.  A  parcel  of  combustible  gas  mixture  must  be 
ignited  in  less  time  than  is  required  to  pass  by  the  recirculation 
region.  This  ignition  time  is  independent  of  flame  holder  size  or 
geometry.  The  experimental  facts  concerning  the  recirculation  region 
are  as  follows: 

1)  The  length  of  the  wake  region  is  proportional  to  the  diameter 
of  a  cone  shaped  flame  holder  and  is  proportional  to  the  square 
root  of  diameter  of  a  cylindrical  flame  holder.1 

2)  The  temperature  behind  the  flame  holder  is  uniform  and  is  about 
90%  of  the  adiabatic  flame  temperature. 

3)  The  region  is  not  well  stirred.2 

4)  The  residence  time  of  the  gas  in  the  recirculation  zone  is  much 
longer  than  the  time  required  for  gas  flow  past  the  recirculation 
zone. 3 


1  Edward  £.  Zukowski  and  Frank  E.  Marble,  Proceedings  Aerothermochemistry 
Gas  Dynamics  Symposium  (1955),  p.  265. 

2  A.  A.  Westenberg,  W.  G.  Berl ,  and  J.  L.  Rice,  ibid,  p.  211. 

3  T.  A.  Bovina,  Seventh  Symposium  (International)  on  Combustion, 
Butterworths ,  London,  p.  692. 
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This  report  proceeds  to: 


a)  determine  the  flow  split  of  the  fuel,  air  mixture  through  and 
around  the  catalytic  monolith. 

b)  using  the  results  of  calculations  by  A.  E.  Cerkanowicz4 5  for 
reaction  in  the  catalytic  monolith,  determine  from  the  predicted 
flow  into  the  recirculation  region  whether  or  not  stable  ignition 
and/or  flame  holding  will  occur. 

c)  program  these  calculations  for  a  digital  computer  so  that  the 
operating  characteristics  of  the  flame  holder  can  be  readily 
determined  as  a  function  of  experimentally  adjustable  parameters. 

Flow  Split 

The  flow  split  Is  determined  on  the  basis  that  pressure  drop 
for  the  flow  around  the  catalytic  monolith  is  the  same  as  for  the  flow 
through  the  catalytic  monolith.  In  other  words,  no  radial  pressure 
gradients  exist  in  the  afterburner.  The  flow  around  the  flame  holder 
is  calculated  from  Bernoulli's  equation  for  the  change  In  pressure 
resulting  from  a  change  In  area  of  the  flow  path. 

(i)  PrP2  *  (p/2)(q/a2)2ci-(a2/a1)21 

where  Q  is  the  volume  flow  rate  and  A^  and  Ag  are  the  areas  of  the  duct. 
The  pressure  loss  due  to  friction  has  been  evaluated  and  found  to  be 
only  about  1%  of  that  due  to  change  in  area.  The  general  equation  for 
pressure  drop  in  a  duct  where  chemical  reaction  and  friction  losses  are 
occurring  i$s 


4  Exxon  Report  No.  7514  955. 

5  Ascher  H.  Shapiro,  "The  Dynamics  and  Thermodynamics  of  Compressible 
Fluid  Flow"  VI,  Ronald  Press  Co.,  N.Y.  (1953). 


(2)  dP/P  -  [kM^ /(l-M^)][-(l/CpT)dhr-(l+(k-l)M^)(2f/0)dx  ♦  (1/W)dw] 

where  k  ■  Cp/Cv,  M  is  the  mach  number,  f  the  friction  coefficient,  D 
the  hydrodynamic  diameter  *  4A/(dAw/dx)  and  W  the  molecular  weight.  The 
existing  Exxon  catalytic  reactor  model  uses  a  simpler  form  in  which 
(k-l)M2  and  dW  are  considered  negligibly  small  and  are  not  included. 

Equations  have  been  developed  by  Cerkanowicz4  which  allow  the 
evaluation  of  the  pressure  drop  through  the  catalyst  tubes  under  conditions 
where  the  reaction  rate  at  the  catalyst  is  much  larger  than  the  transport 
rate  to  the  surface  and  where  the  Pecklet  number  is  large  (>100).  The 
fuel  conversion  (W)  is  first  evaluated  for  no  pressure  loss  and  is 
subsequently  corrected  for  pressure  drop. 

(3)  const.  P  (1+B)lnw  +  B(l-w)  * 

(4)  P  -  (1+A1)P0 

(5)  A,  -  -(Jf/Jc)[HB(l-w)] 

(6)  (1+B)lnw  +  B(l-w)  »  -JD(l+y2) 

The  definitions  of  the  symbols  and  their  values  are  given  in  the  Appendix. 
The  computer  program  calculates  W  from  (3)  using  the  Wegstein  method, 
calculated  A^  from  (5)  and  again  calculates  W  from  (6).  Equations  (5) 
and  (6)  are  solved  cyclically  until  no  further  change  in  W  and  A^. 

The  computer  program  is  provided  with  an  Initial  estimate  of 
the  fraction  of  flow  around  the  catalytic  flame  holder.  It  then  compares 
the  calculated  pressure  drops  for  flow  through  and  around  the  catalytic 
monolith  and  Increments  the  fraction  of  flow  around  the  flame  holder 
until  equal  pressure  drops  are  attained  by  equations  (1)  and  (4). 


Flame  Holding 

Because  the  velocity  of  the  gas  flowing  around  the  catalytic 
flame  holder  Is  greater  than  that  of  the  exit  gas  from  the  flame  holder. 
It  is  supposed  that  entrainment  will  cause  a  recirculation  zone  to  exist 
as  in  the  case  of  the  bluff  body.  This  recirculation  zone  may,  however, 
be  longer  than  in  the  case  of  a  bluff  body  (see  below).  It  is  assumed 
that  the  flame  develops  at  a  distance  XI  from  the  catalyst  as  shown  in 
Figure  1 . 


FIGURE  1 

Xi  is  supposed  equal  to  the  product  of  a  mean  convective  velocity  and 
reaction  time  as  In  the  analysis  of  Kundu6,  et  al. 

(7)  Xi  »  (U0/2)(p0/Mt) 

The  reaction  rate  which  appears  in  (7)  Is  evaluated  at  the  original 
concentration  and  at  recirculation  zone  temperature.  Wt  is  estimated 
from  the  relationship  between  reaction  rate  and  laminar  flame  propagation 
velocity. 

(8)  SL  -  (koWL/Cpo)°-S/0o 
We  therefore  obtain 

(9)  Xi  »  (U0k0/2p0Cp0SL2)[exp(E/RTR)/exp(E/RTf)] 

6  Kundu,  K.  M. ,  Baneriee,  0.,  and  Bhadurl,  0.,  Combustion  Science  and 
Technology  17.,  153  (1977). 
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The  ratio  of  exponentials  occurs  because  the  ignition  source  for  laminar 
flame  propagation  Is  at  flame  temperature  whereas  the  ignition  source 
for  the  flame  holder  is  at  recirculation  temperature.  The  temperature 
of  the  recirculation  zone  must  be  obtained  from  a  heat  balance  equation 
for  the  recirculation  zone. 

The  flow  rate  Into  the  recirculation  zone  from  the  catalyst  Is 
M  and  from  the  turbulent  transport  through  the  recirculation  zone  boundary 
Is  per  unit  length.  It  Is  assumed  that  the  mass  of  gas  entering  the 
recirculation  zone  at  temperature  T  from  the  catalyst  leaves  the 
recirculation  zone  at  the  temperature  Tr. 


Heat  entering  *  M  Cp  T  +  (CpTf  -  CpTr)(L  -  XI) 
Heat  leaving  *  M  Cp  Tr  +  M-j  (CpTr  -  CpTo)(X1) 


r-";  v  t  Mi<cPvcpVa'V'  •  vr  ♦  ivcpvyo:txi) 

« 

Tne  value  of  tt|  can  be  evaluated  from  the  measured  recirculation  zone 
residence  times  for  bluff  bodies.3 

(11)  tr  «  KD/U  where  K  %  104 

0 

2 

(12)  IM  -  o  (volume  of  recirculation  zone)  a,  o  0(0^2)  irLlUQ 


On  a  per  unit  length  basis  we  obtain 
(13)  M1  *  Oir02Uo/132.5 

The  recirculation  zone  length  (L)  Is  estimated  from  the  work  of  Landis 
and  Shapiro7  and  Forstall  and  Shapiro8,  who  investigated  the  turbulent 
mixing  of  concentric  jets  of  gas. 


7  F.  Landis  and  A.  H.  Shapiro,  Proc.  Heat  Transfer  and  Fluid  Mech.  Inst., 
Stanford  U.  Press  (1951). 

3  W.  Forstall  and  A.  H.  Shapiro,  J.  of  Applied  Mechanics  (1950),  p.  381. 


FIGURE  2 


L  Is  the  length  over  which  the  Inner  jet  velocity  (U^)  remains  unchanged 
when  flowing  Into  a  concentric  jet  at  velocity  (U2).  They  found 
(14)  L  *  Add+aUj/U^ 

This  Is  taken  to  be  the  recirculation  zone  length  In  this  model.  For 
the  velocities  typical  of  the  flame  holder,  L  may  be  approximately 
double  that  for  a  bluff  body. 

The  computer  program  solves  equations  (9)  and  (10)  simultaneously 
using  Tf  as  the  Initial  estimate  for  Tr.  XI  is  first  calculated  from 
(9)  and  then  Tr  Is  calculated  from  (10).  The  new  value  of  Tr  from  (10) 

Is  then  used  to  recalculate  XI  from  (9).  The  cycle  repeats  until 
successive  calculations  of  Tr  are  within  a  prescribed  limit.  It  Is 
Important  to  approach  Tr  from  the  high  side  as  otherwise  a  converging 
solution  might  not  be  obtained. 

Ignition  is  assumed  to  take  place  if  XI  calculated  from  (9) 
using  the  catalyst  stream  exit  temperature  Is  less  than  L.  If  XI (T)  is 
greater  than  L  but  XI (Tr)  Is  less  than  L,  a  stable  flame  Is  possible 
buc  may  require  Ignition.  If  X1(Tr)  Is  greater  than  L,  blow  off  occurs. 

A  listing  of  the  computer  program  Is  attached  together  with 
the  results  of  some  typical  calculations.  The  evaluation  of  parameters 
required  as  Input  Is  discussed  In  the  Appendix. 


Results  of  Calculations 


Typical  jet  engine  afterburner  Inlet  conditions  are  T  ■  950°K 

Uo  ■  130  m/sec 


The  values  given  are  In  mks  units. 


>  .25 

To  »  950* 

Table  1. 

*  uo  « 

Variation 
130  m/sec 

of  TO 

02  -  .27 

01  »  .40 

TO 

AP/P 

G 

W 

XI 

RL 

.002 

.0712 

.995 

• 

.003 

.0692 

.985 

0 

.0792 

1.259 

.004 

.0672 

.975 

.005 

.0633 

.960 

.00011 

.1765 

1.57 

.007 

.0557 

.92 

.058 

.256 

2.07 

.009 

Blow  Out 

Table  2. 

Variation 

Of  TL 

-  .007 

To  -  950* 

Uo 

•  130 

02  ■  .27 

01  -  .40 

TL 

AP/P 

S 

W 

XI 

RL 

.25 

.0557 

.92 

.058 

.256 

2.07 

.20 

.0557 

.92 

.090 

.295 

2.05 

.15 

.0484 

.88 

.23 

.637 

2.47 

.10 


Blow  Out 


To  »  950° 
02 


Table  3.  Variation  of  Blockage  Ratio 
UO  -  130  TO  -  .007  TL  »  .25 


D1  *  .40 


02 

AP/P 

S 

W 

XI 

RL 

.27 

.0557 

.92 

.058  . 

.256 

2.07 

.30 

.0859 

.86 

.113 

.4533 

2.49 

.33 

.1311 

.74 

.203 

1.11 

3.07 

.36 

Blow 

Out 

TO 

Table  4. 
-  .007 

Variation 
Uo  «  130 

of  Inlet  Temperature 
TL  -  .25  02  - 

.27 

T 

AP/P 

G 

W 

XI 

RL 

950 

.0557 

.92 

.058 

.256 

2.07 

925 

.0572 

.92 

.063 

.331 

2.06 

900 

.0588 

.92 

.069 

.4439 

2.06 

875 

Blow 

Out 

To  »  950 


Table  5.  Variation  of  Velocity 
TO  -  .007  TL  •  .250  01  -  .40 


D2  »  .27 
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Calibration  of  Model  with  Experimental  Data 

The  parameter  used  In  this  model  whose  estimated  value  is  in 
greatest  error  is  probably  the  laminar  flame  speed  (S^).  Values  of  SL 
have  apparently  not  been  measured  at  the  high  temperatures  (950°) 
characteristic  of  afterburner  inlets.  A  convenient  calibration  can  be 
accomplished  using  an  experimental  value  for  blow  off  velocity.  $L 
can  be  adjusted  so  as  to  make  XI  >  RL  at  the  experimental  blow  off 
velocity  and  with  Tr  *  T  (of  catalyst  flow). 

Conclusions 

The  blow  off  condition  is  approached  rather  explosively  for 
the  variation  of  any  of  the  experimental  parameters.  This  is  because 
an  increase  in  XI  reduces  the  temperature  of  the  recirculation  region 
which  would  in  turn  require  a  larger  XI  for  ignition. 

The  optimum  flame  holder  geometry  cannot  be  reliably  determined 
until  calibration  of  the  model  with  experimental  data  has  been  carried 
out.  For  preliminary  design,  TO  ■  .007  m,  TL  ■  .25  m  appears  favorable. 

A  comparison  of  flame  holder  drag  with  bluff  body  drag  for  the 
same  blockage  may  be  made  for  a  typical  geometry  which  provides  adequate 
flame  holding.  For  TO  «  .007,  TL  »  .25,  01  *  .40,  02  ■  .27  we  note 
that  8%  of  the  flow  goes  through  the  catalyst.  Using  the  equation  for 
dPI  with  g  »  .92  shows  that  a  22%  decrease  in  pressure  drop  as  compared 
with  bluff  body  (g  *  1)  results. 


Appendix 
Program:  CATRE 

INLET  CONDITIONS  (ANY  UNITS  FOR  WHICH  F  -  MA) 

P  (PRESSURE),  TO  (TEMPERATURE),  VO  (VELOCITY),  TD  (TUBE  DIAMETER), 

TL  (TUBE  LENGTH),  EQ  (EQUIVALENCE  RATIO),  R  (GAS  CONSTANT), 

WM  (MOLECULAR  WEIGHT),  H  (HEAT  OF  COMBUSTION  PER  UNIT  MASS), 

CPM  (MEAN  HEAT  CAPACITY  PER  UNIT  MASS),  CP  (HEAT  CAPACITY  FOR  FLAME), 
CPO  (HEAT  CAPACITY  INLET),  D1  (DUCT  DIAMETER),  D2  (MONOLITH  DIAMETER) 
Z  (FRACTIONAL  OPEN  AREA  OF  MONOLITH),  TF  (FLAME  TEMPERATURE), 

G  (FRACTIONAL  FLOW  SPLIT),  DG  (INCREMENT  IN  FLOW  SPLIT), 

Q  (UNITS  CONVERSION  FACTOR  FOR  VELOCITY  »  1.0  FOR  MKS) 

Jf  -  32(TL)/(TD)(Rey) 

Laminar  flow  occurs  In  the  catalyst  tubes. 

JQ  -  (17.03) (TL/TD) ( l+{ . 067 ) (TD/TL ) ( Rey ) ) *  *S/Rey 

B  •  (FO)H/(l+FO)(CPM)(Tmean) 

CPM  Is  mean  heat  capacity. 


The  laminar  flame  speed  for  hydrocarbons  with  more  than  3 
carbons  Is  independent  of  the  number  of  carbons.  It  is  not  especially 
sensitive  to  pressure.  for  an  equivalence  ratio  of  1.0  was  calculated 
from  an  empirical  equation  of  Ougger  and  Simon9;  S  =  10  +  3.42  x  10-4  To2 
cm/sec  for  temperatures  up  to  615°K.  The  dependence  on  $  was  determined 
by  supposing  that  U  *  Uo  +  a  ($-1.0)2. 

The  value  of  a  was  determined  by  fitting  data  of  Spalding10 
for  propane. 

Ko  -  Thermal  Conductivity 

Ko  *  (2. 5) (viscosity) (Cv)  from  kinetic  theory 


9  Gordon  l.  Ougger  and  Dorothy  M.  Simon,  4th  Symposium  on  Combustion, 
p.  336. 

10  0.  B.  Spalding,  "Some  Fundamentals  of  Combustion,"  Butterworths  (1955), 


0.111  0.113  0.1037310 
0.0480 
0.4779 
0.7421 
0.9945 
1.3384 
2.0080 
5.9088 

TO*  950.  H*  0.431E+08  CPR*  0.124E+04  CPO*  0.114E*04  Z«  0.750 

TF*2974.  CO*  0.844E+03  P*  0.233E+05 

UO*  190.0  T 8*0. 007  TL*  0.250  E0>  1.0000 

01*0.400  02*0.270  6*. 900  86*. 010 

BLOU  OUT 

XI*5.9488  RL*2.2457  TR*  >739.1  6*. 900  FRAC  P  LOSS*. 1111 

08.35.20  >«  catr*  fortran 
EDIT  a 

>p  200 

READ<5,*>  R ,S,UN,H,CPH,CP,CPO,Z,CVf E,Q 

2  READ (5(*,END* 150)  P,TO,UO, TB,TL,Eflf 81 ,82,6,16 
A0*3.14«*D1/2.)«*2 

At *3. 14*(02/2. )**2 

A*Z*A1 

A2»A0-A1 

D*P«UN/<R*TO) 

VIS*S0RT(R*T0*UH/3. 14)/(i.0E23*S) 

TC*2.5*VIS*CW 

FR*8*U0*A0 

J*0 

3  FRA*G*FR 
FRT«<1.-6)*FR 
UT0*FRT/(8*A) 

REY»D*TD*UTO/VIS 

P1«P”(FRA**2/(2.*D*A2**2) )*< 1 .-(A2/ (B*AO) )**2) 

DP1*(P-P1)/P 

CF*2.*TL*1 4./(T8*REY) 

CE*P/(8*UT0**2) 

F0*.0583*E0 

CB«17.03*(TL/T8)*< 1 .♦.047*(TD/TL>*REY>**.45/REY 
8«F0*H/( t 1 .  ♦F0)*CPN*1 500. ) 

U*.  I 
NC*1 
n*i 

5  XU*EXP((8*M-C8-8)/(1.*8) ) 

PRINT  112,fl,XU 
HM*1 

CALL  CONV(UfXMf 1 , NC) 

60  T0(4,5),NC 

4  A4«- <  CF/CC ) • ( 1 .♦8*(1  .-U) ) 


fi*  I  wu 

XM1«£XP< <D*U1-C0-A472.-» J/1 1 .*M  > 

PRINT  112. K 
K*K+1 

CALL  C0NV1U1 ,XU1 ,2, NO 
60  TO  <8,7),NC 

P2*11.-(CF/CE)*11.*B*(1.-U1)))«P 
DP2*(P-P2)/P 
PRINT  1 1 3,DP1 ,DP2,XU1 
IF1A8S(DP1-DP2).LT..01)  60  TO  11 
IF1DP1.LT. BP2)  60  TO  9 
IF1J.EQ.1)  60  TO  10 
G«G-DG 
60  TO  3 
Ja1 

6«G+D6 
GO  TO  3 
G»G-DG/2 

T*T0*(1.+B»11.-U1 )) 

UA*6*A0*U0/A2 
UT*(T/T0)*(1 .-6)*(JO*AO/A 
IF (EQ.6T.1 .0)  60  TO  12 
TF«TO*H*FO/CPN 
60  TO  13 

TF*TO>H*FO*12.-EQ)/CPH 

Y1»1.*3.*UT/UA 

RL«T1M.*D2 

SLO*< . 10>.342E_5»T0»*2)*Q 
SL1»SL0*<1 .0-2.623*IEQ-1 .0>»*2) 

EXCR»D*D2*U0/132.3 

TR«TF 

XI«1UA*TC/(2.*B*CPO*SL1*«2)>*EXP<E/1R*TR)>/EXP(E/(R»TF>) 

PRINT  100, XI 
IF1XI.6T.RL)  GO  TO  U 

TRC»(EXCR*<CP*TO<XI+CP*TF*(RL-XI))*FRT*CP*T)/(FRT*CP*EXCR*CP*RL) 
IF  1 ABS(TRC-TR) .LT.50. )  60  TO  15 
TR*TRC 
GO  TO  14 

XI1»1UA*TC/<2.*D*CP0*SL1**2))«EXP<E/(R*T)>/EXP<E/<R»TF)) 

PRINT  105fT0,H,CP«,CP0,Z 
PRINT  106,TF,CV,P,-T 
PRINT  107,U0,TD,TL,EQ 
PRINT  108,01, 02, G,DG 

IF1XI.LT.RL.AND.XI1 .GT.RL)  GO  TO  18 
IF1XI.LTsRL.AND.XI1.LT.RL)  GO  TO  17 
PRINT  110 
60  TO  20 

FORMAT (IX, 15, IX, F9. 7) 

F0R»AT(1X,F4.4> 

PRINT  120 
GO  TO  20 

PRINT  130 


*  (Aittl  UV 

!Q  PRINT  140,XI,RL,TR,G,DP1 

GO  TO  2 

10  FORMAT! IX, 'BLOW  OUT'; 

20  FORMAT ( IX, 'PLANE  WITHOUT  IGNITION') 

30  FORMAT (IX, 'FLAME  REQUIRES  IGNITION') 

40  F0RHAT(1X,'XI*',F4.4,1X,'RL*',F6.4,1X,'TR»',F8.1,1X,'G=',F4.3,1X 

IFRAC  P  L0SS«',F5.4) 

113  FORMAT (1X,F5.3,1X,F5.3,1X,F9.7) 

105  FORMAT ( IX,  'TO*  ' ,F4.0 , 1 X, 'H*' ,E10.3,1 X,  'CPM*' ,E10.3,  < 

11X,'CP0*',E10.3,1X,'Z*',F4.3)  ,,, 

104  FORMAT! IX, 'TF»' ,FS.O, 1X,'CV*' ,E 10.3, 1X,'P»',E) 0.3) 

107  FORMAT ( 'UO*' ,F4. 1 , 1X,'TD*' ,F5.3,1 X, 'TL*',F4.3, IX, 'EQ*'F7.4 ) 

108  F0RHAT<1X,'D1»',F5.3,tX,'D2*',F5.3,1X,'G*',F4.3,1X,'D0»',F4.3> 

ISO  STOP 

END 

SUBROUTINE  CONV(X,Y,NR,NC) 

DIMENSION  XA(10)  ,YA( 10) 

IF ( ABS < X- Y ) .LT. .01 )  GO  TO  4 
IF(NC.LE.l)  GO  TO  5 

XT*(XA(NR)*Y-YA(NR )*X) /<XA(NR)-X*Y“YA(NR) ) 

XA<NR)«X 

YA<NR)«Y 

X*XT 

RETURN 

5  XA(NR)«X 

YA(NR)«Y 
X»Y 
NC*2 
RETURN 
4  X*Y 

NC«1 
RETURN 
END 

EOF: 

> 

EDIT: 

>0 

08.38.18  >lo$ 

10.28  ARU'S,  .21  CONNECT  HRS 
LOGGED  OFF  AT  08.38.27  ON  09AUG7? 


